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As  part  of  the  ongoing  program  in  "Decision  Control  Models  in  Operations 
Research,"  Mr.  J.  Christopher  Mitchell  has  studied  an  inventory  model 
that  uses  a  system-wide  stockout  constraint  rather  than  individual 
shortage  costs.  He  shows  how  to  easily  compute  policies.  The  computa¬ 
tional  experience  demonstrates  that  inventory  investment  can  be  reduced 
more  than  20%  using  this  approach.  The  report  also  suggests  how  to  impl 
ment  the  algorithm  in  a  real-life  system  that  may  contain  thousands  of 
stocked  items.  Other  related  reports  dealing  with  the  program  are  given 
on  the  following  pages. 
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ABSTRACT 


J.  CHRISTOPHER  MITCHELL.  Multi-Item  (s,S)  Inventory  Systems 
with  a  Service  Objective  (Under  the  direction  of  HARVEY  WAGNER.) 

This  paper  considers  a  multi-item  (s,S)  inventory  system. 

The  model  differs  from  standard  treatments  in  that  shortage  costs 
are  replaced  by  stockout  probability  constraints  to  be  satisfied 
in  every  period.  The  value  of  such  a  model  is  that  it  is  often 
easier  to  express  service  objectives  in  terms  of  stockout  prob¬ 
ability  constraints  than  it  is  to  specify  shortage  costs. 
Specifically,  system  service  is  defined  in  terms  of  a  weighted 
average  of  single-item  stockout  probabilities.  An  optimal 
policy  minimizes  system  cost  while  satisfying  a  constraint  on  system 
service.  Using  standard  single-item  approximations,  necessary 
and  sufficient  conditions  for  a  policy  to  be  optimal  (for  the 
approximate  model)  are  derived,  and  a  computationally  efficient 
algorithm,  the  Generalized  Knapsack  Duality  (GKD)  Algorithm,  is 
developed  to  find  such  a  policy.  Computational  experience  on 
inventory  systems  typical  of  many  found  in  the  real  world 
indicates  that  operating  costs  can  be  reduced  significantly  when 
this  model  is  used  rather  than  the  simpler  uniform  service  model 
often  used  by  managers. 

Sensitivity  experience  on  inventory  systems  with  a  struc¬ 
ture  typical  of  many  real-world  inventory  systems  is  reported. 


Specifically,  the  sensitivity  of  the  GKD  Algorithm  to  changes  in 
the  reorder  quantities  and  lower  bounds  on  feasible  policies  is 
reported.  This  experience  suggests  that,  as  in  the  single-item 
case,  very  accurate  specifications  of  the  reorder  quantities  are 
unnecessary.  Recommendations  are  made  to  specify  the  lower  bounds 
as  high  as  possible  while  achieving  significant  cost  savings  below 
that  of  the  Identical  Service  Approach.  This  recommendation  is 
consistent  with  the  objectives  of  many  managers,  and  in  certain 
cases  improves  algorithm  performance. 

This  paper  also  reports  computational  experience  with  sampling 
schemes  of  various  sizes  for  large-scale  inventory  aggregation. 

This  experience  indicates  that  inventory  systems  typical  of  many 
found  in  the  real  world  can  be  well-managed  based  on  decisions  made 
from  a  relatively  small  number  of  items  from  the  systems.  In 
particular,  this  experience  suggests  that  a  central  uniform  sample 
of  about  32  items  is  sufficient  to  make  accurate  decisions  for  the 
entire  inventory  system.  The  paper  concludes  with  a  detailed  im¬ 
plementation  procedure  using  the  GKD  Algorithm  with  sampling  to 
manage  a  large-scale  inventory  system. 
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I.  INTRODUCTION 


This  paper  examines  the  problem  of  specifying  single-item 
service  objectives  in  a  multi-item  inventory  system  subject  to  an 
overall  (or  system-wide)  service-level  constraint.  We  develop  a 
computationally  efficient  algorithm  for  doing  so  when  the  number  of 
items  in  the  inventory  system  is  of  moderate  size.  We  also  investi¬ 
gate  methods  of  aggregating  very  large  inventory  system  so  that  they 
are  computationally  more  manageable.  For  both  of  these  problems 
we  show  by  extensive  numerical  investigations  that  our  methods  can 
result  in  a  substantial  total  cost  savings  over  methods  which  specify 
uniform  service-level  objectives  for  all  items  in  the  inventory  system. 

This  chapter  is  a  non-technical  survey  of  the  literature  on  both 
exact  and  approximate  methods  for  inventory  management.  Section  1 
describes  the  theoretical  and  computational  difficulties  involved  in 
single-item  inventory  management  as  well  as  methods  that  have  been 
used  to  deal  with  these  difficulties.  Certainly  a  multi-item  in¬ 
ventory  system  will  inherit  these  difficulties,  and  so  we  use  such 
single-item  methods  in  our  model  when  appropriate. 

Section  2  discusses  service-level  constraints  in  single-item 
inventory  systems  in  the  same  spirit  as  Section  1.  We  also  survey 
the  methods  that  have  been  used  in  multi-item  inventory  systems  with 
service-level  constraints  (or  similar  models)  in  order  to  motivate 
the  direction  of  our  research. 
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Section  3  surveys  the  literature  on  the  problem  of  large-scale 
inventory  aggregation.  In  particular,  we  discuss  an  empirically- 
observed  structure  that  is  typical  for  many  real-world  inventory  systems. 

Chapter  II  contains  the  mathematical  derivation  of  our  algo¬ 
rithm  to  specify  single-item  service  objectives  in  multi-item  inventory 
systems  with  a  system-wide  service-level  constraint.  There  are  certain 
properties  the  items  in  the  inventory  system  must  satisfy  in  order 
to  guarantee  algorithm  convergence.  We  illustrate  algorithm  performance 
on  some  two-item  inventory  systems  which  contain  only  items  satisfying 
these  properties.  We  also  discuss  our  recommendations  when  there 
are  items  in  the  inventory  system  that  do  not  satisfy  all  these  proper¬ 
ties. 

Chapter  III  contains  a  numerical  investigation  of  some  32-item 
inventory  systems  (with  some  items  that  do  not  satisfy  the  above- 
mentioned  properties)  which  reflect  a  structure  often  observed  in 
practice.  We  show  that  there  is  a  significant  cost  decrease  when 
using  our  algorithm  to  specify  operating  policies  rather  than  using 
the  popular  method  of  specifying  operating  policies  which  give 
uniform  service. 

Chapter  IV  contains  a  numerical  sensitivity  investigation  of 
32- item  and  128- item  inventory  systems  (with  the  same  structure 
mentioned),  and  a  numerical  investigation  of  sampling  schemes  for 
inventory  aggregation. 

Chapter  V  concludes  this  paper  with  a  detailed  implementation 
procedure  of  the  GKD  Algorithm  with  sampling  to  manage  a  large- 
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scale  inventory  system.  We  illustrate  this  procedure  with  a  512- 
item  inventory  system,  and  offer  conclusions  and  directions  for  future 
research. 

1.  INVENTORY  THEORY  WITH  SHORTAGE  COSTS 

1 Model  Formulation  and  Optimal  Policies 

We  consider  the  periodic-review  dynamic  inventory  model.  We 
give  a  non-technical  survey  of  the  literature  on  this  model,  leaving 
relevant  mathematical  descriptions  for  the  next  chapter. 

At  the  beginning  of  each  period  n,  n=l,2,...,  the  inventory 
position  (stock  on  hand  plus  stock  on  order)  is  reviewed,  at  which 
time  a  positive  order  may  be  placed.  An  order  placed  in  period  n  is 
received  and  paid  for  in  period  n+k,  where  k,  the  leadtime,  is  a 
fixed  positive  integer.  There  is  charged  a  fixed  ordering  cost  K  plus 
a  linear  ordering  cost  c.  After  the  inventory  position  is  reviewed 
there  is  a  random  demand  E  ,  where  are  independent  and  iden¬ 

tically  distributed  with  cumulative  distribution  function  (cdf)  $>, 
density  <t>,  mean  y,  and  finite  standard  deviation  o.  After  demand 
realization,  a  linear  holding  cost  h  is  assessed  for  each  unit  of 
inventory  on  hand.  If  demand  exceeds  the  inventory  on  hand,  the  excess 
demand  is  completely  backlogged,  and  a  linear  shortage  cost  p  is 
assessed  for  each  unit  of  backlogged  demand.  Then  period  n+1  is  entered, 
repeating  the  process.  Future  costs  are  discounted  at  the  single¬ 
period  rategwith  0<3<1  (3=0  corresponds  to  a  single-period  model  and 
3=1  corresponds  to  an  undiscounted  model),  and  the  objective  is  to 
minimize  the  total  expected  cost  of  operating  the  system  over  a 
prescribed  horizon. 


This  model  can  be  formulated  as  a  dynamic  program.  Doing  so  and 
using  induction,  Scarf  [1960]  showed  that  under  the  assumptions  of 
finite  horizon  length  T  and  convex  differentiable  expected  single¬ 
period  holding  costs,  there  is  an  optimal  policy  of  the  form  (sn,Sn), 
n=l,...,T.  This  policy  requires  that  if  the  inventory  position  at 
the  beginning  of  period  n  is  less  than  sn,  an  order  is  placed  to  raise 
it  to  S  .  Otherwise,  no  order  is  placed.  Zabel  [1962]  extended  this 
result  to  the  case  when  the  single-period  holding  and  shortage  costs 
are  not  differentiable.  Earlier  Karlin  [1958a]  showed  the  optimality 
of  (sn,Sn)  policies  under  much  stronger  assumptions  than  Scarf;  and 
Veinott  [1966a]  and  Schal  [1976]  extended  Scarf's  result  to  models 
with  more  general  cost  functions.  If,  however,  holding  and  shortage 
costs  are  linear,  as  in  our  model.  Scarf's  result  is  sufficient  to 
guarantee  the  optimality  of  (sn>Sn)  policies. 

A  basic  difficulty  in  implementing  these  policies  is  that  they 
usually  vary  from  period  to  period  and  hence  require  tremendous 
computational  effort.  The  situation  becomes  much  simpler  in  the 
infinite  horizon  case.  In  this  case,  to  minimize  total  expected  cost 
we  must  impose  6<1  or  else  for  almost  any  policy  this  cost  is  infinite 
One  can  consider  an  undiscounted  model  (8=1)  if  the  objective  criterion 
is  changed  to  the  average  expected  cost  per  period.  If  this  is  the 
criterion,  Iglehart  [1963a]  and  [1963b]  proved  that  given  the  cost 
structure  we  have  assumed,  in  both  the  discounted  and  undiscounted  case 
respecti vely,  there  is  an  optimal  policy  that  is  stationary  (s,S). 

Such  a  policy  requires  that  when  the  inventory  position  in  any  period 
falls  below  s,  an  order  is  placed  to  bring  it  up  to  S.  Thus,  there  are 
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only  two  numbers  to  be  computed.  One  must,  of  course,  assume  that  the 
true  horizon  length  is  of  sufficient  length  to  be  reasonably  approxim¬ 
ated  by  an  infinite  horizon  model. 

The  computational  procedures  to  find  exact  stationary  (s,S)  policies 
generally  involve  deriving  a  closed  form  expression  for  the  expected 
infinite  horizon  average  cost  per  period  given  that  the  policy  (s,S)  is 
being  followed  in  every  period,  and  then  minimizing  this  cost  over  all 
(s,S)  policies.  There  have  been  several  methods  used  to  derive  this 
steady-state  cost.  Karlin  [1958a]  derived  it  using  linear  operator 
theory  and  in  [1958b]  using  renewal  theory,  which  was  simplified  by 
Sahin  [1982].  We  mention  that  in  both  references  cited,  Karlin  also 
derived  a  closed  form  expression  for  the  optimal  policy  when  the  demand 
distribution  is  exponential.  Morse  [1959]  derived  the  expected  average 
cost  per  period  using  Markov  processes  theory,  and  Leneman  and  Beutler 
[1969]  used  a  stationary  point  process  approach.  The  first  general  and 
exact  computational  procedure  to  find  an  optimal  stationary  (s,S)  policy 
is  found  in  Veinott  and  Wagner  [1965].  Assuming  that  the  demand  distrib¬ 
ution  is  discrete,  they  derived  the  steady-state  cost  and  then  minimized 
it  using  finite  difference  calculus.  This  algorithm  has  been  programmed 
in  PL-1,  and  the  documentation  can  be  found  in  Kaufman  [1976],  Bell 
[1970]  has  suggested  an  improvement  to  the  Veinott-Wagner  algorithm 
using  optimal  stopping  rule  theory,  and  recently  Federgruen  and  Zipkin 
[1981]  have  described  a  completely  different  approach  using  a  policy 
iteration  technique  for  Markov  decision  processes.  A  method  similar  to 
this  was  suggested  earlier  by  Johnson  [1968].  In  a  somewhat  different 
spirit,  Sivazlian  [1971]  uses  Gaussian  quadrature  methods  to  create 
graphs  to  compute  optimal  (s,S)  policies. 
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There  are  three  significant  difficulties  involved  in  finding 
optimal  (s,S)  policies.  We  discuss  these  in  order  to  introduce  the 
various  approximations  to  optimal  (s,S)  policies  that  have  been 
suggested  in  the  literature. 

First,  the  steady-state  costs  always  involve  the  renewal  function 
of  the  demand  process.  If  the  demand  distribution  is  discrete,  Veinott 
and  Wagner  [1965]  showed  that  the  renewal  function  can  be  evaluated 
recursively,  although  this  can  be  computationally  expensive  when  S  is 
very  large.  If  the  demand  distribution  is  continuous,  there  is  no 
general  computationally  efficient  method  of  even  approximating  the 
renewal  function  (the  exception  being  the  exponential  distribution). 

Second,  although  the  steady-state  cost  is,  in  general,  convex  in 
S,  it  is  not  even  unimodal  in  D=S-s.  Thus,  local  minima  may  not  be 
global  minima.  Veinott  and  Wagner  [1965]  dealt  with  this  by  establish¬ 
ing  bounds  on  the  optimal  (s,S)  policy  and  examining  all  values  of  D 
among  these  policies  to  find  a  global  minimum.  Recently  Sahin  [1982] 
has  shown  that  for  a  class  of  distributions  the  steady-state  cost  is 
pseudo-convex  in  D  and  S,  so  local  minima  are  global  minima.  It  is  not 
clear,  however,  that  this  class  contains  distributions  useful  for  inven¬ 
tory  model s . 

Third,  exact  algorithms  require  that  the  demand  distribution  be  fully 
known.  Many  times  in  practice  only  a  couple  of  moments  will  be  known, 
and  even  these  values  may  be  statistical  estimates. 


1.2  Approximations  to  Optimal  Policies 

In  light  of  the  difficulties  involved  in  finding  optimal  policies, 
numerous  approximation  procedures  have  been  suggested.  We  describe  the 
best  of  these. 

For  an  undiscounted  infinite  horizon  inventory  model  the  ordering 
cost  c  does  not  appear  in  the  cost  functions.  Essentially  this  is 
because  all  demand  must  eventually  be  satisfied,  and  the  cost  of  doing 
so  is  undiscounted.  Thus,  there  are  three  components  of  total  expected 
cost,  one  associated  with  each  of  the  parameters  K,  h  and  p.  They  can 
be  denoted,  respectively,  as  the  expected  replenishment  cost,  expected 
holding  cost  and  expected  shortage  cost.  Another  quantity  of  interest 
is  the  service-level,  the  frequency  of  periods  without  any  backorder 
(roughly,  the  steady-state  probability  of  meeting  all  demand  in  any 
given  period) . 

One  of  the  earliest  and  best  approximations  is  the  Normal  Approxi¬ 
mation  of  Wagner  [1975,  pp. 831-836].  It  is  based  on  asymptotic  renewal 
theory  results  of  Roberts  [1962],  the  empirically-based  heuristic  when 
D=S-s  is  small  discussed  in  Wagner  et  al_  [1965],  and  the  assumption 
that  demand  is  well-enough  approximated  by  a  Normal  distribution. 

There  were  extensive  numerical  investigations  of  this  approximation 
done  in  MacCormick  [1974],  Estey  and  Kaufman  [1975],  MacCormick  et  al_ 
[1977],  and  MacCormick  [1977]  in  which  the  actual  underlying  demand 
distribution  was  assumed  to  be  Negative  Binomial  or  Poisson.  The  Normal 
Approximation  in  general  performed  very  well  even  when  statistical 
estimates  of  the  demand  distribution  mean  and  variance  were  used 
(although  degradation  was  greater  in  this  case).  The  greatest  degrad- 


8 


ation  in  performance  was  observed  when  the  coefficient  of  variation 

-  was  large  (in  which  case  there  is  a  non-negligible  probability  of 

2 

negative  demand  for  a  Normal  random  variable),  and  when  both  — p  and 

^  were  large  (the  latter  corresponds  to  a  high  service  level).  When 
2 

and  jj-  were  large,  the  expected  shortage  cost  and  service-level 

degraded  much  more  than  the  expected  replenishment  and  holding  costs. 

Using  the  general  approach  of  Norman  and  White  [1968]  of  computing 
approximately  optimal  policies  in  Markov  Decision  processes  by  replacing 
the  probability  distributions  by  their  moments,  Porteus  [1979]  developed 
an  algorithm  to  compute  approximately  optimal  (s,S)  policies.  However, 
the  empirical  examination  done  in  Freeland  and  Porteus  [1980]  showed 
that  it  was  not  much  better  than  the  Normal  Approximation,  and  was  a 
lot  more  work  to  compute. 

Using  regression  models  suggested  by  both  asymptotic  analysis  and 
empirical  observations,  Ehrhardt  [1976]  developed  the  Power  Approximation. 
Extensive  numerical  investigations  were  performed  by  Ehrhardt  [1976] 
(compare  Ehrhardt  [1978])  and  Klincewicz  [1976a]  and  [1976b].  These 
investigations  showed  that  not  only  are  the  approximations  especially  easy 
to  compute,  require  only  the  first  two  moments  of  the  demand  distribution 
and  very  useful  for  sensitivity  analysis,  but  they  are  extremely  accurate 
over  a  wide  range  of  parameter  settings.  The  Power  Method's  performance 
was  always  superior  to  the  performance  of  the  Normal  Approximation.  The 
only  cases  when  the  Power  Approximation  did  not  perform  well  occurred  in 
the  expected  shortage  cost  and  service-level  when  statistical  demand 

moments  were  used  and  jj-  was  large  (which  implies  a  high  service-level). 
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Some  technical  improvements  were  made  in  the  Power  Approximation  by 
Hosier  [1981],  and  the  regression  approach  was  used  by  Ehrhardt  [1977] 
(compare  Ehrhardt  [1981])  to  approximate  the  expected  costs  and  other 
quantities  of  interest  when  the  (s,S)  policy  is  given. 

There  are  other  policies  of  a  suboptimal  form  that  have  been 
used  to  control  inventory  systems,  primarily  because  these  policies 
are  considerably  easier  to  compute  and  analyze  than  are  (s,S)  policies. 
Two  of  the  most  popular  are  the  (t,S)  and  (s,Q)  policies  (see  Hadley 
and  Whitin  [1963,  pp. 235-295]) .  A  (t,S)  policy  requires  that  every  t 
periods  an  order  be  placed  to  raise  the  inventory  position  to  S.  An 
(s,Q)  policy  requires  that  whenever  the  inventory  position  falls  below 
s,  an  order  of  size  Q  is  to  be  placed.  Using  these  policies  Naddor 
[1975]  developed  rules  for  approximating  the  optimal  (s,S)  policy. 

These  approximations  were  numerically  compared  in  Kastner  and  Ehrhardt 
[1979]  and  in  Ehrhardt  and  Kastner  [1980].  The  Power  Approximation  was 
always  superior,  although  the  Naddor  Approximation  did  quite  well  for  a 
large  number  of  parameter  settings. 

We  conclude  by  summarizing  how  these  approximately  optimal  policies 

deal  with  the  three  difficulties  described  in  Section  1.1.  Typically, 

the  renewal  function  is  approximated  using  standard  results  from  asymp¬ 
totic  renewal  theory.  Although  we  did  not  discuss  it  in  detail  in 

this  section,  the  complicated  total  expected  cost  dependence  on  0  is 
also  dealt  with  using  asymptotic  renewal  theory,  with  an  empirically- 
based  heuristic  being  implemented  when  D  is  small  (in  other  words,  the 
asymptotic  assumptions  are  unreasonable).  In  order  to  have  approxi¬ 
mations  that  depend  upon  only  a  few  moments  of  the  demand  distribution 
rather  than  the  entire  density,  the  form  of  the  demand  distribution  is 
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assumed  to  be  given  (typically  Normal,  Gamma,  Poisson  or  Negative 
Binomial).  The  Power  Approximation  utilized  a  very  general  distribution 
form  in  which  parameters  were  determined  by  a  regression  fit. 

2.  INVENTORY  THEORY  WITH  SERVICE-LEVEL  CONSTRAINTS 
2.1  Model  Formulation  and  Optimal  Policies 

A  significant  practical  difficulty  inherent  in  all  these  approxi¬ 
mations  (and  exact  algorithms)  is  the  specification  of  shortage  costs. 
These  penalize  backorders,  but  it  can  be  very  difficult  to  measure  the 
cost  impact  of  a  backorder.  It  entails  such  things  as  "loss  of  goodwill" 
or  "customer  dissatisfaction",  or  equally  difficult  quantities  to 
measure.  One  way  to  deal  with  this  difficulty  is  to  replace  the  short¬ 
age  cost  with  a  service-level  constraint.  For  example,  one  can  require 
that  the  frequency  of  backorders  be  no  greater. than  15%.  Such  measures 
can  often  be  easier  to  specify  than  shortage  costs. 

A  question  of  some  interest  is  whether  or  not  there  is  an  optimal 
policy  of  the  (s,S)  form  when  a  service-level  constraint  is  used  rather 
than  a  shortage  cost.  In  general,  there  is  not.  Using  a  fairly  general 
service-level  constraint,  Fromovitz  [1965]  showed  that  even  in  a  single¬ 
period  model  the  optimal  policy  may  be  randomized.  He  does  show  that 
if  the  service-level  constraint  is  convex  there  is  an  optimal  policy 
that  is  non-randomized,  but  most  service-level  constraints  considered 
in  the  literature  are  not  convex.  Beesack  [1967]  has  shown  for  the 
finite-horizon  model  (with  the  constraint  that  the  expected  number  of 
stockouts  be  almost  a  prescribed  fraction  of  expected  demand)  that 
(sn,Sn)  policies  are  still  optimal.  However,  as  most  managers  are  un¬ 
willing  to  implement  a  policy  more  complicated  than  an  (s,S)  policy. 
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most  of  the  research  is  aimed  at  finding  the  best  (s,S)  policy  that 
meets  a  given  service -1 evel . 

2.2  Single-Item  Models 

Roberts  [1962]  showed  that  asymptotically  the  optimal  value  for 
D=S-s  is  independent  of  the  shortage  cost  p,  so  typically  it  is  assumed 
that  D  is  given,  and  then  approximations  for  s  are  derived  such  that 
the  resulting  (s,S)  policy  satisfies  a  prescribed  service-level.  The 
two  servi ce -level s  most  often  considered  are  the  one  we  use  (frequency 
of  backorders)  and  the  one  Beesack  considered  (fraction  of  expected 
demand  met).  In  the  spirit  of  Robert's  work,  Greenberg  [1964]  and 
Schneider  [1978]  derived  closed  form  expressions  for  these  two  service- 
levels  as  functions  of  s  and  D.  Using  asymptotic  renewal  theory  and 
both  a  Normal  and  Gamma  demand  distribution,  Schneider  [1978]  derived 
approximations  to  find  an  approximately  optimal  s  given  D  for  both 
service -level  constraints.  Only  the  approximation  to  the  constraint 
used  by  Beesack  was  investigated  numerically.  The  Gamma  Approximation 
always  performed  well,  and  the  Normal  Approximation  performed  well  as 
long  as  the  coefficient  of  variation  was  not  too  large.  Tijms  and 
Groenevelt  [1982]  extended  these  results  to  the  case  when  the  leadtime 
is  a  random  variable. 

2 . 3  Multi -I tern  Models 

We  now  discuss  the  model  whose  study  is  the  object  of  this  paper. 
We  consider  an  N-item  inventory  system  in  which  each  item  i  is  an 
inventory  model  described  in  Section  1.1  without  shortage  costs 
assigned.  Thus,  for  each  item  i,i=l,...,N,  we  have  a  fixed  leadtime 
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k,j,  fixed  ordering  cost  plus  linear  ordering  cost  c^ ,  and  linear 

holding  cost  .  The  demand  realized  for  item  i  in  any  period  has 

continuous  cdf  4^ ,  density  ,  mean  and  finite  standard  deviation 

o.j •  All  demands  for  all  items  in  all  periods  are  assumed  to  be 

independent.  The  criterion  of  optimality  is  the  minimum  expected 

undiscounted  cost  per  period  over  an  infinite  horizon,  and  the  service- 

level  constraint  to  be  satisfied  in  every  period  is 
N 

T:  [frequency  of  periods  item  i  is  not  backordered]  £  a, 

i=l  N 

where  >  0  and  0  <  a  <  1  are  specified  with  ^  W^=l . 

i=l 

A  feasible  policy  is  a  pair  (s,S),  where  s=(s.| . .  ,sN)  and 

Ss(Sj . ,SN),  and  the  stationary  policy  (s^.S..)  is  followed  for  item 

i  in  every  period. 

This  model  has  not  been  examined  in  the  literature,  but  several 
related  (but  much  simpler)  models  have.  We  briefly  survey  the  results. 

For  a  single-item  inventory  model  with  a  shortage  cost,  if  K=0 
there  is  an  optimal  policy  of  the  form  (S,S),  a  base-stock  policy  (see 
Veinott  [1966b]).  Iglehart  and  Jaquette  [1969]  showed  that  if  the 
shortage  cost  is  replaced  by  a  service-level  constraint,  a  base-stock 
policy  is  still  optimal.  Mitchell  [1982]  considered  this  model  with 
all  leadtimes  k^=0  and  exponential  demands,  and  developed  a  computa¬ 
tionally  efficient  algorithm  to  find  the  optimal  base-stock  policy. 

Evans  [1967]  considered  a  finite-horizon  base-stock  model  in  which  total 
costs  were  minimized  subject  to  a  linear  resource  constraint.  Using 
induction  on  the  dynamic  programming  formulation,  he  derived  the 
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(extremely  complicated)  optimal  policy. 

In  Section  1.2  we  described  suboptimal  policy  forms  which,  because 
of  the  relative  ease  in  analyzing  them,  have  been  used  in  inventory 
control  when  K^>0.  Recall  we  described  the  (s,Q)  policy  (often  called 
a  (Q,r)  policy  in  the  literature),  which  requires  that  whenever  the 
inventory  position  drops  below  s,  an  order  of  size  Q  is  placed.  There 
are  several  papers  that  consider  our  multi-item  model  (with  perhaps  a 
different  constraint)  governed  by  stationary  (s,Q)  rather  than  (s,S) 
pol icies. 

There  are  two  approaches  to  solving  such  problems  that  have  been 
used,  both  of  which  are  discussed  in  Hadley  and  Whitin  [1963,  pp . 21 3- 
219,  304-307,  323-226].  One  is  to  treat  the  problem  as  a  constrained 
nonlinear  program  and  then  solve  the  first-order  Kuhn-Tucker  conditions. 
The  second  is  to  formulate  the  model  as  a  dynamic  program  (a  generalized 
Knapsack  Problem)  and  then  solve  it  stage  by  stage. 

The  method  using  the  Kuhn-Tucker  conditions  is  used  in  Winters 
[1962],  Parker  [1964],  Gerson  and  Brown  [1970],  Presutti  and  Trepp 
[1970],  Schrady  and  Choe  [1971],  and  Schroeder  [1974].  The  method 
using  dynamic  programming  is  used  in  Kaplan  [1978]. 

Winters  [1962]  minimized  total  expected  cost  subject  to  a  product¬ 
ion  smoothing  constraint.  Parker  [1964]  minimized  total  expected  cost 
subject  to  a  service-level  constraint.  Gerson  and  Brown  [1970]  mini¬ 
mized  total  expected  cost  subject  to  a  service-level  constraint  and 
s=0.  Presutti  and  Trepp  [1970]  also  minimized  total  expected  cost 
subject  to  a  service-level  constraint.  Schrady  and  Choe  [1971]  maxi¬ 
mized  a  service-level  criterion  subject  to  investment  budget  constraints. 
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as  did  Schroeder  [1974].  Kaplan  [1972]  considered  a  somewhat  different 
model,  a  finite-horizon  model  with  budget  constraints  that  penalize 
dollars  spent  before  the  end  of  the  horizon.  It  is  a  single-item 
model  with  s=0.  We  note  that  the  models  in  Parker  [1964]  and  Schroeder 
[1974]  assume  a  continuous  rather  than  periodic  review  of  the  inventory 
position.  Such  models  tend  to  produce  policies  that  give  rather  poor 

performance  in  a  periodic-review  setting  (see  Wagner  et  al_  [1965]). 

For  two  reasons,  all  of  these  studies  are  somewhat  unsatisfactory. 

First,  as  mentioned,  (s,Q)  policies  are  not  optimal.  Second,  in  none 
of  them  are  the  policies  produced  compared  with  any  other  method,  making 
it  difficult  to  evaluate  the  performance  of  such  policies  in  a  multi-item 
inventory  system.  In  our  paper  we  use  (s,S)  policies  and  derive  an  algo¬ 
rithm,  the  GKD  (Generalized  Knapsack  Duality)  Algorithm  to  approximately 
solve  our  model,  and  then  compare  the  (s,S)  policies  produced  with  an 
approach  frequently  used  in  practice  by  managers,  the  Identical  Service 

Approach  (compare  Mitchell  [1982]).  This  approach  sets  shortage  costs 
for  each  item  in  an  inventory  system  so  that  the  service-level  for 

each  item  is  the  same,  say  a.  Then  each  item  is  treated  individually, 

with  policies  being  computed  by  some  existing  method.  We  will  show 

that  using  our  model,  which  in  effect  varies  the  individual  item  service- 

levels  while  still  maintaining  an  overall  service-level  of  a  (via  our 

service-level  constraint),  total  expected  costs  can  be  significantly 

reduced  for  many  multi-item  inventory  systems.  This  is  the  first  issue 

wi th  which  we  deal . 

The  second  issue  involves  aggregating  very  large  inventory  systems 
in  such  a  way  that  a  representative  fraction  of  them  can  be  used  to 
determine  the  policies  for  all  of  them,  which  we  discuss  next. 
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3.  MULTI-ITEM  INVENTORY  SYSTEM  AGGREGATION 

Many  real-world  inventory  systems  contain  thousands  of  items. 

It  would  be  cost  prohibitive  to  try  to  include  every  one  of  them  in 
any  model  which  is  to  be  used  to  set  inventory  policies.  Inventory 
aggregation  is  any  method  of  partitioning  such  a  system  into  blocks, 
each  of  which  will  be  dealt  with  essentially  the  same.  For  example, 
a  representative  item  could  be  chosen  from  each  block  to  determine 
the  policy  that  will  be  used  for  all  items  in  the  block.  Thus  the 
problem  is  reduced  to  more  manageable  proportions.  The  question  is, 
of  course,  how  to  aggregate  so  that  the  resulting  policies  are  as 
close  as  possible  (in  a  total  expected  cost  sense)  to  policies  that 
would  have  been  produced  had  every  item  in  the  inventory  system  been 
included  in  the  model . 

Very  little  is  known  about  this  problem  (see  Bitran  and  Hax  [1982], 
which  deals  with  a  finite  horizon,  deterministic  demand  inventory  system). 
However,  extensive  empirical  investigations  of  many  real-world  inventory 
systems  have  suggested  some  remarkable  similarities  in  the  overall  dis¬ 
tribution  of  items  in  inventory  systems.  These  results  are  discussed  in 
Brown  [1959],  [1963],  [1967]  and  [1977],  and  Peterson  and  Silver  [1979, 
pp.  30-37,  71-80].  We  summarize  the  observations  in  Peterson  and  Silver. 

The  annual  dollar  demand  of  an  item  is  defined  to  be  its  annual 
expected  demand  times  its  value  per  unit.  Of  course,  the  "value"  of 
an  item  may  be  difficult  to  determine  exactly,  but  can  often  be  taken 
to  be  its  ordering  or  holding  cost.  Typically,  five  to  ten  percent  of 
the  items  in  an  inventory  system  account  for  about  fifty  percent  of 
the  total  annual  dollar  demand  of  the  entire  inventory. 


These  items  tend  to  be  very  important,  and  are  denoted  Type  A 
items.  On  the  other  hand,  typically  about  thirty  to  fifty  percent 


16 


of  the  items  account  for  only  five  to  ten  percent  of  the  total 
annual  dollar  demand.  These  are  the  relatively  unimportant  items, 
and  are  denoted  Type  C  items.  The  rest  of  the  items  in  the  inventory 
system,  the  moderately  important  ones,  are  denoted  Type  B  items. 

If  the  items  in  such  an  inventory  system  are  ordered  in  descend¬ 
ing  order  by  their  annual  dollar  demand,  and  then  the  cumulative 
percent  of  the  items  is  plotted  against  the  corresponding  cumulative 
percent  of  annual  dollar  demand,  the  graph  would  look  something  like 
Exhibit  1.  Notice  also  that  about  twenty  percent  of  the  items  account 
for  about  eighty  percent  of  the  total  annual  dollar  demand.  This  also 
is  typical  for  many  real-world  inventory  systems. 


£ 


s 

__  V. 
Type  A 


Cumulative  Percent  of  Items 

Exhibit  1 
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Empirical  observations  have  shown  that  this  graph  can  very  often 
be  well -approximated  by  a  lognormal  cumulative  distribution  function. 

We  address  the  second  issue  of  this  paper,  large-scale  inventory 
aggregation,  by  showing  that  sampling  techniques  which  stratify  on  the 
item  types  and  which  maintain  the  structure  illustrated  in  Figure  1  can 
be  used  to  accurately  predict  the  (s,S)  policies  for  all  the  items  in 
the  inventory  system.  We  exploit  the  special  structure  of  our  inventory 
problem  by  using  the  sample  to  estimate  the  Lagrange  constraint  multi¬ 
plier  of  the  inventory  problem,  from  which  we  compute  all  (s,S)  policies. 


II.  THE  GENERALIZED  KNAPSACK  DUALITY  (GKD)  ALGORITHM 


1.  INTRODUCTION 

We  consider  a  multi-item  inventory  system  and  study  the  problem 
of  specifying  both  system-wide  and  individual  service  objectives. 

Most  theoretical  models  considered  in  the  literature  assume  that  hold¬ 
ing  or  shortage  costs  are  applied  to  any  excess  inventory  or  unsatisfied 
demand,  respectively.  A  major  difficulty  in  applying  such  models  in 
practice  is  the  specification  of  shortage  costs.  Frequently,  the 
manager  sets  shortage  cost  parameters  based  on  an  objective  of  satis¬ 
fying  demand  with  at  least  some  minimum  probability.  Such  an  approach 
may  be  preferred  because  subjective  factors  can  more  easily  be  expressed 
as  a  probability  of  satisfying  demand  than  as  a  cost  for  each  shortage 
incurred. 

A  significant  shortcoming  of  this  approach  is  that  it  usually 
entails  the  setting  of  a  probability  of  demand  satisfaction  that  is 
applied  uniformly  to  all  items  in  the  system.  In  this  paper  we  raise 
the  issue  of  specifying  different  service  objectives  for  individual 
items  while  still  satisfying  some  given  system-wide  objective.  The 
value  of  this  approach  is  that  total  system  costs  can  be  reduced  below 
those  of  the  method  that  requires  identical  service  for  all  items. 

We  are  interested  in  the  following  multi-item  inventory  system. 

For  each  item  i  there  is  a  set-up  ordering  cost  K,. ,  a  unit  ordering 
cost  c.j ,  and  a  fixed  leadtime  k,.  between  placement  and  delivery  of 
orders.  These  costs  are  assessed  upon  order  delivery.  There  is  a  unit 
holding  cost  h^  assessed  at  the  end  of  each  period,  and  the  demand 
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realized  for  for  item  i  in  any  period  has  absolutely  continuous  cumula¬ 
tive  distribution  function  (cdf)  £. ,  density  4>^ ,  mean  and  finite 
standard  deviation  All  demands  for  all  items  in  all  periods  are 
mutually  independent.  There  are  N  items,  and  the  criterion  of 
optimality  is  the  minimum  expected  undiscounted  cost  per  period  over 
an  infinite  horizon,  we  do  not  include  shortage  costs  to  penalize 
back-orders.  Rather  we  introduce  the  following  service-level  constraint 
to  be  satisfied  in  every  period: 


N 


W. [frequency  of  periods  item  i  is  not  backordered]  >  at, 


i  =  l 


N 


where 


For  such  an  inventory  model,  an  (s,S)  policy  requires  that  whenever 
the  inventory  position  (inventory  on  hand  plus  on  order)  drops  below 
s,  an  order  is  placed  bringing  the  inventory  position  up  to  S.  An 
(s,S)  policy  is  called  stationary  if  it  is  applied  in  every  period.  We 
only  consider  stationary  (s,S)  policies  because  when  shortage  costs  are 
included  in  the  model,  there  is  an  optimal  policy  which  is  stationary 
and  (s,S)  (Iglehart  [1963]).  Any  stationary  (s.S)  policy  for  our  model 
gives  rise  to  a  unique  and  well-defined  specification  of  shortage  costs. 
If  we  denote  by  P.(s.,S. )  the  frequency  of  periods  that  item  i  is  not 
backordered  when  following  the  policy  (s^,S^),  then,  when  is  contin¬ 
uous  , 


(1) 


where  p.  is  the  shortage  cost  under  which  the  policy  (s-,S-)  is  optimal 
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(Veinott  and  Wagner  [1965]).  We  can  solve  this  equation  for  p.  and  thus 
determine  the  shortage  costs  associated  with  any  stationary  (s,S)  policy. 

For  our  model  a  feasible  policy  is  a  set  |(s,S)}  ,  where 
s=(s^ . . .s^)  and  S=(Sj ,.. . ,$N).  For  item  i  we  follow  the  policy 

(s^,S^)  in  every  period.  We  define  D-j  by  D^=S^-s^,  and  denote  by 

the  (k^+l)-fold  convolution  of  I..  with  itself,  and  by  its  density. 

Also  let  ii.*  and  o*  be  the  mean  and  standard  deviation  of  a  random 
n  i 

variable  with  cdf  In  Mitchell  [1982]  we  considered  this  problem 

when  K^=0  for  all  items  i,  and  developed  a  computationally  efficient 
algorithm  for  the  case  when  all  k^=0  and  all  demands  are  exponentially 
distributed.  We  extend  those  results  in  this  paper,  dropping  all  three 
of  these  restrictions. 

Section  2  describes  the  formulation  of  our  inventory  problem  as 
a  non  linear  program.  Using  the  Kuhn-Tucker  conditions,  we  prove 
necessary  and  sufficient  conditions  for  (possibly  local)  optimality. 

The  main  results  are  found  in  Theorems  1,  2,  6,  7  and  8  (the  reader 
may  choose  to  skip  the  proofs  of  these  results).  There  is  at  the  end 
of  this  chapter  a  list  of  the  notation,  a  summary  of  frequently  used 
relations,  and  a  summary  of  the  theorems  and  lemmas. 

Section  3  is  a  description  of  the  algorithm.  It  is  based  on 
local  duality  theory  of  nonlinear  programming,  and  the  proof  of  con¬ 
vergence  uses  results  from  Section  2.  Section  4  reports  a  numerical 
investigation.  In  Section  5  we  discuss  our  recommendations  when  the 
sufficient  conditions  may  not  hold  for  some  items  in  the  inventory 
system.  Finally,  Section  6  contains  some  conclusions. 
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2.  ANALYSIS 

If  we  denote  by  C^(S^.D^)  the  expected  average  undiscounted  cost 
per  period  over  an  infinite  horizon,  and  by  P^(S^.D^)  the  frequency  of 
periods  that  item  i  is  not  backordered,  we  have  the  following  result. 


Theorem  1 .  For  S^>D^  , 


1 


(i)  C.(S.,D.)  =  - 

1+M.(D.) 


h. 


Di  sry 

//  (Si-y-x)<t>i*(x)mi  (y)dxdy 


o 


A ^Si-x)$i*(x)dx 


4  K. 


(ii)  P.(S.,D.)  =  - - -  /*  5>(sry)m.(y)dy +I.*(S.)  , 

1+W  o 

where  M.  is  the  renewal  function  of  3Lj  and 

Proof.  See  Roberts  [1962],  Greenberg  [1964]  and  Ehrhardt  [1981]  □. 

For  two  reasons  we  restrict  attention  to  (s^S^)  policies  with 
s^>0  (equivalently,  S^>D^).  First,  the  expressions  given  in  Theorem  1 
do  not  hold  when  s^<0.  One  must  replace  with  ( 0 , S - )  the  limits  of 
integration  in  the  first  integral  of  and  the  limits  of  integration 
in  P...  Second,  many  managers  consider  policies  with  s^<0  to  be  unsat¬ 
isfactory  for  practical  implementation. 

Thus  we  want  to  solve  the  following  constrained  nonlinear  program 
(NLP)  to  obtain  an  optimal  policy  for  our  inventory  system: 
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minimize  C(S,D)  =  &‘VV 

N 

subject  to  P(S,D)  =  7:  Ws  [p.(S.,D.)-aj^0 


i  =  l 

,  i  =  l .  ,N  , 


where  0<a<l ,  W^,...,W^>0  and  ^  W^l . 

i  =  l 

There  are  three  serious  difficulties  in  finding  a  solution.  First, 
the  renewal  function  M.  generally  does  not  have  a  closed-form  expression, 
and  usually  is  complicated  when  it  does.  Second,  the  functional  depend¬ 
ence  of  and  P^  on  the  variable  is  complicated,  as  is  easily  seen 
if  one  works  out  the  partial  derivatives.  Third,  the  behavior  of  the 
functions  depends  heavily  on  the  demand  distributions  . 

We  deal  with  the  first  difficulty  by  using  asymptotic  approximations 


for  the  renewal  functions.  In  particular.  Smith  [1954]  has  shown 

(2) 


y  o2.  1 

Mj(y)  =  zr  +  — -  2  +  °(])  as  y—. 

*>  2p. 


We  use  the  following  approximations. 

2 


w  ^  + 


°i  1 

,2'2' 

2ui 


Hi 


HM.(Di)  D.+ty+o^/u.)/^ 


5  Pi 


(3) 


We  approximate  m.(y)=M/(y)  by  a  constant  function,  since  M,.(y) 
is  asymptotically  affine.  One  could  differentiate  (2)  and  approxi- 
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mate  m^(y)  by  1,  but  we  use  the  somewhat  better  approximation  discussed 
in  Ehrhardt  [1981],  namely. 


l-pi 

mi(y)  *  pTUT’  ye(*1.Si). 


(4) 


(This  approximation  approaches  -  as  D 

V  1 

For  the  complicated  functional  dependence  on  ,  we  approximate 
D.j  by  a  constant,  and  regard  and  as  functions  only  of  . 
Roberts  [1962]  proves  that  in  the  single-item  case  where  a  shortage 
cost  p  is  specified,  the  optimal  value  D*  of  D  is  given  by 


+  constant  + 


o(l) 


1/ 

as  K-~c,  p-«r  sUch  that  -  remains  constant. 

Since  D*  is  asymptotically  independent  of  the  shortage  cost 
p,  and  hence  by  (1)  of  the  service-level  a,  it  would  seem  reasonable 


to  approximate  D*  by  the  constant  the  so  called  "economic  order 
quantity." 

An  extensive  empirical  comparison  of  approximately  optimal 
inventory  policies  by  Wagner,  et  al  [1965],  showed  this  to  be  a 
remarkably  good  approximation,  even  for  K  and  p  fairly  small. 


although  it  seems  to  give  better  results  when 


2Ku 


1  •  * 


Moreover, 


the  total  expected  cost  is  typically  quite  flat  near  D*.  so  if  S*  is 
well-approximated,  the  approximation  for  D*  does  not  need  to  be 
especially  accurate  to  have  a  total  cost  close  to  the  cost  associated 
with  the  optimal  policy. 
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Actually,  we  use  the  Power  Approximation  discussed  in  Mosier 
[1981]  for  the  optimal  D.  He  reports  empirical  investigations  that 
indicate  that  this  approximation  is  superior  to  the  economic  order 
quantity.  As  in  the  case  of  the  economic  order  quantity,  this  approx¬ 
imation  tends  to  be  poor  when  it  is  smaller  than  1.5^.  For  the 
derivation  of  our  algorithm,  we  assume  is  at  least  ,  and  thus 
use  the  approximation 


* 

D.j  =  max  • 

,  ,  .494  ( ^ 

’  1  '3yi  \h 

r 

;-Ci 

r] 

.116 

rv  • 

We  note  that  asymptotically  (as  -  °°)  this  approximation  differs 
very  little  from  the  economic  order  quantity. 

To  deal  with  the  third  difficulty,  we  intended  to  assume  that 
the  demand  distribution  in  Normal,  since  Wagner,  et  a^  [1965],  Wagner 
[1975,  ch.  19],  MacCormick  [1974],  Estey  and  Kauffman  [1975],  MacCormick 
et  aj.  [1977],  and  MacCormick  [1977]  all  suggest  that  a  Normal  distribution 
gives  good  approximations  to  the  optimal  policies  even  when  the  true 
distribution  is  skewed  (By  a  good  approximation  we  mean  that  the 
expected  total  cost  of  the  approximation  policy  is  nearly  the  same  as 
for  an  optimal  policy).  But  because  the  Normal  cdf  cannot  be  expressed 
in  a  closed  form,  we  use  a  Logistic  distribution  to  approximate  the  Normal 
distribution.  It  is  an  especially  simple  and  accurate  approximation  to 
the  Normal  distribution,  as  discussed  in  Johnson  and  Kotz  [1970,  ch.  22]. 


They  show  that  if  ?*  is  the  cdf  of  a  Normal  random 
variable  with  mean  p*  and  standard  deviation  a*,  then 
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**(x)  4  Trex'p(-YLx-;j*J)’ 


(6) 


where 


_  15  tt 
16V§a* 


This  is  the  standard  Logistic  distribution  with  parameters  0  and  1 . 


We  introduce  the  chanae  of  variables 


6.  =  e  Y  i 


.  .-Y,<  Vw.*> 


where 


► 


(7) 


y.  =  ]h L 


16^3  0/ 


and  reformulate  the  NLP  using  the  approximations  discussed.  The 
result  is  given  in  Theorem  2  and  (17).  Since  the  transformation 


•  •  •  I 


is  one-to-one,  infinitely  di fferentiable, and  invertible  everywhere, 
the  problem  as  formulated  in  6  is  equivalent  to  the  problem  as 
formulated  in  S.  We  note  that  0^>O  for  any  finite  .  It  is  clear 
from  the  formulation  of  the  problem  that  the  optimal  S  is  finite, 
and  so  the  optimal  6  is  positive.  Thus  we  do  not  need  to  add  the 
constraint  6^0. 

Theorem  2.  If  the  approximations  (2)  throuqh  (6)  and  the  change  of 
variables  (7)  are  used  in  the  expressions  for  C^(S^.D^)  and  P.  (S^.D^), 
and  If  we  denote  by  (6^ )  and  P^ (0^ )  the  resulting  approximations 
(without  policy-independent  terms  and  afunction  in  that  is  o ( 1 )  as 


26 


S.-*«>),  we  have 

i 


(i)  c^e.) 


(i+e.)pi 


e. 

i 


(ii) 


W  =  n-Pi) 


where  6. 


Proof.  We  suppress  the  subscript  i. 

(i)  If  the  approximations  (2)  through  (5)  are  used  in  the  expression 

for  C(S,D),  and  C(S)  denotes  the  resulting  approximation  without  policy- 
independent  terms,  then  we  have 


r  D  S-y  S  -| 

c(s)  =  f> Uf  f  (s-y-x)<t»*(x)^  dxdy  +  h  j  (S-x)$*(x)dx 

L  n  ->  oc  -oo  ^ 


-  (policy-independent  terms). 


Now 


D  S-y 

f  f  ( S-y-x)d>*(x)dxd y 

'  J  D  od  D  ® 

=  f  J (S-y-x)$*(x)dxdy  -f  J  ( S-y-x)<t>*(x)dxdy , 

°  o  S-y 


o  - 


provided  these  integrals  exist. 

We  show  that  they  exist  by  evaluating  them.  First, 


(8) 


27 


To  evaluate  the  other  integral  in  (8)  we  define 


f  f 

C  J  J 


e(S)  c  'J  J  (S-y-x)<f>*(x)dxdy 
o  S-y 


■If 

o  S-y 


(x+y-S  )<}>*(  x)dxdy. 


It  is  straightforward  to  verify  that  on  the  domain  of  integration 


0  <  ( x+y-S )4>*(x)  <  x<f>*(x). 
Since  u*  <  «  and  D  <  °°,  we  have 


0  -00 


x4**(x)  |dxdy  <  <®, 


and  so  the  Dominated  Convergence  Theorem  (Hoffman  [1975, p.  331]) 
justifies  taking  the  derivative  of  c  under  the  integral.  This  gives 

D  00 

e'(S)  =  -J  J  $*(x)dxdy 
os-y 

=  - j  [l-$*(S-y)]dy 
o 

2  ■D0,ny<D{1"r(S‘y)) 

*  ( 

=  -D  [  1 -3E* ( S - D )  ] . 

Obviously  c  1  (S)<  0,  so 

0  >  Jinn  e  '  (S)  *  Jim  {-D[1-E*(S-D)  ]  )  =  0, 


which  implies  J|m  c'(S)  =  0.  Thus  we  have  e($)  =  a+o(l)  as  S  •*  for 
some  constant  a. 

By  assumption  v/e  do  not  include  the  o(l)  function  in  the  approximate 
cost  function,  and  so  we  have  that 


28 


C(S)  =  h ( 1 -p  )S 


*hp/ 


(S-x)$*(x)dx, 


(9) 


where  we  have  dropped  all  policy-independent  terms.  Using  the 
approximation  (6),  we  have  that 


/(S-x)*‘(x)dx  ‘/fce-v* ixvnz 1  dx- 

—  oo  —oo  I 

After  the  change  of  variables  X=e"Y^x-lJ*\  this 


becomes 


where  &  is  given  by  (7).  Then 


(ID 


and 


/j°!X  dx=^§e  ♦  logs  -  log(l+6) , 

J  (1+X)Z  10 

e 


(12) 


and  so  (7),  (9),  (10),  (11),  (12)  and  the  definition  of  C(0)  give 

C(6)  =  h(l-p) (p*  -  lloge) 

-^log  +^log(l+6) 

-  (policy-independent  terms) 


=  Y  log(l+e)  -  £  loge 


(ii)  If  the  approximations  (2)  through  (5)  are  used  in  the  expression 
for  P(S,D),  and  P(S)  denotes  the  resulting  approximation,  then  we  have 
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P(S)  «  P  j  J I*(S-y)  dy  +  **(S)j 

D 

=  ^  yi*(S-y)dy+  p£*(S). 


Using  approximation  (6)  we  have 


(13) 


-/^.e-Yis-y-Vi 


dy 


0  0 


Making  the  change  of  variables  Y  =  ^  this  becomes 

-Y(s-u*) 

1  f  ~ 
y  J  ^1- 


dY 


■y(sv) 


Then 


/  A  =  Y  "  log(l+eY) 
J  l+eT 


so  (7),  (13),  (14),  (15)  and  the  definition  of  P(0)  give 


P(e).^[,D  +  ,o9(J^)]  +  ^ 
[u  ^,03(t^|)] 


=  (1-P) 


+  T+e  D* 


(14) 


(15) 


Let 


-v. (D.-u.*) 
e.u  =  e  1  1  1 


(16) 


so  S^>D.j  is  equivalent  to  0^<0^u.  Then  we  have  formulated  the 
following  NLP,  the  solution  to  which  gives  an  approximately  optimal 
policy  for  our  multi-item  inventory  system  via  the  change  of  vari¬ 
ables  (7); 
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N 

minimize  C(0)  =  7  C^e^) 


i  =  l 
N 


subject  to  P(e)  W^P.te.J-a]  >  0,  k 

i  =  l 

ei"eiU 


(17) 


We  use  nonlinear  programming  techniques  to  solve  (17).  Specific¬ 
ally,  we  derive  first-order  necessary  conditions  for  optimality  using 
the  familiar  Kuhn-Tucker  conditions.  Although  the  objective  function 
is  convex,  the  feasible  region  in  general  is  not.  Therefore  the 
necessary  conditions' may  not  be  sufficient  to  guarantee  even  local  opti¬ 
mality  (Luenberger  [1973,  Sections  6.4,  10.6]).  We  show  that  for 
sufficiently  large,  however,  the  Kuhn-Tucker  conditions  are  sufficient 
for  (at  least)  local  optimality.  We  also  describe  a  computationally 
efficient  algorithm  for  finding  an  optimal  policy  when  it  exists,  and 
compare  the  performance  of  these  policies  with  policies  generated  by  an 
existing  method. 

The  following  formulas  are  needed. 

Lemma  3.  For  i=l,...,N, 


(a)  ci  1  (e^ ) 


(b)  ci"(ei)  = 


(C)  Pi,(e.) 


h. 

1 


h. 

1 


1+(1-P.)e. 

(1+ei)ei  j 


(i-e^e.2  +20^1 

(l+e^e.^ 


rViVj  M H>i) (V_1_) ]  _+  [DjPiY1-6i  +  (1-Pi)(^i-D3ei 
YiDi(i+e.)2(i+6.ei) 
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Proof  (a)  Differentiating  £^(0.)  and  suppressing  i  gives. 


c’(e) 


ii  f_p_ .  ll 
y  [l+e  ej 

.  h  fltlbelel 

y  L  (i+0)e  J  • 


(18) 


(b)  Differentiating  (18)  gives 


C"(8) 


h 

y 

h 

Y 


-(l+e)2  e2_ 

'(i-P)2+2e+il 

(i+e)2e2  J 


(c)  Differentiating  P^(6^)and  supressing  i  gives, 

<’9> 

-  (l-p)(l-3)  ,  _P _ 

'  yD(l+e)(l+50)  ■  (1+0)2  , 

«  .  (l-p)(6-l  )(l+e)  +  Dpyd+se) 
yD(1+6)2(1+60) 

=  .  [Dp>  +  (l-g)(6-l)]  +  DW  +  (!-»)(<-!)]$  Q> 
yD(l+e)2(l+60) 

Theorem  4.  Suppose 

P(0U)  s  0.  (20) 

Then  0*=0U  is  globally  optimal  for  (17). 

Proof.  It  is  straightforward  to  verify  (since  pi  is  decreasing  as  a 
function  of  D^ )  that  (making  use  of  (5)) 

D.  *  u •  implies  0  <  p .  <  .  (21) 

1  1  I 

This  with  Lemma  3(a)  implies  that 

C.’te.)  <  0  whenever  0^  *  0^u, 


(22) 
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and  so  C(&)  is  decreasing  in  each  argument  for  0^©^.  But  by 
hypothesis  0=©u  is  feasible,  so  8=0U  is  a  global  minimum  of  (17)  0. 


In  order  to  investigate  the  case  when  P(6U)<0,  we  make  use  of 
the  Kuhn-Tucker  necessary  and  sufficient  conditions  for  optimality. 

To  use  these,  however,  one  must  verify  that  some  constraint  qualifi¬ 
cation  holds  at  the  optimal  point.  We  show  that  if  P(6U)<0,  then  every 
feasible  point  is  a  regular  point,  that  is,  the  gradients  of  the  active 
constraints  at  the  point  are  1 inearly-independent.  Luenberger  [1973, 
223-227]  shows  that  this  constraint  qualification  is  sufficient  to 
guarantee  that  if  the  Kuhn-Tucker  necessary  and  sufficient  conditions 
hold  at  a  point, then  that  point  is  an  (at  least  local)  optimum. 

Lemma  5.  If  P(8U)<0,  then  any  0  feasible  for  (17)  is  a  regular  point 
of  (17),  that  is,  the  gradients  (at  6)  of  the  constraints  active  at  0 
are  all  1 inearly-independent. 

Proof.  Let  0  be  feasible  for  (17).  Without  loss  of  general ity,  suppose 
that  of  the  constraints  e.*©^,  i=l,...,N,  the  first  j  are  active  and 
the  last  N-j  are  inactive.  Since  by  assumption  6=eu  is  not  feasible,  w 
have  N-j>l.  Let  T  denote  the  matrix  which  has  as  its  columns  the  grad¬ 
ients  (at  e)  of  the  constraints  active  at  0.  If  the  constraint  P(e)^0 
is  inactive,  then 


T  = 


1  0 

j 

Q..1.1 

N-j  0  •  •  •  0 


which  certainly  has  1 i nearly-independent  columns.  If  the  constraint 
P(0)iO  is  active,  then 
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j  1 


1 

0 

WT  P1 ' (01 ) 

j 

T  = 

0 

1 

w2P2'(e2) 

N-j 

;  o... 

0 

WNPN'(0N) 

Certainly  the  first  j  columns  of  T  are  1  inearly-independent.  If  the 
last  column  were  a  linear  combination  of  the  first  j,  then  (since 
N-jil)  WNPN'(eN)=0.  Now  it  is  straightforward  to  verify  that 


D^>0  implies  6.>1, 

so  this  with  Lemma  3(c)  and  (21)  implies  that  PN'(eN)>0. 

And  since  W„  >  0,  we  have  a  contradiction,  and  so  the 
columns  of  T  are  1 inearly-independent,  which  implies  that  6  is  a 
regular  point  of  (17)  □. 


Theorem  6.  Suppose  P(eu)<0  and  for  i=l,...,N,  let 


Ai  *  Ai<Di> 
Bi  * 

Ei  ‘  Ei<Di> 


Dihi[2*5.-p.] 
Djh^+O -»,)(«>!)] 


Fi  ■  w 

Gi  *  Gi<Di> 

Hi  *  w 


Dih,n-pi)«1 

+  (1-fPtij-l)] 

+  O-filWi-D] 


Me-)  = 
1  '  1  1 


VBiVEi*j2tFi*i3 

Gidi+Hi®i2 


(23) 


Then  the  first-order  Kuhn-Tucker  necessary  condition  that  0*  be  a 
minimum  of  (17)  is  the  following: 
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There  exists  a  real  number  A*gO  such  that 

(a)  Xi(0i*)  =  A*  when  e^*<e..u 

A^a.*)  >  A*  when  0^ *=6^ U 

(b)  p(e*)  =  o. 


Proof.  As  proved  in  Luenberger  [1973,  pp. 232-234],  given  the 
constraint  qualification  guaranteed  in  Lemma  5,  a  first-order 
necessary  condition  for  e*  to  be  a  minimum  of  (17)  is  that  e*  be 
feasible  for  (17),  and  that  there  exists  A*>0  and 
5*  =  (t1*,...,cN*)  2  0  such  that 


vL(e*,A*,$*)  =  ?c(e*)  -a*?p($*)  -  4*i  *  o 
x*p(e*)  =  o 
;*T(e*-eu)  =  o, 

where  I  is  the  NxN  identity  matrix,  and  a  superscript  T 
transpose. 

By  assumption  P(eu)<0,  and  clearly 
N 

P(0)  -T^O-g)  >  0; 

i=l 


►  (25) 

* 

denotes 


(26) 


since  it  cannot  cost  less  to  give  a  higher  service -level ,  we  can 


assume  that  at  optimality  P ( 6*}  =  0.  Now  Lemma  3(c),  (21  )  and  (23) 
imply  that  for  all  i , 

P.’  (0..  )<~0  whenever  (27) 

and  so  by  (22)  we  have  that  (25)  is  equivalent  to  the  following 
condi tion: 

There  exists  A*j0  such  that  for  all  i, 
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(a) 


Ci '  (0i*) 

wiPi'(ei*) 


=  \* 


when 


V<  ei 


(b) 


wiPi,(ai*) 

p(©*)  =  o. 


>  X*  when  6..*  =  6..  U 


We  note  that  =  ^'(6^)  -  X* Pi '(©>),  for  each  i. 
proof,  we  need  only  to  show 


N 


►  (28) 


/ 

To  complete  the 


VfV) 

W'  =  w1Pi • (ei*)  ' 

We  suppress  the  subscript  i.  By  Lemma  3(a)  and  (c). 


c‘!e!  -  /-  Jl\ 

ri+O-fJel 

( .vn) 

[(i+e)2(i^)l 

wp'(e)  ‘  \  y  ) 

 (1+6)6  

V 

L  g+hs  J 

After  factoring  (1+0)  out  of  the  numerator  and  denominator,  the 

2 

numerator  becomes  Dh[l+(1 -p)©][1+(6+1 )©+£©  ],  which,  when  multiplied 

?  3  2 

out,  is  A+B0+E©  +F©  ,  and  the  denominator  is  G0+H0  .  Thus  we  have 

ci ’(e.) 

wiPi1(eiT  =  W*  (29) 

as  was  to  be  shown  C. 

Although  the  first-order  Kuhn-Tucker  condition  is  necessary  for 
(at  least  local)  optimality,  in  general  it  is  not  sufficient.  In 
Theorem  7  we  give  a  condition  sufficient  to  guarantee  that  the  Hessian 
of  the  Lagrangian  of  (17)  is  positive-definite,  a  strong  second-order 
sufficient  condition  (see  Luenberger  [1973,  p.235]).  In  Theorem  8 
we  give  conditions  in  terms  of  which  guarantee  that  the  hypotheses 
of  Theorem  7  are  satisfied. 

TheorenW.  Suppose  P(eu)<0  and  that  there  exists  X*>0  so  that  the 
necessary  conditions  (24)  hold  at  0*.  If  also  A.(0.)  is  decreasing 
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for  e.<0.u,  then  the  Hessian  matrix  VAJ1L(e*,X*,^*)  is  positive- 
definite,  which  implies  that  (24)  is  sufficient  for  (at  least  local) 
optimality. 

Proof.  We  first  note  that  Lemma  3(b)  and  (21)  imply  that 

Ci"(6j )  >  0.  (30) 

We  have 

L(e* ,**,?*)  =  tf2c(e*)  -  x*v2p(e*) 

C,  "(*,.*)  -  "(*.*)  0 

o  cN"(0N*)  -  x*wNPN-'(eN*)  . 


Let  Q.  denote  the  i^  diagonal  entry  of  7-^.L .  We  show  that  Q->0 
for  each  i,  as  this  implies  that  V^L  is  positive-definite.  If 
P^'U^JsO,  (30)  implies  Q^O,  since  X*,  W.>0.  If  Pi"(6i*)>o,  (24) 
implies  Q.feC-"(ei*)  -  [X..^.*)]  W^Pi "(0^ *) ,  which  by  (29)  is 

ci’(0i*) 

«i  4  V(V>  -ptT37Tpi"(f,i*K 

or 

Q  >  Pil(ei*)Ci,'(6i*)  ~  cj  ‘  (6i*)Pi  “(Qj*)  # 

1  p.'(e.*) 

By  (27)  the  denominator  is  negative,  so  we  can  conclude  the  proof 
by  showing  that  the  numerator  is  negative.  From  (29)  we  have 

x.  •(..*)  =  pi,(v>ci"(y>  - 

1  1  wi[Pi-(e.*)]2 


(31) 


By  hypothesis  0>X^,(©i*),  so  the  numerator  of  (31)  is  negative. 
Since  W^>0,  we  have  Q^>0,  as  was  to-be  shown. 


Luenberger  [1973,  p.235]  proves  that  given  the  constraint  quali¬ 
fications  guaranteed  in  Lemma  5,  the  first-order  necessary  conditions 
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(24)  are  sufficient  for  (at  least  local)  optimality  provided 


V.  L(0*,A*,c*)  is  positive-definite  on  the  subspace  tangent  to  the 
constraints  tight  at  0*.  We  have  shown  that  V  L(0*,A*,c*)  is  positive- 

ud 

definite,  which  concludes  the  proof.  We  note  that  to  prove  this 
theorem  for  more  general  functions  and  it  is  sufficient  that  C- 
be  convex,  and  P^  and  A..  be  decreasing  □. 

Theorem  8.  For  i=l,...,N,  if  is  sufficiently  large,  then  X^(6^)  is 
convex  and  decreasing  when  0^<0.jU. 

Proof.  As  usual  we  suppress  i.  We  write  X(0)  in  partial  fraction  form 
making  differentiation  easier.  We  then  show  that  for  D  large,  the 
coefficients  of  A(O)  are  such  that  A'(9)<0  for  6<0U.  By  long 
division. 


Me)  -ft  *  Blf  ♦ 

n  Hc 


p-WH-rol  ,  A 

5  .  \  H  / 


H 

e(G+H9) 


X6+A 


XG-AH 


For  any  X,  y  ^  +  g(G+H6)  ’ 


and  so 


where 


>(*)=  +  +  n__ 

H  HZ  G  GHZ(G+Hp) 

n  =  GBH2-  G2(EH-FG)  -AH3  . 


By  (21)  and  (23)  we  have 

A,  B,  E,  F,  G,  H  >  0, 

so  certainly  A(e)>0  whenever  e>0  (the  change  of  variables  (7)  is 
such  that  6>0).  Differentiating  (32)  gives 


(32) 


(33) 
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and 


n 

GH(G+HO)2 


>•"(«) 


_2A  +  2n  . 
G63  G(G+H6)3 


We  first  show  that  ,V'(0)  >  0  for  D  sufficiently  large  by  showing 
r,  >  0  for  D  sufficiently  large  (since  (33)  implies  all  other  quantities 
are  positive).  We  use  the  following  notation. 
f(x)  ~g(x)  if  lim  =  1. 

To  examine  the  behavior  of  r\  for  large  D,  we  examine  the  behavior 
of  A,  B,  E,  F,  G  and  H  for  large  D.  The  following  are  easily  inferred 
from  their  definitions: 


lim  6  =  <>■-> 

limp  =  0 

or** r 


lim  Dp 

cr° 


v- 


(34) 
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They  imply 

A  ~  Dh 
B  -  Dh$ 

E  ~  2Dh$ 

F  -  Dhtf 
G  ~  W$ 

H  ~  W(|i-y+l  )£  , 

and  so  n  ~  (W$)(Dh*)l.'2(HY+1)262 

-  (Dh)(W3(pT+l)V). 


W262[(2Dh5)(VJ(nY+l)^)  -  (Dh*)(W$)] 


(35) 


The  limiting  behavior  of  n  is  determined  by  its  highest  order  term, 


which  is  D6  ,  SO 


n  ~  [W3h(-J>yf1  )2-2W3h(nY+l )  +  W3h]  D$4 

•3  2  4 

=  W^h^ly)  D5  . 

Since  W3h(vy)2  >0,  (34)  implies  lim  w3h(yy)2D$4  =  »,  so  1  in  n  =  a’  al  so . 

D-W  Q->cc 

Thus,  for  D  sufficiently  large,  n  >  0,  and  so  A"(6)  >  0,  as  was  to  be 
shown.  Hence  A(0)  is  convpx .  To  show  that  A'(e)  <  0  for  9  <  6U, 
it  is  sufficient  to  show  X ' ( 0U )  <  0.  We  have 


x '  ( e u )  *  £  - 


- -  * 

GH(G+H0U) 


f 


where 


A _ 

G(nu)2  ’ 


Since  for  n>0  large  D,  we  need  only  show  that  £<0  for  large  D  to 
conclude  X'(e)  <  0  for  6  <  Qu.  Relations  (35)  imply  that 


DhjT  Dh 

'  ~  W(l*Y+l  )~6  ’  JW* 
W6%— 
6 
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-  Dh  l  1  6  \ 

W  \^y+l  ’  e2»* )  • 

Relations  (34)  imply  the  quantity  in  parentheses  tends  to  -°°  as  D-**,  and  so 


lim  t,  -  lim 


-oc. 


Thus,  C,  <  0  for  D  sufficiently  large,  and  so  the  proof  is  complete  □. 

Thus,  Theorems  7  and  8  imply  that  if  is  large  enough  that  X^ 
is  decreasing,  then  the  first-order  conditions  (24)  are  sufficient 
for  (at  least  local)  optimality.  We  exploit  this  in  our  algorithm, 
which  is  described  in  the  next  section. 


3.  GKD  ALGORITHM 


Our  algorithm  is  based  upon  local  duality  theory  of  nonlinear 
programming  (see  Luenberger  [1973,  pp. 312-320]) .  Suppose  x*  is  a  local 
minimum  (at  which  the  constraint  qualification  described  in  Lemma  5 
holds)  of  the  nonlinear  program 

minimize  f(x) 


(36) 


subject  to  g(x)=0, 

with  associated  Kuhn-Tucker  constraint  multiplier  X*.  Writing  the 
Lagrangian  L(x,>)  =  f(x)-Xg(x),  the  NLP  dual  to  (36)  is  defined  as 


where 


maximize  i/(\) , 

X  1 

/ 

M\)  =  minimum  L  (x,X) . 
x  near  x* 


(37) 


(We  mention  that  (36)  and  (37)  can  be  written,  respectively,  as 

min  max  L(x,X)  and  max  min  L(x,X),  which  illustrates  their  dual  nature), 
x  X  X  x 
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If  L(x,X*)  is  twice  continuously-differentiable  and  convex  at 
x*  (and  hence  convex  in  a  neighborhood  of  x*),  then  the  following  can 
be  shown  to  be  true: 

in  *(  A)  is  twice  continuously-differentiable  and  convex  at  X*. 

(2)  X*  uniquely  solves  (37). 

(3)  If  we  denote  by  x(X)  the  (unique)  value  of  x  minimizing  the  right- 
hand  side  of  4>(X)  in  (37),  then  x*=x(X*). 

Thus  rather  than  solving  (36)  one  could  solve  (37),  and  do  this  by 

solving 

^fr(X)=0, 

which  can  be  shown  to  be  equivalent  to  solving 

g(x(X) )=0.  (38) 

Generally  this  approach  is  not  useful  computationally  because  the 

function  x(X)  is  often  difficult  to  compute.  This  method  may  be  quite 

efficient,  however,  if  x(X)  is  reasonably  easy  to  compute  and  the  dimension 

of  X  is  smaller  than  that  of  x.  One  problem  with  such  structure  is  the 

"Generalized  Knapsack  Problem,"  which  is 

N 


mi 


inimize  f(x)  =  2Zf--(x*) 
i=l  1 


subject  to  g(x)  =  Xj  g^Xj)  =  0. 


►  (39) 


/ 


Clearly  (17)  is  a  Generalized  Knapsack  Problem  with  upper-bounded  variables. 
These  bounds  do  not  change  the  essential  structure  of  the  problem  or  its 
solution  by  solving  its  dual,  and  the  details  are  left  to  Theorem  9. 

For  (39),  x(X)  is  found  by  solving  the  system  of  equations 

VxL(x,X)=0, 
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where 

x,L  -  x, <*,“). 

For  convenience,  we  assume  that  the  items  are  reordered,  if  necessary, 
so  that 

V  s  >2L  -•••*  hl.  (4i) 

Theorem  9.  If  either  P(eu)  ^  0  or  all  are  sufficiently  large  so  that 
the  hypotheses  of  Theorem  8  hold,  then  the  following  algorithm 
is  well-defined  and  converges  to  a  (possibly  local)  minimum  of  (17). 

GKD  Algorithm. 

(1)  If  P(0U)  i  0  set  i =N ,  and  go  to  step  (6).  Otherwise  set  i=0. 

(2)  Set  XL  -  X^,  (•  maxjx,!;, . xf(L},by  (41)). 

(3)  Define  the  function  P.  on  [a  ,+«>)  by 

VM  =  p(e1U2...,eiu,ti+1(A),...,tN(A)). 

(4)  If  P.(AL)  <  0,  go  to  step  (5). 


(5) 

(6) 


Otherwise  set  i-i+1,  go  to  step  (2). 
Solve  P.(A)  =  0  for  A*c[AL,  »). 


Set  6j*  - 


6; 


tj (a*) 


S  .*  =  y . 

J  V3 


*  -  log0.* 


s.*  =  S.* 
J  J 


1 

sr. 

j 

D.. 

J 


j=l 

j=i+l,...,N  , 

V 


Proof.  Steps  (1)  and  (2)  are  certainly  well-defined.  By  definition  of 
AL,  A  5-  >L  is  in  the  domain  of  t.+1  ( A ) ,  — .  so  these  are  all  well- 

defined.  Hence  so  is  P(A),  and  so  step  (3)  us  well-defined.  Step  (4)  is 
well-defined,  except  for  the  possibility  that  i  may  become  larger  than  N. 
But  if  i=N,  then  we  have  P(6^u,. . . ,0^u)  >  0,  which  is  equivalent 
to  P($u)  >  0.  But  step  (1)  guarantees  P(6U)  < 0,  so  i  will  never  exceed  N. 
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A  | 

Step  (5)  requires  justification.  We  show  that  '  ( X ) >0  on  [A  ,+°°) 

A 

and  P.j (+°°)>0.  These  with  step  (4)  imply  that  there  is  a  unique  solution 
in  [X*- ,+°°)  to  P.(A)=0.  We  have 

P,(>)  =  P(B,u,...,0^,  ti+1(X) . tN  (X)). 


so 

N 

VW  5  HWjPj'(ej>  tj'O).  wher®  Bj  *  tj(X), 
j=i+l 

A 

by  the  chain  rule.  Relation  (27)  implies  that  P.'(e.)<0,  and  clearly 

-  J  J  i 

tj'(X)<0,  so  indeed  P. '(\)>0.  Stage  i  is  the  first  time  that  Pi(AL)<0. 
If  i=0,  then 

p.H  =  p(t1H . tN(»)) 

=  P(0, . . . ,  0)  >  0, 
by  (26).  And  if  i  >1,  then  Pi_1(At")  >  0,  so 

o  <  P.j (al)  =  P(e1u,...,ei^,ti,...,tN) 

P(H1  »•••>?.,•  »t.j+i , . . .  ,t^) , 

since  at  stage  i-1,  AL=A.jL,  and  t(X.L)  =  6-U. 

Hence,  0  <  . . . ,0i U,ti+1 , . . . ,tN) . 

We  have 

9  P(o1u,...,eiu,t.+1(»),...,tN(«')) 

=  p(o1u,...,Piu,o,...,o). 

A  A  . 

By  (27)  P  is  decreasing  in  each  argument,  so  P^ (~)  >  P.  ^(\L)  >  0,  as 
was  to  be  shown. 

Step  (6)  and  (7)  are  well-defined,  so  the  entire  algorithm  is  veri¬ 
fied.  If  P(6U)^0  at  step  (1),  then  by  Theorem  4,  e*=eu  is  globally 
optimal  for  (17).  If  P(0U)<O  at  step  (1),  the  fact  that  the  functions 

Aj  and  tj  are  inverses  with  the  fact  that  A.L=A.(e.u)  imply  that  the 

J  J  J 
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algorithm  terminates  with  0<6*<eu  and  A*>0  satisfying 

(a)  \j(y)  =  Aj1,  j-l.. ...i 

AjOj*)  =  A*  ,  j=i+l . N 

(b)  p(e*)  =  o. 

But  (41)  implies  X,(©.*)  >  AL  for  j=l . ,1 ,  and  certainly  XL  *  A*. 
Moreover,  since  6^*=©^,  for  j=l,...,N,  e*  satisfies  (24).  By  Theorems 
6,  7  and  8,  these  equations  imply  the  local  optimality  of  6*  for  (17). 

Recall  that  We  reformulated  the  multi-item  inventory  problem  using 
the  one-to-one  change  of  variables  given  by  (7).  Our  GKD  Algorithm 
finds  an  optimal  6*  for  (17),  and  this  gives  rise  to  an  optimal  S* 
via  the  inverse  of  the  change  of  variables  (17),  which  is 

logGi 

S.  =  u* - 1* 

i  y. 

as  given  in  step  (6)  □. 

The  algorithm  is  computationally  efficient  in  that  it  only  requires 

A  /-s 

solving  P1-(A)=0,  where  P^(A)  is  a  real-valued  function  of  one  real 

variable.  The  computation  of  the  functions  t.(A)  must,  however,  be 

J 

done  numerically  (the  details  are  in  the  Appendix  to  this  chapter). 
Moreover,  the  hypotheses  of  Theorem  8  are  easy  to  verify  and  guarantee 
that  the  algorithm  is  well-defined  and  converges  to  a  minimum  of  (17). 

The  proof  of  Theorem  8  shows  that  is  sufficient  to  check  that 
A. ' (9.u)  <  0  for  each  i.  The  extension  of  the  algorithm  when 
A.j'(©.jU)  >  0  is  discussed  in  Section  5. 
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4.  POLICY  PERFORMANCE 

In  order  to  illustrate  the  performance  of  the  (s,S)  policies 
generated  by  our  algorithm,  we  performed  the  following  test.  We  examined 
fifteen  2-item  inventory  systems  with  the  parameter  settings 
K1  =  K2  =  32 

^1  =  ^2  =  ^ 


h1  =  .05,  =  1  ,  10,  25  ,  50,  100 

u  =  3  *1  =  1.  1 .  3 
'  V,  3 
2.  2 
!i-  _  g 

1*1  ^  ' 

W1  =  W2  =  1  * 


We  chose  these  parameter  settings  for  three  reasons 

1.  Most  of  them  give  rise  to  large  D  (much  larger  than  1.5y). 

Since  our  approximations  and  theorems  generally  hold  for  large 
D,  it  should  perform  well  on  these  systems. 

2.  These  parameters  are  realistic,  in  that  they  reflect  costs 
associated  with  many  actual  inventory  systems. 

3.  As  in  Mitchell  [1982J,  our  algorithm  should  perform  effectively 

when  there  are  both  "expensive"  and  "inexpensive"  items 
in  the  system.  This  holds  for  those  systems  we  considered 
^2 

with  r—  large. 
nl 

The  comparison  we  used  assumes  the  underlying  distribution 
is  Negative  Binomial  (the  discrete  version  of  a  Gamma  distribution, 
which  is  skewed  for  our  parameter  settings),  specifies  a  service- 
level  a,  and  uses  the  Veinott-Wagner  exact  algorithm  (Veinott 
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and  Wagner  [1965])  to  compute  optimal  policies  for  both  items  in  each 
system  using  the  same  service-level  a.  We  ran  our  algorithm  using 
a.  Then  we  compared  the  holding  plus  set-up  costs  generated  by  the 
Veinott-Wagner  algorithm  and  our  algorithm.  The  results  are  displayed 
in  Table  1 .  , 

We  make  the  following  observations. 

1.  The  service-level  a  is  not  the  same  for  all  runs,  but  varies  from 
.89  to  .97.  The  average  service-level  for  the  fifteen  runs  is  .93. 
There  are  two  reasons  for  this.  First,  for  some  of  these  items, 
the  policy  with  the  smallest  service-level,  (0,0),  has  a  very 
high  service-level,  so  that  the  system-wide  service-level  is  high. 
For  other  items,  the  smallest  service-level  is  small.  It  turns 
out  that  to  specify  a  uniform  service- level  is  to  specify  a=.97 
(since  one  of  the  systems  has  this  as  the  smallest  service-level 
for  both  items),  which  is  unreasonably  large.  The  second  reason 
is  that  the  logistic  approximation  we  use  in  our  algorithm  does 
not  always  generate  policies  with  service- level  a  when  the  actual 
distribution  is  Negative  Binomial.  However,  preliminary  tests 
on  larger  systems  (32  items)  seem  to  indicate  that  this  problem 
is  less  severe  for  larger  systems. 

2.  Our  approximation  for  is  almost  exact  for  every  single 
item.  This  suggests  that  the  asymptotic  approximations  we 
are  using  should  be  fairly  accurate.  For  all  items  is 
sufficiently  large  to  guarantee  that  the  hypotheses  of 
Theorem  7  are  satisfied,  so  all  policies  are (possibly  local) 
minima  (not  saddle  points). 


48 


3.  Over  half  the  systems  had  an  item  with  s  <  0  for  the  policy 

*  generated  by  the  Veinott-Wagner  algorithm.  Our  algorithm 

can  only  consider  policies  with  s  >  0,  but  its  policies  still 
compared  very  favorably  against  those  generated  by  the  Veinott- 
Wagner  algorithm. 

4.  As  expected,  there  was  a  cost  reduction  for  all  systems,  except 

the  ones  with  identical  holding  costs.  Moreover,  the  greatest 

h 

cost  decreases  were  for  those  systems  with  large. 

This  initial  investigation  demonstrates  that  inventory-operating 
costs  can  be  reduced  significantly  when  a  system-wide  service  level  is 
specified,  rather  than  specifying  an  identical  service-level  for  each 
i  tern . 
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50 


25 


10 


40 


26 


33 


23 


14 


Table  1 


31 


25 


19 


10 


W, 


Percent  decrease  in  holding  +  set-up  costs 
Assumptions : 

K]  =  K2=32  ,  k i  =  k 2=0  ,  h^.05 
Z  Z 

!i  =  °-l.  9  H.l  n  =3 

Ml  V2  ’  “l^2  2  •  pl 
all  demands  Negative  Binomial. 
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5.  GKD  ALGORITHM  MODIFICATION  WHEN  D* IS  SMALL 
Theorem  8  of  Section  2  contains  the  most  general  condition  under 
which  we  were  able  to  prove  (see  Theorem  7)  that  the  GKD  Algorithm 
converges  to  a  (possibly  local)  minimum  of  (17).  This  condition  is 
that  for  all  items  i=l,...,N  the  functions  X^e^)  are  decreasing  on 
(0,  0,jU).  We  give  our  recommendations  when  this  condition  may  not 
hold  for  some  items  in  the  inventory  system. 

Recall  from  the  proof  of  Theorem  8  that  A. (0..)  is  convex  on  (0,  6^u) 
when  n.j>0  and  is  decreasing  on  (0,  0^u)  when  both  n^>0  and  £.j<0.  More¬ 
over,  lim  r).=o*  and  lin  £•=-'*-,  so  X.(0.)  is  decreasing  on  (0,  0.u)  when 
d.-«x  1  iii  i 

is  sufficiently  large.  The  GKD  Algorithm  (given  in  Theorem  9)  exploits 
this  property  by  assuming  the  existence  of  the  inverse  function  of  X.. , 
which  is  denoted  by  t^ . 

The  question  is  what  to  do  when  D^  is  not  large  enough  to  guarantee 
^>0  and  £.<0.  In  either  case  it  is  impossible  to  guarantee  that  is 
decreasing,  and  hence  invertible,  so  step  (3)  of  the  GKD  Algorithm 
may  not  be  well-defined.  The  algorithm  may  converge,  but  there  is 
no  guarantee  that  it  will  converge  to  a  local  minimum. 

We  recommend  the  following  modifications  when  n^<0  or  ^>0.  If 
^.>0  but  n^>0,  then  X..  is  still  convex  but  not  decreasing  on  (0,  0^u). 

It  is  always  true  (and  easy  to  verify)  that  X *  ( 0 )=-«',  so  this  case  is 
illustrated  in  Figure  1. 

In  this  case  we  recommend  that  the  domain  of  X^ (0^)  be  shortened 
to  (0,  e.jU),  where  e^u  is  the  largest  value  of  0^  for  which  X^  ( 0^ )  is 
decreasing.  This  has  the  effect  of  introducing  the  additional  constraint 
S^>sA  ,  where  Since  the  are  lower  bounded  anyway  by  , 
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this  constraint  does  not  change  the  structure  of  the  problem  at  all. 
In  particular,  the  6KD  Algorithm  will  still  converge  to  a  (possibly 
local)  minimum. 


Computationally,  ft.u  can  be  found  by  increasing  the  given  lower 
bound  D.j  on  S/  a  unit  at  a  time  until  6^u  is  found.  This  procedure  is 
described  in  detail  at  the  end  of  this  section,  and  computational 
experience  on  typical  32-item  inventory  systems  is  reported  in  the 
next  chapter. 

Although  we  have  in  this  case  a  recommendation  that  still  assures 
that  the  GKD  Algorithm  will  converge,  the  case  when  n^<0  poses  greater 
difficulty.  The  may  not  be  convex,  in  which  case  we  have  been 
unable  to  find  a  modification  of  the  problem  that  guarantees  convergence 
to  a  local  minimum.  We  essentially  ignore  this  case  and  use  the  algor¬ 
ithm  as  it  is  (making  the  recommended  modification  when  C->0).  Our 
computational  experience  reported  in  the  next  chapter  seems  to  indicate 
that  this  is  not  a  problem  in  many  cases.  Since  is  easy  to  compute 
prior  to  using  the  GKD  Algorithm,  we  recommend  doing  so  in  order  to  be 
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alerted  to  possible  problems  when  n^<0. 

In  light  of  all  this,  the  only  recommended  change  involves  the 
computation  of  6^°.  The  6KD  Algorithm  is  just  as  given  In  Theorem  9, 
with  the  following  step  inserted  between  steps  (1)  and  (2): 

(1.1).  For  j  =1 . N: 

Set  SjL=Dj 

(1.1.1) .  Set  e.jU  =  e_>j^Sj  _lJj  K 

J 

If  x. ,(aiU)<0,  go  to  step  (1.1.2). 

J  J 

Otherwise,  set  S.  *■=$.  ^+1. 

J  J 

Go  to  step  (1.1.1). 

(1.1.2) .  Next  j . 

6.  CONCLUSIONS 

We  have  formulated  a  constrained  NLP  to  solve  an  (s,S)  inventory 
model  with  a  service-level  constraint  rather  than  shortage  costs. 

This  model  has  the  advantage  that  a  manager  can  more  easily  specify 
service  objectives  than  shortage  costs.  We  derived  necessary  and  suffi¬ 
cient  conditions  for  an  optimal  (s,S)  policy.  Using  asymptotic  approxi¬ 
mations  (where  D=S-s  is  large)  and  assuming  Logistically  distributed 
demands,  we  derived  a  computationally  efficient  algorithm  to  find 
approximately  optimal  (s,S)  policies.  We  showed  that  for  some  inventory 
systems  there  is  a  substantial  cost  savings  when  using  this  algorithm 
instead  of  the  Independent  and  Identical  Service  approach  that  is  often 
used  in  practice.  We  also  gave  recommendations  to  modify  our  algorithm 
when  D  is  not  sufficiently  large  to  guarantee  its  convergence. 
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its  unit  holding  cost  h).  This  system  is  described  in  Figure  2, 
where  the  items  are  ordered  in  decreasing  order  by  hp.  Essentially 
the  first  seven  items  are  high-value  items  (have  a  very  high  annual 
dollar  demand),  and  the  other  twenty-five  items  are  low-valued  items. 

As  a  base-case,  we  specified  inventory  parameters  as  follows.  For 
i-1 , . . . ,N, 

Ki  =  24 

wi  =  J2 

k.  =  4 
1 


and 

^1 » ^5 » ^9  ♦  • • *  *  ^29  =  ^ 

^2’v6,u10’ * ' ' ’  p30  =  8 
V13*'J7  *U1 1  ’  •  •  •  ’  Ti^i  =  ^ 

^4  ’  “8’’'l  2  ’  ‘ ‘  ’  u32  = 

Of  course,  h.  is  determined  by  p.  from  the  value  (hpj^h.p.. 

We  are  interested  primarily  in  the  algorithm  performance  at  a 
servi ce- 1 e vel  of  850,  one  managers  often  prefer.  Our  criterion  of 
comparison  is  to  assume  that  the  demand  distribution  is  Negative  Bincmial 
and  compare  the  total  expected  cost  associated  with  the  (s,S)  policies 
that  the  GKD  Algorithm  produces  with  the  total  expected  cost  of  the 
policies  given  by  optimally  setting  all  the  (s,S)  policies  so  that 
every  item  has  a  service-level  of  850.  We  call  this  the  Identical 
Service  Approach,  and  it  is  an  approach  often  used  by  managers. 

We  intended  to  use  the  Veinott-Wagner  [1965]  Algorithm  (also  Kaufman 
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[1976])  to  compute  the  optimal  policies,  but  the  amount  of  computation 
required  was  impractical  because  the  optimal  policies  often  had  Si  too 
large. 

Instead  we  used  the  Power  Approximation  of  Ehrhardt  11976]  and 
[1979]  to  compute  approximately  optimal  policies.  We  chose  this 
approximation  because  extensive  numerical  investigations  in  Ehrhardt 
[1976],  and  Klincewicz  fl976a]  and  [1976b]  showed  this  to  be  both 
the  most  accurate  and  most  easily  computed  approximately  optimal 
policy  available.  We  use  the  version  of  the  Power  Approximation 
given  in  Hosier  [1981]  because  of  the  technical  improvements  made. 

The  Power  Approximation  requires  the  specification  of  shortage  costs 
rather  than  the  service-level.  Letting  a  denote  the  service-level  speci¬ 
fied  for  the  inventory  system  (using  the  Identical  Service  Approach),  we 
specified  the  shortage  costs  p.  using  the  expression  (see  Ehrhardt  [1977,  p.8]) 

P-  +  0.0695h. 

*  =  - -  . 

Pi+h, 

Ms  expression  is  an  empirical  improvement  over  the  well-known 
relationship  (see  Veinott  and  Wagner  [1965]) 


Because  both  the  GKD  Algorithm  and  the  Power  Approximation  are 
approximations,  it  was  impossible  to  realize  the  exact  service-level 
specified.  We  therefore  specified  various  service-levels,  and  inter¬ 
polated  intermediate  service-levels.  Since  we  used  the  Power 
Approximation  in  the  GKD  Algorithm  to  compute  the  values  for  D,  we 
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only  compared  the  expected  holding  costs  for  this  algorithm  and  the 

Power  Approximation,  the  expected  replenishment  costs  being  identical 

(Typically  the  expected  replenishment  cost  accounts  for  about  half  of 

the  total  expected  cost).  The  result  of  this  experiment  is  displayed  in 
Figure  3.  There  is  a  significant  decrease  in  expected  holding  costs  at 
and  near  a  service-level  of  85%. 

In  order  to  investigate  the  sensitivity  of  expected  holding  cost 
to  the  various  parameters,  we  performed  four  additional  experiments 
in  which  we  modified  one  parameter  (or  parameter  group)  at  a  time  in 
our  base-case  model.  The  same  experimental  dpsign  was  used  for  these 
experiments  as  for  the  base-case. 

Figure  4  describes  the  result  when  the  means  were  changed  to 

p2»U4» • • • *^32  ~  6* 

As  can  be  seen,  there  was  a  very  slight  degradation  in  the  GKD 
Algorithm  performance.  The  cost  savings  are  still  substantial. 

Figure  5  describes  the  result  when  the  base-case  was  altered  so 

that 

1^=8,  i  =  l , . . .  ,32. 

Again  there  was  only  a  slight  degradation  in  cost. 

Figure  6  is  the  base-case  with 

2 


The  cost  degradation  is  greater  than  in  previous  experiments,  but  the 
cost  savings  remai ns  substantial . 
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Figure  7  describes  the  change 

Uj .p2» • • • =  16 

li17*wl8’****|J32  =  8- 

The  cost  savings  are  still  substantial,  but  not  as  great  as  in  previous 
experiments . 

As  a  final  experiment,  we  did  a  worst-case  type  experiment,  in 
which  all  of  the  following  changes  were  made: 


2 


Vi»-“»Vi6  16 

Jj17,  . . .  »V32  =  8. 

This  experiment  is  described  in  Figure  8.  The  degradation  is  substan¬ 
tial,  but  the  cost  savings  are  still  large.  The  behavior  of  the  GKD 
Algorithm  policies  for  service-levels  larger  than  about  87%,  where 
the  total  cost  is  greater,  seem  to  be  caused  by  the  error  in  approximat¬ 
ing  the  Negative  Binomial  demand  by  a  Logistic  distribution  (Section  2 
of  the  previous  chapter).  We  address  this  issue  at  the  end  of  this 
chapter. 

We  make  the  following  general  observations  concerning  the  perform¬ 
ance  of  the  GKD  Algorithm  for  these  32-item  inventory  systems. 

1.  The  total  holding  cost  savings  are  substantial  at  and  near  a  service 
level  of  85%.  The  algorithm  is  an  improvement  over  the  Identical  Service 
Approach  often  used  by  managers,  and  significantly  reduces  the  overall 
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cost  of  maintaining  an  inventory  system  subject  to  a  service-level 
constraint. 

2.  In  most  cases  the  policies  produced  by  the  GKD  Algorithm 
degrade  sharply  around  a  service-level  of  86?I  or  87%,  as  is  evidenced 
by  the  sudden  sharp  increase  of  the  holding  cost.  This  suggests  that 
the  GKD  Algorithm  may  not  produce  substantial  cost  savings  when  a 
high  overall  service-level  is  specified. 

3.  The  GKD  Algorithm  seems  to  perform  best  when  there  is  substan¬ 
tial  variation  in  demands  means,  when  the  leadtime  is  not  too  large, 
and  when  the  variance  to  mean  is  not  too  small.  The  degradation  in 
performance  is  due  to  the  error  in  approximating  a  Negative  Binomial 
demand  distribution  by  a  Logistic  distribution.  This  is  especially  a 
problem  when  there  are  certain  patterns  in  the  demand  means  and  large 
leadtimes.  This  issue  is  discussed  in  detail  in  the  next  chapter. 


3.  SINGLE-ITEM  SERVICE-LEVELS 

For  the  six  experiments  described  above,  we  also  examined  the 
single-item  service-levels  of  the  policies  produced  by  the  GKD  Algorithm. 
For  each  inventory  system  we  chose  that  experiment  that  had  the  average 
service-level  closest  to  85%  and  plotted  the  items  (in  decreasing  order 
of  their  annual  dollar  demand,  hy)  against  their  optimal  service-level. 
The  results  are  displayed  in  Figures  9,  10,  11,  12,  13  and  14.  We  make 
the  following  observations  concerning  them. 

1.  Typically  the  items  in  the  inventory  system  are  skewed,  in. 
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that  the  seven  high-value  items  (i=l,...,7)  have  a  service-level 
considerably  below  85%,  and  most  of  the  twenty-five  low-value  items 
have  a  service-level  considerably  above  85%.  In  some  cases  there 
are  a  few  items  with  service-levels  near  85%. 

2.  There  is  an  obvious  negative  correlation  between  the  value 

h.ij-  and  the  service-level  for  item  i.  This  indicates  that  low-value 

i  i 

items  should  be  stocked  high  and  high-value  items  should  be  stocked 
much  lower,  rather  than  stocking  all  items  the  same. 

3.  There  is  an  obvious  cyclical  pattern  corresponding  to  the 
cyclical  pattern  of  the  means.  For  the  high-value  items,  the  service- 
levels  are  higher  for  those  items  with  smaller  means,  and  for  the  low- 
value  items,  the  service-levels  are  lower  for  those  items  with  the 
smaller  means.  This  suggests  that  it  may  be  possible  to  wel 1 -approxi - 
mate  the  ordering  of  the  items  by  their  service-levels  using  only  h^ 
and  . 


4.  ITEMS  WITH  D.  TOO  SMALL 

The  discussion  in  Section  5  of  the  previous  chapter  indicates 
that  a  sufficient,  but  not  necessary,  condition  that  the  GKD  Algorithm 
converge  to  a  (possibly  local)  minimum  (17)  is  that  £.<0  and  n^>0. 

These  are  guaranteed  if  is  sufficiently  large.  If  is  not  large 
enough  to  guarantee  that  ^<0,  it  is  sufficient  to  in  effect  add  the 
additional  constraint  S^iS.*"  in  order  to  still  guarantee  that  the 
algorithm  will  converge  to  a  local  minimum.  If  ri^>0,  no 
recommendation  could  be  given  to  guarantee  algorithm  convergence  to 
a  local  minimum.  In  this  chapter  we  report  on  those  items  in  the  six 


64 


inventory  systems  examined  that  violate  either  ^<0  or  r^>0. 

Interestingly  enough,  there  is  not  a  single  item  in  any  inventory 
system  with  n^<0,  and  so  the  GKD  Algorithm  converged  to  a  (possibly 
local)  minimum  in  every  experiment.  Each  system  contained  items  with 
£^>0,  however,  necessitating  adding  constraints  ,  as  discussed 

in  the  previous  chapter.  This  is  reported  in  Table  2. 
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5.  CONCLUSIONS 

We  have  found  that  for  inventory  systems  typical  of  ones  found 
in  the  real-world,  there  is  a  significant  cost  savings  when  the  GKD 
Algorithm  is  used  rather  than  the  Identical  Service  Approach  often 
used  by  managers.  There  are  two  issues  that  still  need  to  be  addressed. 

First,  the  behavior  observed  in  Figure  8  (the  identical  service 
approach  performed  better  than  the  GKD  Algorithm)  needs  to  be  better 
understood.  This  is  especially  necessary  in  order  that  practical 
recommendations  may  be  given  as  to  what  types  of  inventory  systems  one 
can  expect  substantial  cost  reductions  when  using  the  GKD  Algorithm  rather 
than  the  identical -service  approach.  Moreover,  insights  into  the  behavior 
observed  in  Figure  8  may  suggest  modifications  to  the  GKD  Algorithm 
that  would  enable  it  to  perform  better  in  situations  like  those  which 
prevail  in  Figure  8.  This  issue  is  addressed  in  the  next  chapter. 

Second,  methods  of  large-scale  inventory  aggregation  need  to  be 
explored  so  that  the  GKD  Algorithm  can  be  used  in  an  inventory  system 
with  perhaps  many  thousands  of  items.  Only  then  can  the  algorithm  be 
of  practical  value.  This  issue  is  also  addressed  in  the  next  chapter. 


66 


INI 

16 

50.0 


12 

37.5 


25. 

1  tem/curnulatf ve  percent  of 


Ej  percent 
Figure  "2" 


total  inventory 


20  .  .  . 
62.5  ;  ; ] 

■  -r  h  +  i  I  i  fI 


24 

\15.Q  ] 

i  a  j 

t  *  - 


32  | 
100.0 


12. 


Uoo 


W- 

i  tem 

cumulate  ve 

|  I  “  - 


67 


68 


69 


70 


optimal  GKD  service-level 


73 


1* 


G£> 


-•—4- 


4  111 


12 

■4—4  >  4 


16 

H — i — I — ► 


20 


24 

"»-!  i 


23 


32 


..95 ; 


.90 


85 


-*—i  t  t 


80 


75 


.70 


65 


60 


55 


50 


45 


40 


-1 — t- 


4  +  -+ 


4—1  -+H 


■+“!*- 4- 


i  : 


K.-24, 


o. 

rH~ 


v 


Dj  usi 


K. 


P..U 


'■  I 


12 


16 


-4 — - 


20 


“4 — t  t  4 
24 


■4— 4  — r— 4* 


tem  (decreasing  order  by1  hp)  ;-j 


3? 


,  Av 


:rage  Servic 


16.  8 


4. 


ig  Power  Metpod 

!  T  1 


'  *  1  1  r  *  1  1  ]  !  I  '  ’ 


L-Leve1=.856,  k^=4,' 


2.....1 


!  '  ‘ 


5.  8,  4,  2  . 

.  :  -  j  ;!  .  : 


1  ’ 


Figure  9 ;  :  :  ; ; 


32 


74 


optimal  GKD  service-level 


75 


1.00 


-l — *- 


12 

i->-»  * 


16 


■*  -*  <  » 


20 

■*—»  *  * 


24 

<  4  »  i 


2$ 


32 

*  ♦  *• 


.95 


.90 


.85 


-.80 


.75 


.70 


.65 


.60 


.55 

.50 

.45 

.40 


1  4  — + — I — * — t  -*■ 

4 


°i 


24, 


W 


h 


V 


1 

3? 


+  — t — h- 
12 


“4  ^4— +— t- 


+  »  t 


16 


20 


— H — t— 1 — H 

24 


,  Av 


item  (decreasing  order  by  hyj 
*  '.852,  k>4, 


2rage  Servic 


£-Level 


16,  8 


4,  .2,.,. 


.15 


8,  4 


.  l_! 


1  l-jt 


32 


Di  usi 


P9 


Power  Met 

I  i  i 


tod  : : ; : 


Figure ;  1,1 


;  :  t  • 


optimal  GKD  service-level 


76 


77 


1.00 

v95 


.90 


-» — *- 


-I — I - I - H 


4  • 


12 

4  4  k  ♦ 


♦  ♦  • 


4  ft 


16 

■4  ♦  <  < 


-20 

■4—4 — I — h 


24 

— t  <  »  I 


23 


♦  •  • 


•  • 


32 


it  k  4 


.85 


.80 


(D 

> 

CD 

I 

QJ 

U 

> 

Sw 

CD 

</> 

Q 

v: 

o 


05 

E 

4~> 

CL 

o 


75 


70 


.65 

.60 


.55 


.50 


.45 


»  • 


.40 


— rr 


T“T 


l  * 

'1— : 


4 — 4 


t  f 
t 


•!! 

i  i 


— 4  — t — f — 4- 

4 


4-  +=4 


-+  *4-4— |— r 


=24, 


W, 


nr 


D-  usi 

.  1  f  .  t 


12 


-4—4—4-+- 


=  32.  Avp 


16 

item  (decreas 

rage  Service-level  * k.=4,  |  ! 


20 

'ing  order  by 


9*  u . :  16,.. 

:>!:  M  ’  ■  : 

|  i  i  ;  ’  i  :  '  • 

ng  Power  Me t|i 

mm; 


,16,  B,..., 
od 


B'  : 


1  1  !  t  t  :  ’  •  t 


t  I  " 
Mill' 


:  !  !  I 

MM 


-4—1-4 

24  . 

Wl- 

iii  t* 


+  +  +  +  I  + 


28 


Figur 


*  t't  :  1 


;  i  i 


i  •? 


!  I!  I  '  t 


4 


!  I 


e  13.  !  ;  | : 

r  7'M  *  t 


32 


optimal  GKD  service-level 


78 


79 


<Sj 

o 


JO 

m 


Cvj 

co 


N 

r~> 


-ft 


IV.  SENSITIVITY  EXPERIENCE  WITH  THE  GKD  ALGORITHM 


1.  INTRODUCTION 

In  Chapter  II  we  considered  the  problem  of  specifying  single-item 
service  objectives  in  a  multi-item  (s,S)  inventory  system  subject  to  a 
system-wide  service  objective.  We  formulated  this  problem  as  a 
constrained  NLP  and  developed  the  Generalized  Knapsack  Duality  (GKD) 
Algorithm  to  solve  it.  In  Chapter  III  we  reported  our  computational 
experience  with  several  32-item  inventory  systems  with  a  structure  typical 
of  many  real-world  inventory  systems.  In  this  chapter  we  report  sensitiv¬ 
ity  experience  with  the  base-case  and  worst-case  inventory  systems  used 
in  Chapter  III.  The  criterion  of  comparison  is  to  compare  the  expected 
holding  cost  (computed  exactly)  associated  with  the  (s,S)  policies 
produced  by  the  GKD  Algorithm  for  a  given  system-wide  service  level  a 
with  the  expected  holding  cost  of  the  policies  given  by  the  Identical 
Service  Method,  in  which  all  the  (s,S)  policies  are  set  by  the  Power 
Approximation  given  in  Mosier  [1981]  so  that  every  item  has  a  service- 
level  a  (as  discussed  in  Chapter  III  the  exact  algorithm  of  Veinott  and 
Wagner  [1965]  is  computationally  unwieldy  for  our  inventory  system, 
and  the  Power  Approximation  is  the  most  accurate  approximation  available). 
In  computing  the  expected  costs  we  assume  the  underlying  demand  distri¬ 
bution  is  Negative  Binomial,  and  we  are  primarily  interested  in  the  sensi¬ 
tivity  of  the  GKD  Algorithm  at  service-levels  of  85£and88%.  We 
are  interested  in  85%  because  this  is  a  service-level  managers  often 
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prefer,  and  we  are  interested  in  881  because  of  the  difficulties 
reported  in  Chapter  II I  with  the  worst-case  inventory  system  at  this 
servi ce-level . 

In  Section  2  of  this  chapter  we  report  sensitivity  experience 
for  D^=S.j-s.j,  which  is  a  constant  input  parameter  to  the  GKD  Algorithm. 
In  Section  3  we  report  sensitivity  experience  for  the  lower  bounds  on 
S.j ,  which  are  input  to  the  GKD  Algorithm,  but  may  be  increased  by  the 
GKD  Algorithm  (Chapter  II  Section  5).  In  Section  4  we  consider  several 
sampling  schemes  for  large-scale  inventory  aggregation,  where  a  sample 
is  used  with  the  GKD  Algorithm  to  specify  (s,S) f pol i c ies  for  every  item 
in  the  inventory  system.  In  Section  5  we  investigate  the  sensitivity 
of  these  sampling  schemes  to  the  number  of  items  in  the  scheme,  and  in 
Section  6  we  offer  conclusions  and  recommendations  concerning  specifying 
the  D.j ,  the  lower  bounds  on  S^  and  the  size  and  type  of  inventory  scheme 
to  use  for  large-scale  inventory  aggregation.  All  figures  are  at  the  end 
of  the  chapter. 

2.  SENSITIVITY  EXPERIENCE  FOR  D. 

Recall  from  Section  2  of  Chapterll  that  D^=S--s^  is  approximated 
by  a  constant  (the  Power  Approximation  of  D. )  which  is  input  to  the 
GKD  Algorithm.  This  is  justified  because  in  the  single-item  case  the 
optimal  D^  is  asymptotically  independent  of  the  service-level.  For  any 
given  item  i  in  the  inventory  system,  however,  the  policy  (s^,S^) 
produced  by  the  GKD  Algorithm  with  associated  service-level  may  not 
be  an  optimal  policy  when  the  item  i  is  considered  as  a  single-item 


82 


inventory  system.  And  for  a  non-optima 1  policy  (s^.S^),  is  not 
necessarily  independent  of  the  service-level  a... 

Thus  we  investigated  the  impact  of  varying  the  input  values  of 
the  D.j ,  to  determine  if  the  GKD  Algorithm  performance  degraded  (or 
perhaps  even  improved).  As  stated  in  the  Introduction,  we  used  the 
32-item  base-case  inventory  system  of  Chapter  III.  This  system  has  the 
structure  that  twenty  percent  of  the  items  constitute  eighty  percent 
of  the  total  value  (where  value  is  taken  to  be  hy)  of  the 
inventory  system.  The  parameters  are  as  follows  for  1*1,..., N: 


K.  =  24 


yi 

v*l  •  •  • 

■,  p29  = 

*  16 

Who**’ 

*’  u30 

=  8 

V3.p7.pn... 

...  P31 

=  4 

P4.P8.P12’- • 

■ P32 

=  2. 

The  h^  are  determined  through  the  value  h..y..  Additional 
details  are  given  in  Chapter  III. 

Because  of  the  structure  of  the  inventory  system,  the  items  1-25  are 
low-value  items  and  items  26-32  are  high-value  items.  Typically,  as 
discussed  in  Section  3  of  Chapter  III,  the  GKD  Algorithm  produces  policies 
so  that  the  service-levels  of  the  high-value  items  are  very  low  and  the 
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service-level s  of  the  low-value  items  are  very  high.  These  observations 
suggest  two  experiments. 

In  the  first,  we  increased  the  value  of  (given  by  the  Power 
Approximation)  25%  for  the  low-value  items  and  decreased  the  value  of 
25,.  for  the  high-value  items.  In  the  second  experiment,  we  decreased 
the  value  of  D^.  by  25%  for  the  low- value  items  and  increased  it  25%  for 
the  high-value  items. 

These  are  severe  tests  of  robustness.  In  the  experiments  in  Chapter 
III  we  compared  the  expected  holding  costs  of  the  GKD  Algorithm  policies 
and  the  Power  Approximation  policies  (which  give  identical  service  to 
all  items)  because,  since  we  used  the  Power  Approximation  Di  in  the  GKD 
Algorithm,  the  expected  replenishment  costs  are  the  same.  Such  is  not  the 
case  in  the  experiments  described  in  this  section,  so  we  compare  the  total 
expteced  costs,  namely,  the  expected  holding  plus  replenishment  costs. 

These  experiments  and  their  outcomes  are  summarized  in  Table  3  and  Figures 
15,  16  and  17.  In  Experiment  2,  for  example,  in  which  the  Power 
Approximation  D..  is  increased  25%  for  the  low-value  items  and  decreased  25% 
for  the  high-value  items,  the  expected  total  cost  for  the  GKD  policies 
at  a  service-level  of  85%  is  approximately  198,  and  the  expected  total 
costs  for  the  Identical  Service  Method  at  that  same  service-level  is 
approximately  276.  This  represents  a  28%  decrease  in  expected  cost  when 
using  the  GKD  policies  rather  than  the  Identical  Service  policies. 
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Table  3.  Sensitivity  Experience  for  dy.  Comparison  of  GKD  Algorithm 
Policies  with  Identical  Service  Policies. 


Experiment 

Low-Value 
Items  (1-25) 

High-Val ue 
Items  (26-32) 

%  4  in  cost 
at  Service- 
Level  85% 

%  4  in  cost 
at  Service- 
Level  88% 

See 

Fi gure 

1 

No  change 

No  change 

29 

26 

15 

2 

D.  t  25% 

Di  4  25% 

26 

22 

16 

3 

D.  4  25% 

l 

D.  t  25% 

28 

27 

17 

We  make  the  following  observations  and  conclusions: 

(1)  The  change  in  the  performance  of  the  GKD  Algorithm 
in  Experiment  2  and  3  is  slight,  especially  considering 

the  severity  of  the  test.  This  suggests  that  not  only  is  an 
especially  accurate  specification  of  D  unnecessary  for  a  single¬ 
item  inventory  system,  it  is  also  unnecessary  for  a  multi-item 
inventory  system  of  the  general  structure  of  our  base-case.  Of 
course,  additional  experiments  on  variations  of  this  base-case 
would  help  confirm  or  refute  this  conjecture. 

(2)  Although  it  may  be  possible  to  specify  the  so  that  greater 
cost  savings  are  realized  than  when  using  the  Power  Approximation 
D.j ,  we  recommend  using  the  Power  Approximation,  since  obvious  ways 
of  changing  do  not  improve  the  GKD  Algorithm  performance. 
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3.  SENSITIVITY  EXPERIENCE  FOR  LOWER  BOUNDS  ON  Si 
Recall  from  Chapter  II  that  we  required  S^D.  in  the  GKD  Algorithm 
with  that  lower  bound  on  S^  possibly  being  increased  by  the  algorithm 
to  guarantee  convergence  to  a  (possibly  local)  minimum  (see  Section  5). 

The  main  reason  given  to  restrict  attention  to  policies  with  S^>Di  is 
that  (aside  from  the  mathematical  difficulties  involved)  many  managers 
consider  policies  with  (equivalently,  s^<0)  unsatisfactory  for 

practical  implementation.  There  are  two  reasons,  however,  that  even 
higher  lower  bounds  on  S^  might  be  imposed. 

First,  many  managers  prefer  not  only  that  s^>0,  but  that  s^>p-*  = 
OVUVr  so  the  safety  stock  (s^  is  always  at  least  large  enough  to 
cover  the  expected  leadtime  demand. 

Second,  in  Chapter  III  the  Identical  Service  Approach  performed  better 
than  the  GKD  Algorithm  in  the  worst-case  experiment  at  a  service-level 
at  and  near  88%. 

We  offer  a  conjecture  to  explain  this  behavior  and  show  that  by 
modifying  the  lower  bounds  on  S^  the  GKD  Algorithm  will  in  the  case 
just  described  perform  considerably  better  than  the  Identical  Service 
Approach  at  both  service  levels  85%  and  88%. 

The  worst-case  inventory  system  has  the  same  80/20  structure 
(twenty  percent  of  the  items  represent  eighty  percent  of  the  total  system 
value)  with  the  following  parameters  for  i=l,...,N: 

K.  =  24 
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2 


and 

W'j *1^2 •  •  •  • »  = 

y17*y18 . y32  =  8' 

There  is  a  notable  difference  in  the  single-item  service  levels 
determined  by  the  GKD  Algorithm  at  a  service-level  of  only  85\  for  both 
the  base-case  and  the  worst-case  inventory  systems,  as  reported  in 
Figures  9  and  14  of  Chapter  1 1 1, respectively.  For  the  high-value  items 
(items  1-7)  in  both  cases  the  (s,S)  policies  given  by  the  GKD  Algorithm 
have  all  the  at  their  lower  bounds.  The  service-levels  associated 
with  the  high-value  items  for  the  worst-case  are  considerably  lower  than 
those  for  the  base-case,  and  so  the  service-levels  for  the  low-value 
items  for  the  worst-case  are  considerably  higher  than  those  for  the 
base-case.  The  GKD  Algorithm  assumes  the  demand  distribution  is  Logistic, 
whereas  in  these  experiments  we  assume  the  demand  distribution  is  Negative 
Binomial.  The  approximation  of  a  Negative  Binomial  distribution  by  a 
Logistic  distribution  is  poorest  in  the  tail,  which  corresponds  to  very 
high  individual  item  service-levels.  Thus  we  conjecture  that  the  poor 
performance  of  the  GKD  Algorithm  in  the  worst-case  is  due  to  the  poorness 
of  the  Logistic  demand  approximation.  We  test  this  conjecture  by 
imposing  higher  values  of  the  lower  bounds  on  the  .  This  forces  the 
high-value  items  to  be  stocked  at  a  higher  service-level,  so  the  low- 
value  items  are  stocked  at  a  lower  service  level,  as  in  the  base-case. 

We  note  that  Figures  9  and  14  indicate  that  the  higher  the  demand  mean, 
the  lower  the  service- level  of  the  (s,S)  policy  with  S  at  its  lower  bound. 
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perform  as  well  at  a  service-level  of  88jt  as  it  does  at  91)6  with  r=l. 
These  experiments  and  their  outcomes  for  the  base-case  are  summarized 
in  Table  6.  The  results  are  similar  to  those  in  Experiments  2,  3,  and 
4,  and  so  we  do  not  report  the  expected  cost  graphs. 


Table__6.  Sensitivity  Experience  for  Lower  Bounds  on  Si:  Base-Case; 

Comparison  of  GKD  Algorithm  Policies  with  Identical  Service 
Pol icies. 


Experiment 

r,  where 

vw 

Minimum 

Overal 1 
Service-Level 

%  4  in  cost 
at  Service- 
Level  88?; 

6 

.8 

.860 

33 

7 

.9 

.885 

— 

These  experiments  and  their  outcomes  for  the  worst-case  are 
summarized  in  Table  7  and  Figure  22.  Since  theresults  of  these  experi¬ 
ments  are  similar  to  each  other  we  report  the  expected  cost  graph  only 
for  Experiment  6. 


91 


Table  7.  Sensitivity  Experience  for  Lower  Bounds  on  S • :  Worst-Casei 
Comparison  of  GKD  Algorithm  Policies  with  Identical  Service 
Policies. 


Minimum 

%  4-  in  cost  %  4-  in  cost 

r,  where 

Overall 

at  Service-  at  Service- 

See 

Experiment 

SjfcD.+rpj* 

Service-Level 

Level  85% 

Level  88% 

Figure 

6 

.8 

.811 

26 

23 

22 

7 

.9 

.878 

20 

-- 

We  make  two  additional  observations: 

(5)  The  previous  remarks  for  the  base-case  also  hold  for 
r=.8,  .9. 

(6)  For  the  worst-case  system  there  is  a  dramatic  improvement 
in  algorithm  performance  when  r=.8.  In  this  case  there  is  a 
substantial  cost  savings  at  both  service-levels  85%  and  88%. 

We  offer  the  following  conclusions. 

(1)  The  behavior  at  a  service-level  of  88%  for  the  worst-case 
inventory  system  does  seem  to  be  caused  (at  least  in  part)  by  the 
poorness  of  approximating  a  Negative  Binomial  distribution  by  a 
Logistic  distribution. 

(2)  The  algorithm  performs  best  when  there  is  a  lower  bound  on  S 
just  large  enough  to  prevent  too  many  of  the  low-value  items  to  be 
stocked  extremely  high.  We  recommend  that  a  user  experiment  with 
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various  lower  bounds  (in  terms  of  the  single  input  value  r)  to 
find  the  best  one. 

(3)  As  most  managers  prefer  as  high  a  lower  bound  as  possible  on 
S  (as  long  as  average  system  service  is  not  too  large),  we  reco¬ 
mmend  specifying  r  as  large  as  possible  to  include  the  desired 
targeted  service-level . Our  experience  with  both  the  best  and  worst- 
case  inventory  systems  suggests  that  this  approach  will  yield  a 
significant  cost  reduction  over  the  identical  service  approach. 

4.  SENSITIVITY  EXPERIENCE  FOR  SAMPLING  SCHEMES  FOR 
INVENTORY  AGGREGATION 

Many  real-world  inventory  systems  contain  tens  of  thousands  of 
items.  Although  the  GKD  Algorithm  is  quite  efficient,  it  could  not 
handle  inventory  systems  of  this  size.  This  is  not  a  problem  unique  to 
our  algorithm,  and  most  managers  deal  with  it  using  inventory  aggregation. 
This  is  any  method  of  partitioning  an  inventory  system  into  groups  of 
items  so  that  every  item  in  the  group  is  controlled  essentially  the 
same.  A  representative  item  is  chosen  from  each  group  and  used  to 
specify  (s,S)  policies  for  all  items  in  its  group. 

Our  method  of  aggregation  exploits  the  special  structure  of  our 
problem  (as  formulated  as  an  NLP)  and  of  the  GKD  Algorithm,  and  certain 
empirical  observations  in  Chapter  III.  Recall  from  Section  3  of  Chapter  II 
that  the  GKD  Algorithm  exploits  the  "Knapsack"  structure  of  the  NLP  by 
solving  its  dual  NLP.  The  dual  has  only  one  variable,  X,  and  to  find 
the  optimal  S*  of  the  original  NLP  it  is  sufficient  to  find  the  optimal 
dual  variable  X*.  This  is  reflected  in  the  GKD  Algorithm  (Chapter  II 
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Section  3)  by  solving  P(x)=0  for  X*,  from  which  S*  is  computed.  Our 
method  of  aggregation  is  to  use  a  subset  of  the  inventory  system  to 
estimate  X*,  and  then  compute  S*  for  all  the  items  in  the  inventory 
system  from  this  estimate. 

The  inventory  system  we  constructed  to  test  various  aggregation 
schemes  is  a  128-item  system  with  the  same  essential  structure  as  the 
32-item  base-case  system  described  in  Section  2.  In  particular,  twenty 
percent  of  the  items  (items  101-128)  constitute  eighty  percent  of  the 
total  value  of  the  inventory  system,  and  the  parameters  are  as  follows, 
for  i  =  l ,  2, .  . . ,  128: 


Ki  =  24 


and 

W***’  u125  = 
p2,u6,...,  U] 26  = 
u3,u7,. ,  Pi 27 
^4  ,'J8*  *  *  *  *  ^  1 28  = 

To  examine  the  cost  savings  in  using  the  GKD  Algorithm  rather  than 
the  Identical  Service  Method  (via  the  Power  Approximation)  described  in 
Section  1,  we  ran  the  GKD  Algorithm  with  the  entire  system.  As  described 
in  Section  1,  we  compared  the  expected  holding  costs  with  those  of  the 
Identical  Service  Method,  and  the  result  is  displayed  in  Figure  23. 


16 

8 

4 

2. 
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divided  by  the  total  value  of  the  sample  from  that  strata.  This  was  done  for 

both  strata,  so  each  maintained  its  total  value,  thus  preserving  the 
80/20  structure.  And  to  maintain  the  given  constraint  weights  W.. ,  they 
were  each  recomputed  as  the  total  weight  of  the  strata  divided  by  the 
total  weight  of  the  sample  from  that  strata. 

Initially  we  tested  three  central  nonrandom  stratified  uniform, 
sampling  schemes.  If  one  wants  to  choose  M  items  from  a  stratum  consist¬ 
ing  of  the  N  items  1,  2,...,  N  by  such  a  scheme,  then  the  items 

N  3N  5N  ( 2K-1 )N 

2K’  2M’  2M .  2M 

are  taken,  where  the  fractions  ^  are  rounded  to  the  nearest  integer. 

For  this  sample,  the  80/20  value  structure  and  the  constraint-weight  structure 

is  maintained  by  multiplying  tv  and  1C  by  .£.  Ivy./  .1.  h^,  and  ^-j 

system  sample 

by  N/M  for  each  item  in  the  sample.  The  three  schemes  differed  in  the 
number  of  items  taken  from  each  of  the  two  strata.  They  are  described 
in  Table  8.  Notice  that  Scheme  1  is  a  nonstratified  scheme  that  roughly 
takes  every  4th  item  from  the  inventory  system  (since  the  means  repeat 
themselves  every  four  items,  to  take  every  4th  item  would  create  an 
inventory  system  with  every  mean  the  same,  which  seemed  objectionable). 
Scheme  2  creates  an  inventory  subsystem  of  half  high-value  and  half  low- 
value  items,  and  that  Scheme  3  takes  all  the  high-value  items  (it  seems 
reasonable  that  A*  is  more  sensitive  to  them  than  to  the  very  low  value 
items,  and  so  this  scheme  seemed  plausible). 
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Table  8.  Central  Nonrandom  Stratified  Uniform  Sampling. 


Scheme 

Number  Chosen  from 
Low-Value  (1-100) 

Number  Chosen  from 
High-Value  (101-128) 

1 

24 

8 

2 

16 

16 

3 

4 

28 

The  experiment  we  performed  to  test  these  schemes  is  as  follows: 

For  a  variety  of  specified  service-levels  we  ran  the  GKO  Algorithm  using 
the  inventory  subsystem  (determined  by  the  sampling  scheme)  to  compute 


\*,  the  optimal  multiplier  for  the  subsystem.  Each  value  of  X*  gives 
rise  to  a  service-level  for  the  subsystem  (0(32)  and  for  the  entire 
system  (ai 33 ^ ’  and  t0  the  exPected  holding  cost  at  these  service-levels 
(H^2  and  H^g,  respectively).  The  criterion  of  comparison  is  how  close 

A 

Q^anda^g  and  Hg.,  and  H^g  (for  a  given  X)  are.  A  "perfect"  subsystem 
would  give  rise  *°  a32=Qi28*  ^32=^128’  ^s>lts  °Ptimal  multiplier 

A 

X*=>.*,  the  optimal  multiplier  for  the  entire  system.  The  service-level 
graphs  for  the  three  schemes  are  given  in  Figures  24,  25  and  26,  and  the 
expected  holding  cost  graphs  are  given  in  Figures  27,  28  and  29. 

We  make  the  following  observations  concerning  these  sampling  schemes 
(1)  They  are  all  good.  Although  each  one  gives  rise  to  a  misspeci- 
fication  of  a^g  and  H^g,  the  total  cost  savings  will  still  be 


97 


substantial  over  the  Identical  Service  Approach  (compare  Figure  23). 
(2)  Scheme  2  seems  to  be  superior  to  Schemes  1  and  3.  Scheme  1 
has  a  preponderance  of  low-value  items.  Scheme  3  has  a  preponder¬ 
ance  of  high-value  items,  and  Scheme  2  has  the  same  number  of  each. 

It  seems  that  the  proportion  of  items  that  are  high-value/low-value 
at  least  in  part  affects  the  performance  of  the  scheme. 

The  efficiency  of  Scheme  2  may  have  something  to  do  with  the  fact  that 
the  demand  means  of  our  base-case  inventory  system  follow  the  pattern 
2,  4,  8,  16,...,  2,  4,  8,  16.  To  investigate  this,  we  used  a  central 
random  stratified  uniform  sampling  scheme  with  16  expensive  and  16  inex¬ 
pensive  items  (like  Scheme  2).  In  other  words,  we  chose  16  items  uniformly 
randomly  (using  the  IMSL  Routine  GGUD)  from  both  the  high-value  and  low- 
value  items.  We  performed  six  such  tests,  one  with  a  different  initial 
seed.  The  details  of  the  subsystems  chosen  can  be  found  in  the  Appendix 
to  this  chapter,  and  the  scheme  itself.  Scheme  4,  is  described  in  Table  9. 


Table  9. 

Central  Random  Stratified 

Uniform  Sampling 

Scheme 

Number  Chosen  from 
Low-Value  (1-100) 

Number  Chosen  from 
High-Value  (101-128) 

4 

16 

16 

We  performed  the  same  experiment  as  we  did  for  the  first  three 
schemes,  and  the  results  of  this  experiment  for  the  first  seed  is 
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reported  in  Figures  30  and  31  (results  for  all  of  the  seeds  were  similar). 

We  make  the  following  observations. 

(1)  As  for  the  three  nonrandom  schemes,  the  random  scheme  is  good. 

(2)  The  random  scheme  is  not  as  good  as  Scheme  2,  the  nonrandom 
scheme  with  the  same  proportion  of  high-value/low-value  items  as 
this  scheme.  It  would  appear  that  the  sequence  of  demand  means  is 
important,  and  that  further  research  in  stratifying  not  only  on  the 
value  but  on  the  demand  means  could  give  even  better  sampling 
schemes. 

(3)  The  method  of  sampling  seems  relatively  unimportant  as  long  as 
the  80/20  value  structure  of  the  inventory  system  is  maintained. 

We  recommend,  therefore,  that  the  scheme  be  chosen  on  the  basis  of 
its  ease  in  implementation.  Perhaps  the  nonrandom  Scheme  1  is 
easiest,  as  it  requires  only  one  stratum  (it  represents  a  central 
sampling  of  every  4th  item). 

As  a  final  experiment  we  tested  whether  Scheme  4  with  the  seed  of  the 
first  experiment  was  accurate  when  lower  bounds  are  imposed  on  the  Sy  We 
chose  Scheme  4  because  it  gave  rise  to  neither  the  best  nor  the  worst 
performance  among  the  schemes.  We  chose  lower  bounds  corresponding  to 
r=.75  (see  Section  3)  because  this  was  the  largest  value  of  r  which  allowed 
an  overall  service-level  of  at  least  85%.  This  experiment  is  described 
in  Table  10  and  the  outcome  is  reported  in  Figures  32,  33,  and  34,  where 
Fiqure  32  reports  the  expected  cost  when  using  all  128  items  in  the  GKD 
Algorithm. 
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Table  10. 

Central  Random  Stratified 
S.*D.+.75p.*. 

Uniform  Sampling: 

Scheme 

Number  Chosen  from 
Low-Value  (1-100) 

Number  Chosen  from 
High-Value  (101-128) 

4 

16 

16 

We  make  the  following  observations: 

(1)  There  is  a  significant  cost  savings  at  service-levels  852; 
and  88%. 

(2)  The  sampling  scheme  is  very  accurate.  It  seems  that  our 
previous  conclusions  about  the  accuracy  of  these  sampling  schemes 
also  hold  true  when  lower  bounds  are  imposed  on  the  . 


5.  SENSITIVITY  EXPERIENCE  FOR  THE  SAMPLE  SIZE 
FOR  INVENTORY  AGGREGATION 

In  order  to  investigate  how  the  reliability  of  inventory  sampling 
schemes  depends  on  the  size  of  the  scheme,  we  constructed  a  512-item 
base-case  inventory  system  with  the  80/20  value  structure  described  in 
Section  2.  Not  only  is  this  structure  typical  of  many  real-world  inven¬ 
tory  systems,  but  if  the  items  in  the  system  are  ordered  in  descending 
order  by  their  value  (h^.)  and  the  cumulative  percent  of  the  items  is 
plotted  against  the  corresponding  cumulative  value,  then  the  resulting 
graph  is  frequently  wel 1 -approximated  by  a  Lognormal  cumulative  distribu¬ 
tion  function  (Peterson  and  Silver  [1979,  pp.  30-37]).  We  constructed 


101 


(Nonrandom)  Uniform  Sampling  Scheme  (compare  Scheme  1  of  the  previous 
section).  The  last  experiment  (Scheme  4)  of  Section  4  suggests  that 
the  approach  of  using  a  sample  to  investigate  the  sensitivity  of  the 
system  to  r  is  a  reliable  procedure.  We  considered  r=0,  .6,  .7,  .8, 
.9,  1,  with  a  targeted  service-level  of  85%  (a  service-level  popular 
with  many  managers),  and  we  want  to  include  for  consideration  service- 
levels  83%  to  88%. 

This  experiment  and  its  result  is  reported  in  Table  11  and  Figure 
35.  Since  we  want  to  include  for  consideration  a  service-level  of  83* 
we  chose  r=.6.  The  cost  savings  forr  =  .6  at  a  service-level  of  85*  is 
about  51 r ,  which  is  very  signficant. 


Table  11 .  128-Item  Central  Nonrandom  Uniform  Sample  from  the  512- 
Item  Base-Case  System:  Comparison  of  GKD  Algorithm 
Policies  with  Identical  Service  Policies. 


where 

w 

Minimum 

Overal 1 
Service-Level 

%  4-  in  cost 
at  Service- 
Level  85« 

0 

<  .8100 

50 

.6 

.8140 

51 

.7 

.8439 

47 

.8 

.8719 

-- 

.9 

.8951 

-- 

1 

.9163 
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As  an  additional  check  of  the  reliability  of  using  this  sampling 
method  to  chose  r,  we  repeated  the  experiment  on  the  entire  512-item 
inventory  system.  The  results  are  summarised  in  Table  12  and  Figure  36. 
The  sampling  method  of  choosing  r  is  seen  to  be  very  reliable,  as  was 
the  case  in  Section  4. 


Table  12.  512-Item 
Pol icies 

Base-Case  System:  Comparison  of 
with  Identical  Service  Policies. 

GKD  Algorithm 

Minimum 

%  4-  in  cost 

r,  where 

Overal 1 

at  Service- 

Si>Di+rpi* 

Service-Level 

Level  85% 

0 

<  .8100 

49 

.6 

.8138 

50 

.7 

.8441 

46 

.8 

.8720 

— 

.9 

.8961 

— 

1 

.9165 

-- 

There  were  two  sampling  schemes  considered  in  Section  4  which  seemed  to 


perform  the  best.  The  first  is  the  recommended  scheme,  the  Central 
Nonrandom  Uniform  Sampling  Scheme.  The  second  is  the  Central  Nonrandom 
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Stratified  Uniform  Sampling  Scheme  (compare  Scheme  2)  in  which  half  the 
sample  is  taken  from  the  high-value  items  and  the  other  half  is  taken 
from  the  low-value  items.  Because  of  the  80/20  value  structure  of  our 
512-item  base-case  system,  items  1-102  are  high-value  and  items  103-512 
are  low-value.  The  sensitivity  experimental  design  is  identical  to 
that  described  in  Section  4,  here  with  r=.6.  We  chose  sample  sizes  of 
32,  64  and  128  for  both  schemes  described,  and  the  six  experiments  are 
summarized  in  Table  13.  For  each  sample  size,  the  results  of  these  exper¬ 
iments  were  similar  for  both  the  uniform  and  stratified  uniform  schemes, 
and  so  the  comparative  graphs  are  reported  only  for  the  Central  Nonrandom 
Uniform  Sampling  Scheme  (Figures  37  through  42). 

Table  13.  Comparison  of  GKD  Algorithm  Policies  with  Identical  Service 
Policies. 


Central  Nonrandom  Uniform  Sampling  from  51 2- 1  tern  Base-Case  System, 
S.>D.+.6y.*. 


Experiment 

Sample 

Size 

Figure  Reporting 
Service-Levels 

Figure  Reporting 
Expected  Costs 

1 

32 

37 

40 

2 

64 

38 

41 

3 

128 

39 

42 

Central  Nonrandom  Stratified  Uniform  Sampling 
System,  S.2D^+. 6y.*. 

from  512-Item  Base-Case 

Experiment 

Number  from 
High-Value  Items 

Number  from 
Low-Value  Items 

4 

16 

16 

5 

32 

32 

6 


64 


64 
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We  make  the  following  observations: 

(1)  All  six  schemes  are  reliable  enough  for  practical  implement¬ 
ation.  The  128-item  schemes  are  superior  to  the  64-item  schemes 
which  are  superior  to  the  32-item  schemes. 

(2)  There  does  not  seem  to  be  any  real  difference  between  the  uni¬ 
form  and  stratified  schemes. 

(3)  An  examination  of  the  experiments  in  Section  5.  in  which  32- 
item  samples  were  taken  from  a  128- item  inventory  system,  indicates 
that  the  32-item  samples  are  just  as  reliable  in  the  512-item  system 
as  in  the  128-item  system.  This  is  extremely  significant  because  it 
strongly  suggests  that  32-item  schemes  are  sufficiently  reliable  for 
practical  implementation  no  matter  how  large  the  inventory  system  is, 
as  long  as  it  has  the  80/20  value  structure  described. 

(4)  We  recommend  using  the  Central  Uniform  Sampling  Scheme  with 
around  32  items  for  any  inventory  system  with  an  80/20  value 
structure.  Although  larger  samples  may  be  more  reliable,  this 
scheme  is  reliable  enough.  Moreover,  the  computational  work  necess¬ 
ary  to  use  the  GKD  Algorithm  on  such  a  sample  is  very  small,  making 
this  scheme  suitable  for  frequent  use  in  a  real-world  setting. 

6.  CONCLUSIONS 

In  this  chapter  we  performed  sensitivity  tests  on  the  GKD  Algorithm 
with  inventory  systems  with  a  structure  typical  of  real-world  systems, 
that  twenty  percent  of  the  items  represent  80  percent  of  the  value  of  the 
system.  We  concluded  the  following: 
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(1)  The  total  cost  savings  are  not  affected  very  much  by  changes 
in  . 

(2)  Increasing  the  lower  bounds  on  S..  may  improve  the  algorithm 
performance,  especially  when  the  high-value  items  also  tend  to  have 
high  expected  demands.  Even  when  increasing  these  lower  bounds 
degrades  the  performance,  however,  the  cost  savings  are  substantial. 
We  recommend  increasing  these  lower  bounds  as  high  as  possible 
(while  still  maintaining  a  cost  savings  below  that  of  the  Identical 
Service  Method)  as  this  appraoch  is  most  consistent  with  managerial 
goals. 

(3)  Uniform  sampling  with  about  32  items  is  a  reliable  and  computa¬ 
tionally  efficient  method  of  inventory  aggregation,  giving  accurate 
predictions  via  the  GKD  Algorithm  of  overall  inventory  service  and 
expected  operating  costs. 
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Figure  40.  512  -  Item  Base  -  Case  with  Central  Nonrandom  Uniform 

Sampling  Scheme  (32-Item  Sample) 

S1  >  Di  +  .6}^.* 
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Expected  Holding  Cost:  Subsystem 
Figure  41.  512  -  Item  Base  -  Case  with  Central  Nonrandom  Uniform 
Sampling  Scheme  (64  -  Item  Sample) 
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Figure  42.  512  -  Item  Base-Case  with  Central  Nonrandom  Uniform 
Sampling  Scheme  (128- 1  tern  Sample) 

si  *  Di  +  -6  “i* 
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APPENDIX  TO  CHAPTER  IV 

This  appendix  reports  the  seeds  used  in  the  IMSL  Routine  GGUD 
to  create  the  six  random  samples  of  Scheme  4  of  Section  4,  and  a 
listing  of  the  512-item  inventory  system  of  Section  5.  The  seed 
used  in  the  IMSL  Routine  GGUD  to  generate  the  y..  for  this  system  is 
872245376. 


Table  14. 

Seeds  for  Six  Tests  of  Central  Random 
Stratified  Uniform  Sampling  Scheme  4. 

Test 

Seed 

1 

96665107 

2 

227164307 

3 

929811759 

4 

2058317624 

5 

75813808 

6 

1743717576 

NOTE:  Experiment  for  Test  1  reported  in  Figures  29  and  30. 


51 2- Item  Inventory  System  of  Chapter  IV 


For  all  items  i  =  1,2,. ..,512: 


K.  =  24 


k. 


l 


=  4 


9 


i 

hi 

Di 

Vi 

i 

i.  ei45 

9.0000 

9.0000 

17.5010 

1.  9446 

2 

2. 0363 

8.0000 

8. 4853 

15.7557 

1.9695 

3 

2.  7002 

6.0000 

7. 3485 

12. 1972 

2.0329 

4 

1.2363 

13. 0000 

10. 8167 

24. 6235 

1.  8941 

5 

2.6514 

6. 0000 

7.3485 

12.3104 

2.0517 

6 

1. 2090 

13. 0000 

10. 8167 

24.9025 

1.9156 

7 

0. 9691 

16. 0000 

12.0000 

30.3035 

1. 8940 

8 

3.8186 

4. 0000 

6.0000 

8.7403 

2. 1851 

9 

1.0020 

15. 0000 

11.6190 

29.0211 

1.9347 

10 

1. 8469 

8.0000 

8. 4853 

16.  5535 

2.0692 

1 1 

1.  1163 

13. 0000 

10.8167 

25.9282 

1.9945 

12 

1. C957 

13. 0000 

10. 8167 

26. 1746 

2.0134 

13 

1. 0747 

13.0000 

10.8167 

26. 4318 

2.  0332 

14 

1.  5218 

9. 0000 

9.0000 

19. 1301 

2.  1256 

15 

0. 8388 

16. 0000 

12.0000 

32.6002 

2.0375 

16 

2.6290 

5. 0000 

6.7082 

1 1.  5135 

2.3027 

17 

0.9  193 

14. 0000 

11.2250 

29.4759 

2. 1054 

18 

2.  C997 

6. 0000 

7.3485 

13.8526 

2.3088 

19 

1.  541  1 

8.0000 

8.4853 

18.  1413 

2. 2677 

20 

1.2063 

10.  0000 

9.4868 

22.  4390 

2.2439 

21 

0.  9077 

13.0000 

10. 8167 

28.7900 

2. 2146 

22 

0.9618 

12. 0000 

10.3923 

27.0711 

2. 2559 

23 

2. 2574 

5. 0000 

6.7082 

12.4364 

2.4873 

24 

2.7593 

4. 0000 

6. 0000 

10. 3021 

2.5755 

25 

1.  0792 

10. 0000 

9.4868 

23. 7387 

2.  3739 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

*2 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 
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0. 8793 
0. 8597 
3.  3623 
1. 4088 
0.8035 
4. 7132 
0.7680 
3.  0037 
0. 881  1 
2.8719 
0. 5266 
0.9155 
0. 6715 
0. 4926 
1.9272 
0.6855 
0.5269 
0. 4511 
0.7646 
1.7274 
1. 3523 
2.2058 
6. 4770 
0.7925 
0. 6206 
0. 5064 
0. 3718 
1.1651 
1.  1410 
0. 3492 
0.4975 
0.4874 
1.  7507 
0.6432 
0. 4584 
0. 4492 
4. 6424 
0. 2966 
0. 3877 
0.  2850 
0.7452 
0. 3372 
4. 2988 
0. 4216 
0. 2953 
1.  3518 
0.6630 
1.9513 
0.  7658 
0. 3757 


12. 0000 
12.0000 
3. 0000 
7. 0000 
12.0000 
2. 0000 
12. 0000 
3.0000 
10.0000 
3. 0000 
16.0000 
9. 0000 
12.0000 
16. 0000 
4.0000 
1  1.  0000 
14. 0000 
16. 0000 
9.0000 
4. 0000 
5. 0000 
3. 0000 
1.  0000 
8.  0000 
1  0.0000 
12. 0000 
16. 0000 
5.0000 
5.0000 
16. 0000 
1  1.  0000 
1  1.  0000 
3.0000 
8. 0000 
1  1.  0000 
1 1. 0000 
1. 0000 
16. 0000 
12.0000 
16.0000 
6. 0000 
13. 0000 
1. 0000 
10.  0000 
14. 0000 
3. 0000 
6.0000 
2. 0000 
5.0000 
10.0000 


10.3923 

10.3923 
5.  1962 
7.9373 

10. 3923 
4.2426 

10.3923 
5.  1962 

9.4868 
5.  1962 
12.0000 
9.0000 

10.3923 
12.0000 
6.0000 

9.9499 
11.  2250 
12.0000 
9.0000 
6.0000 

6.7082 
5. 1962 
3.0000 
8.4853 

9. 4868 

10.3923 

12.0000 

6.7082 

6.7082 

12.0000 

9.9499 

9.9499 
5.  1962 
8.4853 

9.9499 

9.9499 
3.0000 
12.0000 

10.3923 
12.0000 
7. 3485 
10.  8167 
3.0000 

9.4868 
11. 2250 
5.  1962 
7.3485 
4.2426 

6.7082 

9.4868 


28.3268 
28.6517 
8.3412 
18.0101 
29.6504 
6.0171 
30.3348 
8.8310 
26.  3044 
9.0338 
41. 2601 
24.7401 
32.4664 
42.6775 
12. 3539 
31.0273 
39.0649 
44.6212 
26.7486 
13.0572 
16.1173 
10.3242 
3. 9324 
25.3993 
31.4075 
37.4529 
49.2023 
17.3793 
17.5642 
50.7896 
36.4898 
36.8728 
1 1.6048 
28.  2280 
38.0353 
38.4273 
4.  5559 
55. 1623 
42.8717 
56. 2865 
23. 3980 
47.  5148 
4. 8388 
38.  1969 
52.3615 
13.  2271 
24. 8234 
9.4012 
21.  4917 
40.4908 


2.3606 
2.3876 
2.7804 
2. 5729 
2.4709 
3.0085 
2.5279 
2.9437 
2.6304 
3. 0113 
2.5788 
2. 7489 
2.7057 
2.6673 
3.0885 
2.8207 
2.7904 
2.7888 
2.9721 
3.  2643 
3.  2235 
3. 4414 
3.9324 
3. 1749 
3.  1407 
3. 1211 
3.0751 
3.4759 
3.5128 
3. 1744 
3.3173 
3.3521 
3.8683 
3.  5285 
3.4578 
3.  4934 
4.5559 
3.4476 
3.5726 
3.5179 
3.8997 
3.6550 
4. 8388 
3. 8197 
3.7401 
4.4090 
4.1372 
4.7006 
4.2983 
4. 0491 


76 

77 

70 

79 

00 

81 

62 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

113 

114 

115 

1 16 

117 

1 18 

119 

120 

121 

122 

123 

124 

125 


0.9216 
0.3289 
0. 8876 
0.2178 
0.4276 
0.5597 

1.  0990 
1. 6187 
3.  1780 
0. 2081 
0.7665 
0.6022 
0. 4930 
0.  3632 
0.2854 
0. 4006 
0. 2296 
2.7077 
0.  3801 
2.6151 
0.  1606 
0.  2297 
1.2417 
0. 3488 
0.2 182 
0.5900 
2.3206 
0.1426 
0.7481 
0.5518 
0.3101 

2.  1354 
0.3001 
0. 1879 
0.1452 
0. 4001 
0. 1789 
0.4842 
0.3177 
0.9379 
0.2637 
0. 1398 
1.7887 
0.  1956 
0.8667 
0.2438 
0.4201 
0.8272 
0.4073 
0.2292 


4. 0000 
11.0000 
4.0000 
16. 0000 
8.0000 
6. 0000 
3.0000 
2.0000 
1. 0000 
15.0000 
4. 0000 
5. 0000 
6. 0000 
8. 0000 
10.  0000 
7. 0000 
12. 0000 
1.0000 
7.0000 
1.0000 
16. 0000 
1  1.0000 
2. 0000 
7. 0000 
11.0000 
4.0000 
1.0000 
16.0000 
3.0000 
4.0000 
7.0000 
1.0000 
7.0000 
11.  0000 
14.0000 
5.0000 
11.0000 
4.0000 
6.0000 
2. 0000 
7.0000 
13.0000 
1.0000 
9.0000 
2.0000 
7.  0000 
4.0000 
2.  0000 
4. 0000 
7.0000 


6.0000 

9.9499 
6.0000 
12.0000 
8. 4853 
7.3485 
5. 1962 

4.2426 
3.0000 
11.6190 
6.0000 
6.7082 
7. 3485 
8.4853 
9.4868 

7.9373 
10. 3923 
3.0000 

7.9373 
3.0000 
12.0000 

9.9499 

4.2426 

7.9373 

9.9499 
6. 0000 
3.0000 
12.  0000 
5.  1962 
6. 0000 

7.9373 
3.0000 

7.9373 

9. 9499 
11. 2250 
6.7082 

9.9499 
6. 0000 
7.  3485 

4. 2426 

7.9373 
10. 8167 

3.  0000 

9.0000 

4. 2426 

7.9373 

6. 0000 

4.2426 

6. 0000 

7.9373 


17.9437 
44. 9928 
18.2878 
64.4933 
34. 7056 
27.0443 
14.6877 
10. 3336 
5.6372 
64.2834 
19.6980 
24.2694 
28.8395 
37.6953 
46.5287 
34. 0286 
55.8810 
6. 1138 
34.9453 
6. 2224 
75.2350 
53.9564 
11.  8168 
36.5007 
55.3723 
22. 4861 
6.6102 
79.9024 
17.  8439 
23.  2612 
38.7357 
6. 8944 
39.  3867 
59.7316 
74.9896 
29.8499 
61.  2218 
24. 8512 
36. 0207 

13.  6199 
42. 0454 
74.  1927 

7.  5410 
54. 0145 
14. 1748 
43.  7518 
26. 7037 

14.  5132 
27. 1245 
45.  1417 


4.4859 

4.0903 

4.5719 

4.0308 

4.3302 

4.5074 

4.8959 

5.  1668 
5.6372 
4. 2856 

4.  9245 
4.8539 
4.8066 
4. 7119 
4.6529 
4. 8612 
4.6567 

6.  1138 
4.9922 
6. 2224 
4.7022 
4.9051 
5.9084 

5.  2144 
5.0338 
5.6215 
6.6102 
4.9939 
5.9480 
5.  8153 
5.5337 
6. 8944 
5.6267 
5. 4301 
5.3564 
5.9700 
5.5656 
6.2128 
6.0035 
6. 8099 
6.0065 
5.7071 
7. 5410 
6.0016 
7. 0874 
6. 2503 
6.6759 
7. 2566 
6.7811 
6.4488 


126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 

161 

162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 


138 


0. 1436 

11.0000 

9. 

0. 1729 

9.0000 

9. 

0.7663 

2.0000 

4. 

1.5096 

1.0000 

3. 

0.1859 

8.0000 

8. 

0.1465 

10.0000 

9. 

0. 1203 

12.0000 

10. 

0. 1580 

9.0000 

9. 

0.2001 

7.0000 

7. 

0. 2760 

5.0000 

6. 

0.  1700 

8. 0000 

6. 

0. 1031 

13.0000 

10. 

0.0825 

16.0000 

12. 

0. 0868 

15.0000 

11. 

1.2828 

1.0000 

3. 

0. 1580 

8.0000 

8. 

1.2462 

1.0000 

3. 

0.0945 

13.0000 

10. 

1.2109 

1.0000 

3. 

0.3979 

3.0000 

5. 

0. 1681 

7.0000 

7. 

0.0829 

14.0000 

11. 

0.3814 

3.0000 

5. 

0.2256 

5.0000 

6. 

0.0856 

13.0000 

10. 

0.0686 

16.0000 

12. 

0. 1545 

7.0000 

7. 

0.0821 

13.0000 

10. 

0.  1315 

8.0000 

8. 

0. 0741 

14.0000 

11. 

0.0640 

16.0000 

12. 

0.0918 

11.0000 

9. 

0.0905 

11.0000 

9. 

0.0893 

11.0000 

9. 

0.2423 

4.0000 

6. 

0. 0869 

11.0000 

9. 

0.2359 

4.0000 

6. 

0. 1551 

6.0000 

7. 

0.0918 

10.0000 

9. 

0.  0697 

13.  0000 

10. 

0.2981 

3.0000 

5. 

0. 1261 

7.0000 

7. 

0. 0968 

9.0000 

9. 

0. 0661 

13.0000 

10. 

0.0530 

16.0000 

12. 

0. 0523 

16.  0000 

12. 

0. 1033 

8.0000 

8. 

0.0906 

9.0000 

9. 

0.0503 

16. 0000 

12. 

0.3975 

2.0000 

4. 

9499 

68.4225 

6. 2202 

0000 

57.5011 

6.3890 

2426 

15.0862 

7.5431 

0000 

8.2170 

8.2170 

4853 

52.9044 

6.6130 

4868 

65.  2118 

6.5212 

3923 

77.  5194 

6.4599 

0000 

60.  1889 

6.6877 

9373 

48. 3499 

6.9071 

7082 

36. 0159 

7. 2032 

4853 

55.  3493 

6.9187 

8167 

86.5519 

6.6578 

0000 

105.  3780 

6.  5861 

6190 

100. 0795 

6.6720 

0000 

6.9225 

8.9225 

4853 

57.4299 

7.  1787 

0000 

9.0540 

9.0540 

8167 

90.4492 

6.9576 

0000 

9. 1865 

9.  1865 

1962 

24.5588 

8. 1863 

9373 

52.  8028 

7.5433 

2250 

99.5967 

7.1141 

1962 

25.  0930 

8.3643 

7082 

39.8857 

7.9771 

8167 

95. 1075 

7.3160 

0000 

115.7561 

7.2348 

9373 

55.  1047 

7.8721 

8167 

97.  1407 

7. 4724 

4853 

63.  0241 

7.8780 

2250 

105.3869 

7.5276 

0000 

119. 8814 

7. 4926 

9499 

85. 8229 

7.8021 

9499 

66. 4180 

7. 8562 

9499 

87. 0152 

7.9105 

0000 

35. 2760 

8.8190 

9499 

88.2160 

8. 0196 

0000 

35.  7612 

8. 9403 

3485 

51.7659 

8.6276 

4668 

82. 5855 

8. 2586 

8167 

105. 4965 

8. 1 151 

1962 

28. 4234 

9.4745 

9373 

61.  0821 

8.7260 

0000 

77.  1265 

8.  5696 

8167 

108. 3613 

8.3355 

0000 

131.8308 

8. 2394 

0000 

132. 7069 

8. 2942 

4853 

71.  2145 

8.9018 

0000 

79.7237 

8.8582 

0000 

135. 3535 

8. 4596 

2426 

21.0302 

10.5151 

176 

177 

178 

179 

180 

181 

182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 

201 

202 

203 

204 

205 

206 

207 

208 

209 

210 

211 

212 

213 

214 

215 

216 

217 

218 

2  19 

220 

221 

222 

223 

224 

225 


139 


0.  0872 

9.0000 

0.1937 

4.0000 

0.0588 

13.0000 

0.2517 

3.0000 

0. 1864 

4.0000 

0.7362 

1.0000 

0.  1817 

4.0000 

0.0897 

8.0000 

0.0709 

10.0000 

0.  1 167 

6.0000 

0.  0532 

13.0000 

0.2276 

3.0000 

0.0674 

10.0000 

0.0476 

14.0000 

0. 1316 

5.0000 

0.0406 

16.0000 

0.0642 

10.0000 

0.3170 

2.0000 

0.3132 

2.0000 

0. 0413 

15.0000 

0.0470 

13.0000 

0.0755 

8.0000 

0. 1 193 

5.0000 

0.0536 

11.0000 

0.1165 

5.0000 

0.0576 

10.0000 

0.0948 

6.0000 

0.0511 

11.0000 

0.0505 

11.0000 

0.0343 

16.0000 

0.0904 

6. 0000 

0.0383 

14.0000 

0.0408 

13.0000 

0.0655 

8.0000 

0.0324 

16.0000 

0.0853 

6.0000 

0. 1686 

3.  0000 

0.0333 

15.0000 

0.0353 

14. 0000 

0. 0326 

15.  0000 

0.0805 

6.0000 

0. 0434 

11.0000 

0.0674 

7. 0000 

0.4667 

1.0000 

0.0288 

16.0000 

0.4562 

1.  0000 

0.0322 

14. 0000 

0.0496 

9.0000 

0.1103 

4. 0000 

0.0727 

6.0000 

9.0000 

81.3036 

6. 0000 

39.5102 

10.8167 

114.9538 

5. 1962 

30.9631 

6.0000 

40.  2834 

3.0000 

11. 8171 

6.0000 

40.8033 

8.  4853 

76. 4770 

9.4868 

94. 1518 

7.3485 

59.7951 

10.  8167 

120. 9852 

5.  1962 

32.5827 

9. 4868 

96.5606 

11. 2250 

131. 8974 

6.7082 

52.3992 

12.0000 

150.8694 

9. 4868 

99.  0014 

4.2426 

23.  5789 

4.2426 

23.7250 

11.6190 

145.7934 

10.8167 

128.7533 

8.4853 

83.4607 

6.7082 

55. 0499 

9.9499 

112. 6709 

6.  7082 

55.7234 

9. 4868 

104. 61 14 

7. 3485 

66.4175 

9.9499 

115.  4288 

9.9499 

116. 1238 

12.0000 

164. 3162 

7. 3485 

68.  0294 

11.2250 

147.2022 

10. 8167 

138.4140 

8.4853 

89.7047 

12.0000 

169. 2668 

7.3485 

70.0733 

5.  1962 

37. 9261 

11.6190 

162.  4096 

11.  2250 

153.3959 

11.6190 

164.  3150 

7.3485 

72. 1495 

9.9499 

125. 3620 

7.  9373 

83.  8334 

3.0000 

14. 8825 

12.  0000 

179.4040 

3. 0000 

15. 0543 

11. 2250 

160.6420 

9. 0000 

108.2049 

6. 0000 

52.  5402 

7. 3485 

75.9687 

9.0337 
9.8775 
6.8426 
10.3210 
10.0709 
11.8171 
10.2008 
9.5596 
9.4152 
9.9659 
9.3066 
10.8609 
9.6561 
9.4212 
10.4798 
9.4293 
9.9001 
11.7895 
11.8625 
9.7196 
9.9041 
10.4326 
11.0100 
10.2428 
11. 1447 
10.4611 
11.0696 
10.4935 
10.5567 
10.2698 
11.3382 
10.5144 
10.6472 
11. 2131 
10.5792 
1  1.  6789 
12.6420 
10.8273 
10.9568 
10.9543 
12.0249 
1  1.  3965 
1  1.  9762 
14.8825 
11.2127 
15. 0543 
1  1.4744 
12.0228 
13. 1351 
12.6615 


I 


326 

327 

328 

329 

330 

331 

332 

333 

334 

335 

336 

337 

338 

339 

340 

341 

342 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

373 

374 

375 


142 


0.0145 
0.0105 
0.0131 
0.0389 
0.0193 
0.0218 
0.0509 
0.  0167 
0.0744 
0.0164 
0.0365 
0.0091 
0.0160 
0.0178 
0.0101 
0.0087 
0.0116 
0.0275 
0.0682 
0.0270 
0.0191 
0.0083 
0.0329 
0.0109 
0.0216 
0.0428 
0.0085 
0.0126 
0. 0313 
0.0207 
0.0123 
0.  0203 
0.0086 
0. 0080 
0.0132 
0.0079 
0. 0146 
0.  0145 
0.0384 
0.0114 
0. 0087 
0.0140 
0.  0074 
0. 0276 
0.0073 
0.0181 
0.0098 
0.0  153 
0.0177 
0.  0131 


11.0000 
15.0000 
12.0000 
4.0000 
8.  0000 
7.0000 
3.0000 
9.0000 
2.0000 
9.0000 
4.0000 
16.0000 
9.0000 
8.0000 
14. 0000 
16.0000 
12.0000 
5.0000 
2.0000 
5.0000 
7.0000 
16.0000 
4.0000 
12.0000 
6.0000 
3.0000 
15.0000 
10.0000 
4.0000 
6.0000 
10.0000 
6.0000 
14.0000 
15.0000 
9.  0000 
15.0000 
8.0000 
8.0000 
3.0000 
10. 0000 
13.  0000 
8.0000 
15.0000 
4.0000 
15.0000 
6.0000 
11.0000 
7.0000 
6.0000 
8.0000 


9.9499 
11.6190 
10.3923 
6. 0000 
8.  4853 
7.9373 
5.  1962 
9.0000 
4. 2426 
9. 0000 
6. 0000 
12.  0000 
9.0000 
8. 4853 
11.  2250 
12.0000 
10. 3923 
6.7082 
4. 2426 
6.7082 
7.  9373 
12.0000 
6.0000 
10.3923 
7.3485 
5. 1962 
11.6190 
9.4868 
6.0000 
7.3485 
9. 4868 

7.  3485 
11.  2250 
11.6190 

9.0000 
11.6190 
8. 4853 
8.4853 
5. 1962 
9.4868 
10. 8167 
8.4853 
11. 6190 
6. 0000 
11.6190 
7.3485 
9.9499 
7.9373 
7. 3485 

8.  4853 


218.2586 
290.7210 
238.3285 
89.0698 
166.6244 
148. 4071 
69.8341 
187. 0061 
49. 1140 
189. 4807 
91.8888 
322.  4315 
192.0082 
173.  4085 
289. 2186 
328. 1514 
253.5781 
115.6913 
51.3213 
116.7033 
158. 5338 
336. 8383 
96.8375 
261.  4045 
140.4336 
75.  8083 
324.4671 
225.4156 
99. 3673 
143.4757 
228. 3220 
144.7030 
312.5748 
334.2946 
211. 1276 
337. 1332 
191.4476 
192.  2550 
80. 1962 
237. 1690 
302. 1882 
195. 5042 
347. 1765 
105. 8552 
350. 0770 
153. 4614 
266. 2042 
177.  7266 
155. 3766 
202. 0963 


19.8417 
19. 3814 
19.8607 

22.  2674 
20.8280 
21.2010 
23.2780 

20.  8673 
24.5570 

21.  0534 
22.9722 
20. 1520 
21.3342 
21.6761 
20.6585 
20.5095 
21.  1315 

23.  1383 
25.6606 
23. 3407 
22.6477 
21.0524 
24. 2094 
21.7837 
23.4056 
25. 2961 
21.6311 
22. 5416 

24.8418 
23.9126 
22.8322 
24. 1172 
22.3268 
22.2863 
23.4586 
22.4755 
23.9310 
24. 0319 
26.7321 
23. 7169 
23.  2452 
24.4380 
23. 1451 
26. 4638 
23.3365 
25.  5769 
24. 2004 
25.3895 
25.8961 
25. 2620 


376 

377 

378 

379 

380 

381 

382 

383 

384 

385 

386 

387 

388 

389 

390 

391 

392 

393 

394 

395 

396 

397 

398 

399 

400 

401 

402 

403 

4  04 

405 

406 

407 

408 

409 

410 

411 

412 

413 

414 

415 

416 

417 

418 

419 

420 

421 

422 

423 

424 

425 


0.0075 
0.0069 
0.0103 
0. 1018 
0.0084 
0.0072 
0.0083 
0.0197 
0.0082 
0.0061 
0.0963 
0.0478 
0.0316 
0.0085 
0.0117 
0.0116 
0.0131 
0.0304 
0.0056 
0.0100 
0.0074 
0.0088 
0.0097 
0.0217 
0.0108 
0.0285 
0.0065 
0.0065 
0.0167 
0.0064 
0.0206 
0.0068 
0.0135 
0.0804 
0.0160 
0.0079 
0.0112 
0.0156 
0.0387 
0.0064 
0.0085 
0.  0095 
0. 0375 
0.0106 
0. 0739 
0. 0046 
0. 0091 
0.0045 
0.0055 
0.0089 


14.0000 
15.0000 
10.  0000 
1.0000 
12.0000 
14.0000 
12.0000 
5.0000 
12.  0000 
16.  0000 
1.0000 
2.0000 
3.0000 
11.0000 
8.  0000 
8.0000 
7.0000 
3.0000 
16.0000 
9.0000 
12.0000 
10.0000 
9.0000 
4.0000 
8.0000 
3.0000 
13.0000 
13.0000 
5.  0000 
13.0000 
4.0000 
12.0000 
6.0000 
1. 0000 
5.0000 
10.0000 
7. 0000 
5.0000 
2.0000 
12.0000 
9.0000 
8.0000 
2.0000 
7.0000 
1.0000 
16.0000 
8.0000 
16.0000 
13.  0000 
8.0000 


11. 2250 
11.6190 
9.  4868 
3.0000 
10. 3923 
11. 2250 
10.3923 
6.7082 
10. 3923 
12.0000 
3.0000 
4.2426 
5. 1962 
9.9499 
8. 4853 
8.4853 
7.9373 
5.  1962 
12.0000 
9.0000 
10.3923 
9.4868 
9.0000 
6. 0000 
8.4853 
5. 1962 
10. 8167 
10. 8167 
6. 7082 
10. 8167 
6.0000 
10.3923 
7.3485 
3. 0000 
6.7082 
9.4868 
7.9373 
6.7082 
4.2426 
10. 3923 
9.0000 
8. 4853 
4. 2426 
7.9373 
3.0000 
12.  0000 
8. 4853 
12.  0000 
10. 8167 
8.4853 


336.9713 
360.3380 
250. 2889 
32.  1526 
297.7060 
343.  9351 
300. 1412 
136.8934 
302.  5878 
394.8866 
33.0801 
6  1.  4488 
88.5092 
285. 3096 
214.  7961 
215.6586 
191.9905 
90. 3006 
409.4144 
243.7023 
317.5083 
270.  1981 
246.6200 
119.5422 
223.5112 
93.  2100 
349.6675 
351.0436 
148. 8506 
353.8052 
122. 8745 
331.5509 
178.0443 
36.  2354 
152.  3736 
285.4337 
207.  6954 
154. 1530 
68. 3488 
341. 9855 
264. 5501 
238.7892 
69. 4074 
213. 3618 
37.8036 
454. 8307 
243. 3946 
458.3055 
380. 6686 
246. 1825 
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24.0694 

24.0225 

25.0289 

32. 1526 

24. 8088 

24. 5668 

25.0118 

27.3767 

25. 2156 

24.6804 

33.0801 

30.7244 

29.5031 

25.9372 

26. 8495 

26.9573 

27.4272 

30.  1002 

25.5884 

27. 0780 

26.4590 

27.0198 

27.4022 

29.  8855 

27.9389 

31.0700 

26.8975 

27.0034 

29.7701 

27.2158 

30.7186 

27.6292 

29.6741 

36. 2354 

30.4747 

28. 5434 

29.6708 

30.8306 

34.  1744 

28.  4988 

29.3945 

29.8486 

34.  7037 

30.4803 

37.8036 

28. 4269 

30.4243 

28.6441 

29. 2822 

30.7728 


*26 

*27 

*28 

*29 

*30 

*31 

*32 

*33 

*34 

435 

436 

*37 

*38 

*39 

440 

*41 

*42 

*43 

*44 

*45 

446 

447 

448 

449 

*50 

451 

*52 

*53 

454 

455 

456 

457 

*58 

459 

*60 

461 

462 

463 

464 

*65 

466 

*67 

468 

469 

470 

471 

472 

473 

474 

*75 
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0.00*4 
0.0702 
0.0696 
0.0173 
0.0049 
0.0085 
0.0085 
0.0075 
0.0061 
0.0132 
0.0073 
0.0130 
0. 0065 
0.  0054 
0.0046 
0.0211 
0.0314 
0.0125 
0.0062 
0.0077 
0.0153 
0.0038 
0. 0120 
0.0050 
0.0037 
0.0042 
0.0053 
0.0097 
0.0041 
0.0071 
0.  0063 
0.0113 
0.0080 
0.0278 
0.0050 
0.0046 
0.0060 
0.0090 
0.0077 
0. 0089 
0.0035 
0. 0263 
0.0065 
0. 0065 
0.0086 
0.0032 
0.0507 
0.0056 
0.  0250 
0.0033 


16.0000 

1.0000 

1.0000 

4.0000 

14.0000 

8.0000 

8.0000 

9.0000 

11.0000 

5.  0000 
9.0000 
5.0000 

10.0000 
12.0000 
14.0000 
3.0000 
2.0000 
5.0000 
10. 0000 
8.0000 
4.0000 
16.0000 
5.0000 
12.0000 
16. 0000 
14.  0000 
11.0000 
6.0000 
14.0000 
8.0000 
9.0000 
5.0000 
7.0000 
2.0000 
11.0000 
12.0000 
9.0000 
6.0000 
7.0000 

6.  0000 
15.0000 

2.0000 
8.0000 
8.0000 
6. 0000 
16. 0000 
1.0000 
9.0000 
2.0000 
15.0000 


12.0000 
3. 0000 
3.0000 
6. 0000 
11.  2250 
8.4853 

8.  *853 
9.0000 
9.9*99 

6.  7082 
9. 0000 
6.7082 
9. 4868 

10. 3923 
11.2*30 
5.  1962 
«. 2426 
6.7082 

9.  *868 
8. *853 
6.0000 

12.0000 
6.7082 
10.3923 
12.0000 
11.  2250 
9.9*99 

7.  3485 
11.2250 

8.4853 
9.0000 
6.7082 
7.9373 
*. 2*26 
9.9*99 
10. 3923 
9.0000 
7.  3485 
7.9373 

7.  3485 
11.6190 

4.2*26 

8.  4853 
8. 4853 
7. 3*85 

12. 0000 
3.0000 
9. 0000 
4.  2426 
1 1.  61*0 


*63. 5*66 
38.  8216 
38.968* 
134.1998 
*16.5789 
25 1 . 8 1 37 
252.7595 
282. 1588 
339.6170 
167.5711 
285.3350 
168. 8233 
316. 1968 
374.42*0 
*32.40*3 
108.5672 
75. 9585 
172.  6129 
323.2835 
265.2435 
142.9468 
501.212* 
175.  8086 
388.4*89 
506.73*1 
*50.2121 
362.8951 
210. 83*0 
*55. 1421 
275.0872 
307.  0439 
181.6476 
2*6. 5741 
80.  8092 
373.5322 
*05.6792 
313.7392 
218.5636 
251.93*8 
220.  1292 
506. 0466 
83. 1531 
288. 2001 
289.2236 
224.0720 
546. 3857 
45.7480 
326. 2352 
85.  2364 
S??. 4059 


28.9717 
36.8216 
38.9684 
33.5*99 
29.7556 
31.4767 
31.59*9 
31.3510 
30.8743 
33. 5142 
31.7039 
33.7647 
31.6197 
31.  2020 
30.8860 
36. 1891 
37.9793 
34.5226 
32.328* 
33. 1554 
35.7367 
31.3258 
35. 1617 
32.3707 
31.6709 
32.1580 
32.9905 
35.  1390 
32.5101 
3*. 3859 
34. 1160 
36.3295 
35.2249 
*0.4046 
33.9575 
33.8066 
34. 8599 
36.4273 
35.9907 
36.6882 
33.7364 
*1.5765 
36. 0250 
36.1530 
37. 3453 
34.1491 
45.7480 
36.248* 
*2.  6182 
34.8271 


476 

477 

478 

479 

480 

481 

482 

483 

484 

485 

486 

487 

488 

489 

490 

491 

492 

493 

494 

495 

496 

497 

498 

499 

500 

501 

502 

503 

504 

505 

506 

507 

508 

509 

510 

511 

512 
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0.0070 
0.0035 
0.0243 
0.0035 
0.0060 
0.0034 
0.0032 
0.0067 
0.0029 
0.0232 
0.0461 
0.0030 
0.0076 
0.0075 
0.0034 
0.0032 
0.  0147 
0.0146 
0.0029 
0.0036 
0.0054 
0.0071 
0. 0047 
0.0026 
0.  0419 
0.0035 
0.0034 
0.0103 
0.0031 
0. 0029 
0.0067 
0.0400 
0.0031 
0.  0028 
0.0098 
0.0028 
0. 0129 


7.0000 

7.9373 

14.0000 

11.  2250 

2.0000 

4.2426 

14.0000 

11.2250 

8.0000 

8.4853 

14.0000 

11.  2250 

15.0000 

11.6190 

7.0000 

7.9373 

16. 0000 

12.0000 

2.0000 

4.2426 

1.0000 

3.0000 

15.0000 

11.6190 

6.0000 

7.3485 

6.0000 

7.  3485 

13.0000 

10.  8167 

14.0000 

11.  2250 

3.0000 

5. 1962 

3.0000 

5. 1962 

15.0000 

11.6190 

12.0000 

10. 3923 

8.0000 

8.4853 

6.0000 

7. 3485 

9.0000 

9.0000 

16. 0000 

12.0000 

1.0000 

3.0000 

12.0000 

10.3923 

12.0000 

10.3923 

4. 0000 

6. 0000 

13.0000 

10. 8167 

14.0000 

11. 2250 

6.0000 

7.3485 

1.0000 

3.0000 

13.0000 

10. 8167 

14.  0000 

11.  2250 

4.0000 

6.0000 

14.0000 

11.  2250 

3.0000 

5. 1962 

262. 8604 

37.5515 

493.9922 

35.  2852 

86.4406 

43.2203 

497.4593 

35.5328 

300.6247 

37.5781 

500.9408 

35.7815 

535. 3425 

35.6895 

269. 3603 

38.4800 

571.8179 

35.7386 

88.5719 

44. 2860 

48. 0402 

48. 0402 

544.6977 

36.3132 

238.6094 

39.7682 

239.4329 

39.9055 

483.0563 

37. 1582 

518.5639 

37.0403 

130.  1654 

43.3885 

130.6121 

43.5374 

557.9564 

37. 1971 

456. 9109 

38.0759 

317.6751 

39.7094 

246.0816 

41.0136 

355.7140 

39.5238 

602.0170 

37.6261 

50.  3990 

50.  3990 

466. 3276 

38.8606 

467.9082 

38.9923 

174.6764 

43.6691 

506.6400 

38.9723 

543.8452 

38.8461 

253.  6908 

42.  2816 

51.6036 

51.6036 

513.5010 

39. 5001 

551.2000 

39.3714 

178.8342 

44.7085 

554.8994 

39.6357 

139.2730 

46. 4243 

V.  CONCLUSION:  IMPLEMENTATION  OF  THE  GKD  ALGORITHM 


WITH  SAMPLING  FOR  A  LARGE-SCALE  INVENTORY  SYSTEM 

1.  INTRODUCTION 

This  chapter  concludes  the  studies  in  Chapter  II.  Ill  and  IV, 
in  which  we  examined  the  problem  of  specifying  single-item  service 
objectives  in  a  multi-item  inventory  system  subject  to  an  overall 
service  objective.  A  list  of  the  notation  we  use  may  be  found  at 
the  end  of  Chapter  II. 

In  that  chapter  we  formulated  this  inventory  problem  as  a  con¬ 
strained  nonlinear  program  (NLP)  and  developed  the  Generalized 
Knapsack  Duality  (GKD)  Algorithm  to  solve  it.  This  alqorithm  finds 
the  Lagrange  constraint  multiplier  of  the  NLP  and  from  it  computes 
the  (s,S)  policies  for  all  items  in  the  inventory  system.  In  Chapter 
III  we  reported  computational  experience  with  several  32-item  inven¬ 
tory  systems  which  have  a  structure  typical  of  many  real-world 
inventory  systems.  This  structure  is  called  an  80/20  value  structure, 
and  is  one  in  which  twenty  percent  of  the  items  represent  eighty 
percent  of  the  value  of  the  systems,  where  value  is  taken  to  be  hy. 
These  items  are  the  high-value  items,  and  the  other  eighty  percent 
are  the  low-value  items.  In  order  to  investigate  the  performance  of 
the  GKD  Algorithm,  we  assumed  that  the  demand  distribution  is  Negative 
Binomial  and  compared  the  (exact)  expected  holding  cost  of  the  (s,S) 
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policies  produced  by  the  GKD  Algorithm  with  the  (s,S)  policies  pro¬ 
duced  by  a  method  popular  with  managers,  the  Identical  Service 
Approach  (The  replenishment  costs  for  these  methods  are  the  same, 
and  represent  about  half  of  the  total  expected  cost  ).  Given  a  tar¬ 
geted  overall  service-level  a  (we  targeted  85%),  the  Identical 
Service  Approach  sets  (s,S)  policies  by  the  Power  Approximation  of 
Ehrhardt  (Mosier  [1981])  so  that  every  item  has  a  service-level  a. 

(The  GKD  Algorithm  varies  individual  service-levels  while  still 
maintaining  an  average  overall  service-level  a.)  We  showed  that 
there  is  a  significant  cost  savings  when  using  the  GKD  Algorithm 
rather  than  the  Identical  Service  Approach  to  manage  these  inventory 
systems. 

In  Chapter  IV  we  reported  several  sensitivity  tests  performed 
on  the  base-case  and  worst-case  inventory  systems  of  Chapter  III. 

In  particular,  we  performed  sensitivity  tests  on  the  D^  and  on  the  lower 
bounds  for  S^.  We  also  examined  several  sampling  schemes  and  sizes 
for  inventory  aggregation  for  larger  (128-  and  512-item)  systems. 

Since  most  real-world  inventory  systems  contain  thousands  of  items, 
typically  managers  sample  from  the  system  and  make  decisions  about 
policies  for  all  items  in  the  system  based  on  the  sample.  We  ex¬ 
amined  sampling  schemes  (both  random  and  nonrandom)  of  several  sizes 
which  maintained  the  80/20  value  structure  and  the  overall  service 
constraint  weights.  We  summarize  the  conclusions  and  recommendations 
given  in  that  chapter  as  follows: 
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1.  Specify  all  D.)  using  Ehrhardt's  Power  Method  (Mosier  [1981]). 

2.  Specify  the  lower  bounds  on  Si  as  D.+rp>,  where  r  is  taken  as 
large  as  possible  while  still  including  as  feasible  the  targeted 
service-level  and  while  the  expected  holding  cost  is  significantly 
smaller  than  that  of  the  Identical  Service  Approach. 

3.  Sample  about  32  items  using  a  central  nonrandom  uniform  sampling 

scheme  to  create  the  inventory  subsystem  for  aggregation  (sample 

items  from  a  stratum  containing  items  1,2,...,N  by  taking  items 

N/2M,  3N/2M,  5N/2M, . . . , (2M-1 )N/2M,  where  the  fractions  are  rounded 

to  the  nearest  integer).  Recompute  the  holding  and  replenishment 

costs  to  maintain  the  value  structure  of  the  system,  and  recompute 

the  constraint  weights  to  maintain  the  constraint-weight  structure 

of  the  system  (multiply  the  costs  by  I  h-p./.Z  h.y.  and  the 

i  in  1  1  ip  1 
system  sample 

constraint  weights  by  N/M  for  each  of  the  M  items  sampled  from  1,2.... h)* 
In  Section  2  of  this  chapter  we  prescribe  a  detailed  implemen¬ 
tation  procedure  of  the  6KD  Algorithm  with  sampling  for  a  large-scale 
inventory  system  with  an  overall  service  objective.  In  Section  3  we 
report  the  results  of  implementing  this  procedure  for  the  512-item 
system  described  in  Chapter  IV.  Section  4  contains  conclusions  and 
directions  for  future  research,  and  the  Appendix  contains  a  listing 
of  the  FORTRAN  code  of  the  GKD  Algorithm  and  related  programs. 

2.  IMPLEMENTATION  PROCEDURE 

Based  on  the  conclusions  and  recommendations  of  Chapters  II, 

III.  and  IV,  we  offer  the  following  procedure  to  specify  (s,S) 
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policies  in  large-scale  inventory  systems.  Our  empirical  evidence 
suggests  that  this  procedure  is  a  significant  improvement  over  the 
Identical  Service  Approach  often  used  by  managers,  capable  of  sub¬ 
stantially  reducing  operating  costs  in  an  inventory  system  with  an 
overall  service  objective. 

The  procedure  is  as  follows: 

1.  Determine  the  desired  range  of  service-levels  and  the  targeted 
service-level.  This  is  a  managerial  decision,  but  typically  the 
range  is  80%  to  90%  with  a  target  of  85%. 

2.  Determine  r  as  follows.  For  r  =  0, .5, .6, .7, .8, .9  and  1,  input 
the  following  values  to  the  program  LSMISS  (see  the  Appendix 
to  this  chapter): 

TITLE  =  user-supplied  title 

LBAL  =  smallest  service-level  to  consider 

UBAL  =  largest  service-level  to  consider 

INCRAL  =  service-level  increment  between  LBAL  and  UBAL 

NSIG  =  8 

MAXFN  =  100 

NI  =  number  of  items  in  inventory  system 

SCHEME  =  1 

ENTIRE  =  0 

STRA(l)  =  NI 

NUMB(l)  =  about  32. 

Using  the  output  from  LSMISS,  graph  the  expected  holding  costs 
for  the  various  service-levels  for  each  value  of  r  and  for  the 
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Identical  Service  Approach  (compare  Figure  34).  Take  as  r  the 
largest  value  of  r  such  that  LBAL  is  achievable  and  such  that 
there  is  an  acceptable  cost  savings  below  that  of  the  Identical 
Service  Approach. 

NOTE  1: 

The  GKD  Algorithm  assumes  that  the  demand  distribution  is 
Logistic,  an  approximation  to  the  Normal  distribution  (Section 
2  of  Chapter  II).  If  the  actual  distribution  is  not  well- 
approximated  by  a  Normal  distribution,  then  the  service-level 
of  the  (s,S)  policies  produced  by  the  GKD  Algorithm  may  be  dif¬ 
ferent  from  the  specified  service-level.  The  user  may  there¬ 
fore  have  to  try  various  LBAL,  UBAL  and  INCRAL.  If  the  service- 
level  specifications  are  as  described  in  Step  (1)  and  the  demand 
distribution  is  well -approximated  by  a  Gamma  or  Negative  Binomial 
distribution,  our  experience  suggests  using 
LBAL  =  .75 
UBAL  =  .85 
INCRAL  =  ,01. 

NOTE  2: 

We  recommend  a  sample  size  ( NUMB ( 1 ) )  of  about  32  because  in  our 
experience  larger  samples  do  not  significantly  improve  the 
performance  of  the  sampling  scheme. 

NOTE  3: 

The  parameters  NSIG  and  MAXFN  are  used  to  terminate  IMSL  sub¬ 
routines  in  the  GKD  Algorithm.  We  have  found  the  recommended 
values  to  be  satisfactory  in  every  experiment  performed  in  this  paper. 
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3.  Using  the  output  from  LSMISS  and  linear  interpolation,  determine  A* 
for  a  targeted  service-level.  Input  the  following  values  to 
the  program  MISSNLSI  (see  the  Appendix  to  this  chapter): 

TITLE  =  user-supplied  title 
LS  =  A* 

NSIG  =  8 
MAXFN  =  100 

NI  =  number  of  items  in  system 

0  write  output  only  to  printer 

0UT  =  , 

1  write  output  to  both  printer  and  an  external  file 

The  output  of  MISSNLSI  contains  all  single-item  policy  specifi¬ 
cations. 

NOTE: 

The  parameter  OUT  allows  the  user  to  save  the  output  of  MISSNLSI 
in  an  external  file.  See  the  Appendix  to  this  chapter  for  details. 

3.  EXAMPLE 

In  the  previous  chapter  we  performed  Steps  (1)  and  (2)  of 
this  procedure  for  the  512-item  base-case  inventory  system,  and 
decided  to  use  r=.6.  We  now  perform  Step  (3)  for  this  system.  We 
used  the  output  from  LSMISS  for  the  scheme  described  in  Step  (2)  of 
the  procedure,  which  corresponds  to  Experiment  1  of  Table  3.  For 
this  experiment.  Figure  43  reports  the  values  of  A*  for  the  service- 
levels  near  85%.  our  targeted  service-level.  Linear  interpolation 
gives  A*  =  4207.  The  output  from  the  program  MISSNLSI  is  summarized 
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in  Figures  44,  45,  and  46.  The  operating  characteristics  of  the 
individual  items  follow  the  same  pattern  noted  in  Chapter  III.  in 
that  the  high-value  items  are  stocked  at  low  service-levels  and  the 
low-value  items  are  stocked  at  high  service-levels. 

4.  CONCLUSIONS,  FUTURE  RESEARCH  DIRECTIONS 
Using  the  GKD  Algorithm  of  Chapter  II  and  the  extensive  com¬ 
putational  and  sensitivity  experience  reported  in  Chapters  III  and 
IV,  we  have  described  a  detailed  implementation  procedure  of  the  GKD 
Algorithm  using  sampling  to  manage  a  large-scale  inventory  system  with 
an  overall  service  objective.  We  have  shown  that  for  inventory 
systems  with  a  structure  typical  of  that  found  in  many  real-world 
inventory  systems,  this  procedure  is  practical  and  results  in  a 
significant  cost  savings  below  that  of  the  Identical  Service  Approach 
often  used  by  managers. 

We  recommend  the  following  future  research  directions: 

1.  This  procedure  should  be  tested  in  a  statistical  environment  in 
which  the  moments  of  the  demand  distribution  are  estimated  from 
recent  demand  history.  In  light  of  the  results  reported  in 
Ehrhardt  [1976],  Estey  and  Kaufman  [1975],  Klincewicz  [1976a], 
MacCormick  [1974,  1977],  and  MacCormick  et  aj_.  [1977],  the 
performance  of  the  procedure  will  probably  degrade  somewhat, 
and  this  should  be  investigated.  Considering  the  cost  savings, 
however,  we  conjecture  that  in  a  statistical  environment  our 
procedure  will  still  result  in  significant  cost  savings  below 
that  of  the  Identical  Service  Approach. 
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2.  The  variable  r  should  be  included  in  the  6KD  optimization 
routine.  This  would  improve  accuracy  and  make  an  implementation 
procedure  simpler  in  that  Step  (2)  of  the  procedure  would  be 

el iminated. 

3.  The  GKD  Algorithm  assumes  the  demand  distribution  is  Logistic, 
while  demand  distributions  are  typically  skewed  like  a  Gamma  or 
Negative  Binomial  distribution.  As  reported  in  Section  3  of  Chapter 
III,  this  can  lead  to  significant  degradation  in  algorithm  per¬ 
formance  (although  we  showed  it  is  possible  to  improve  such 
performance).  Moreover,  as  discussed  in  Note  1  under  Step  (2) 

of  our  implementation  procedure,  the  actual  service-level  of 
the  policies  produced  by  the  GKD  Algorithm  may  differ  from  that 
specified.  We  recommend  that  the  Logistic  distribution  be  replaced 
by  a  distribution  which  better  approximates  a  Gamma  (or  Negative 
Binomial)  distribution.  The  distribution  to  use  is  not  obvious; 
an  examination  of  the  proof  of  Theorem  8  of  Chapter  II  suggests 
that  the  distribution  function  must  have  a  fairly  simple  form  in 
order  to  prove  algorithm  convergence  in  the  way  we  did  (A(0)  needs 
to  be  rational  in  6).  We  have  been  unable  to  find  a  distribution 
that  better  approximates  a  Gamma  distribution  and  yet  is  simple 
enough  to  allow  an  extension  of  the  proof  of  Theorem  8. 

4.  Inventory  systems  with  a  value  structure  different  than  80/20 

should  be  studied.  Based  on  previous  but  unreported  experiments, 
we  conjecture  that  the  algorithm  performance  will  degrade  for 
X/20  systems  as  X  decreases. 
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Figure  44.  51 2- Item  Inventory  System 
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APPENDIX  TO  CHAPTER  V 

This  appendix  contains  a  listing  of  the  FORTRAN  computer  program 
for  the  programs  LSMISS  and  MISSNLSI.  The  subroutines  EXACT  and  POWAPP 
are  the  same  for  both  of  these  programs,  and  the  subroutine  COMPAR 
is  slightly  different  for  the  two  programs. 


o  n 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PPOGBAH  LSBISS(NL)  :  LABGE  -  SCALE 

HULTI  -  ITEH  (S,S) 

NOBHAL  D EH  AND  DISTRIBUTION  APP20X IB ATION 
LOWER  BOUNDS  ON  BIG  S 

VBBSION  DATED  C7/20/83 

THE  INPUT  IS  PRC  B  CABDS 

TITLE  =  USEB  -  SUPPLIED  TITLE  (LE  36  ALPHA  CHARACTERS) 

NI  =  NUBBEB  Or  ITEBS  (NI  LE  512) 

K  =  SET-UP  ORDERING  COST 

L  =  LEADTIHE  i  1 

HO  =  UNIT  HOLDING  COST 

ED  =  BEAN  DEBAND 

SD  =  DEBAND  STANDARD  DEVIATION 

NT  =  CONSTRAINT  WEIGHT 

NSIG  =  NUHBER  OF  SIGNIFICANT  DIGITS  DESIRED  FROB  THE  IMSL 
BOUTIN  ES  ZBRENT  AND  ZFALSE 
EPS  =  1/(10  **  NSIG) 

EMFN  =  HAXIBUH  NUBBEB  OF  ITERATIONS  ALLOWABLE  BT  USER  FOR 
THE  IBSL  ROUTINES  ZEBENT  AND  ZFALSE 
I8R  =  (OUTPUT)  EREOB  PAR  AH  ETEB  FBOB  IBSL  ROUTINES  AND  EXACT 

THE  FOLLOWING  ARE  INPUT  TO  OB  COBPUTED  BY  GKD 

LITS  *  LITTLE  S 
BIGS  =  BIG  S 

D  =  REORDER  QUANTITY  BIG  S  -  LITTLE  S 
BSL  =  R  ♦  L  *  HU 

=  LOWER  BOUND  FOR  BIG  S  (HAY  BE  INCREASED  BY  GKD) 

AL  *  OVERALL  SE  R  V I  C  E- LEV  EL  (0  <  AL  <  1) 

THIS  PROGRAH  COHPUTES  POLICES  POR 

AL  =  LB  AL,  LBAL-UNCRAL,  LB  AL+2*INCH  AL,  .  .  .  ,OBAL 

THE  FOLLOWING  ARE  CCHFUTED  BY  EHRHARDT * S  POWER  HETHOD 

FLITS  =  LITTLE  S 
P BIGS  =  BIG  S 

THE  FOLLOWING  APF  CCHFUTED  EXACTLY 

HGKD  =  EXPECTED  HOLDING  COST  FOR  GKD  POLICIES 

1GKD  =  EXPECTED  TOTAL  COST  FOR  GKD  POLICIES 

AVSL  =  WEIGHTED-AVFFAGE  SERVICE-LEVEL  FOB  GKD  POLICIES 

HPOi  =  EXPECTED  HOLDING  CCST  FOR  POWER  POLICIES 
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C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


TPOf  *  EXPECTED  TOTAL  COST  FOR  POKER  POLICIES 

PAVSL  =  WEIGHTED” AVERAGE  SERVICE- LEVEL  POR  POKER  POLICIES 

TUB  PCLLOKI MG  PARAMETERS  ARE  OSED  IN  THE  STRATIFIED  SAMPLING 
SUBROUTINES  STSABP  AMD  BNDSTS 


SCHEME  - 


SAMP 
BUMS! 
ST  A (I) 
NUMB  (I) 


0  DC  ACT  DSE  AMI  SAMPLING  SCHEME 

1  DSE  5TBAT1FED  SAMPLING  SCHEME 

2  DSE  BANDOM  STRATIFED  SAMPLING  SCHEME 

NUHBEB  OF  ITEMS  IN  THE  SAMPLE 
NUMBER  OF  STBATA  TO  SAMPLE  FROM 
LAST  ITEM  IN  STBATA  I 

BUHBEB  CF  ITEMS  FBON  STATA  I  IN  TBE  SAMPLE 


I 

(  ON 
|  INPUT 


THE  FOLLOWING  V 1B2  ABLE  IS  NEEDED  ON  INPUT  IF  SCHEME  =  2 


S  F  ED  =  SEED  FOB  PS EU CO- BAN DC M  NUMBER  GENERATOR  GGUD  (IHSL) 

BEQUIEE  SEED  IN  (1,2147483674) 

ITEM  ( 1)  ,  .  .  ITEM  (SAMP)  ARE  THE  ITEMS  IN  THE  SAMPLE  ]  OUTPUT 
ENTIRE  =  0  OTHERWISE 

1  EVALUATE  EXPECTED  COSTS, ETC. ,  FOB  ENTIRE  INVENTORY 
SYSTEM  (THIS  IS  FOB  TESTING) 

tt**t******************* ****************************************** 


c 

INPUT  THE  SYSTEM  FAB 

AHFTERS  AS  FOLLOWS: 

c 

c 

TITLE 

c 

LB  AL 

UB  AL  INCRAL 

NS IG  HAIPN 

NX 

c 

NT  ( 1 ) 

D(l)  K(1) 

1(1)  HO  ( 1) 

HU  (  1)  SD  ( 1 ) 

R  (1) 

c 

NT  (2) 

D  (2)  K(2) 

L  (2)  HO  (2) 

HU  (2)  SD  (2) 

R  (2) 

c 

0 

•  0 

0  0 

0  0 

0 

c 

• 

0  0 

0  • 

0  0 

• 

c 

• 

0  0 

0  0 

0  0 

0 

c 

WT  (N) 

D(N)  K(N) 

L  (N)  HO  (N) 

MU  (N)  SD  (N) 

R(N) 

c 

SCHEME 

ENTIRE 

c 

SAMP 

NUHST 

_ 

c 

STB  A  (  1) 

STB  A  (2)  .  . 

.  STRA (NUHST)  | 

16  TO  A 

c 

NUMB  ( 1) 

NUMB  (2)  .  . 

.  NUMB ( NUHST)  _J 

LINE 

c 

SEED 

)  ONLY  IF  SCHEME  =  2 

L 

c 

THE  INPUT  FOBMAT  IS: 

c 

C 

XIXXYXXIXIXXXIXXXXXXXXXXXXXXXXXX 

c 

X. IXXXX 

X. XXXII  X.XXXXX  XXXXXX  XXXXXX  XXXXXX 

c 

XXI. XXIX  XXX. XXXX 

0  0  0 

XXX.XXXX  XXX.XXXX 

c 

XXX.XXXX  XXX. XXXX 

0  •  • 

XXX.XXXX  XXX.XXXX 
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• 

m 

•  • 

• 

•  • 

m 

• 

•  • 

• 

•  • 

• 

• 

•  • 

• 

•  • 

XXX.XXXX  XXX.XXXX 

•  * 

• 

XXX.XXXX  XXX.XXXX 

XX  XX 

XXXX 

XXXX 

1 

XXXX 

XXXX 

.  XXXX 

1 

NOT  NEEDED  IF 

XXXX 

XXXX 

.  XXXX 

1 

SCHEME  *  0 

C  1XXXXXXXXX  _l 

c 

C  ALL  DATA  STAB7S  IN  COLUMN  1  (EXCEPT  TITLE  STARTS  IN  COLUMN  2), 

C  AND  THBRB  IS  A  SINGLE  SPACE  BETWEEN  EACH  DATA  ENTBT 
C 

C  . . . . . 

c 

R  BA  L*  9  UTS  (512)  , FIGS (5 1 2)  ,BSL(512)  ,D(512> 

P  SAL*9  bSLl  (512)  t  D 1  (512)  ,BSLB  (512) 

B  SAL  *8  R  (512)  ,1(512)  ,HO(5l2)  .MU  (512)  ,  SD  (5  12)  ,HT(S12> 

REALMS  K1  (512)  ,L1  (512)  ,H01  (512)  #MU1  (512)  ,SDl  (512)  ,WT1  (512) 
B3IL*9  AL,1PS,LB A  L  ,  U  fcA  L,  I  NCR  A  L,  LS  (20)  , LSI  ,5 LI # SUM  , S E E D, B 
C 

I  ST  EGER  SI,  V,  H,N5IC, HAIFA, TITLE  (0)  ,  J,  HUM,  ITEM  (512)  ,  ENTIRE 
I NT  EG  ER  SCHEME  ,  SAMP  ,  BUB ST, STB  A  15 1 2) , «U MB ( 5 1 2 ) ,1SEE0 
C 

BEAD  (1,3)  (TITLE  (J)  ,J=1,8) 

3  FOBfl  AT  (8  A4) 

B BAD  (1, 5)  LBAL,UBAL,INCRAL,NSIG,HAXFN,NI 
5  FORMAT  (3(F7,5,1X),3  (16  ,  IX) ) 

EPS  =  10.  **  (-HSIG) 

C 

DO  10  I  =  1,NI 

READ  (1,8)  HT  (I)  ,D  (I)  ,  K  (I)  ,  L  (I)  ,  HO  (I)  ,MU  (I)  ,  S  D  (I)  ,B 
BSL(I)  =  D  (I )  -)  E* L  (I)  *MU(I) 

8  FORMAT(8  (F8. «,1X) ) 

10  CONTINUE 
C 

READ  (1,11)  SCHEME, ENTIRE 

11  FORMAT (12,11,12) 

IP  (SCHEME.  EQ.  0)  GO  TO  19 
F  EAD ( 1, 1 5)  SAMP, NUMST 
C 

C  NUMST,  STRA  AND  NUMB  MOST  EE  CHANGED  AS  FOLLOWS  FOR  INPUT  TO 

C  THE  SUBROUTINES  STSAMP  AND  RNDSTS SEE  THE  DOCUMENTATION  FOR 

C  THESE  SUBROUTINES 
C 

NUMST  =  NUMST-)  1 
srpA(i)  =  o 

NUMB  { 1 )  =  0 

R  EA  D  ( 1,  1  5)  (STRA  (I)  ,1=2,  NUMST) 

R  E  A  D  (  1 ,  1  5)  (NUMD  (I)  ,1=2,  NUMST) 


non 
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15  FORMAT  (16  (14,  II)  ) 

I  SEED  «=  0 

IP  (SCHEME.  EQ.  2)  BEAD(1,16)  ISEED 

16  FORM  AT (110) 

SBED  =  ISEED 

HAKE  THE  WEIGHTS  HI  INTO  EQUIVALENT  CONVEX  WEIGHTS  WT 

19  SUM  =  0. 

DO  20  I  '  1  » NI 

20  SUB  SUH-)WT  (I) 

SUB  =  1.0D0/SUB 
DO  25  I  =  1 #  N I 

25  WT  (I)  =  SUH*WT  (I) 

B  =  2+(UBAL-LEAL)/INCRAL 
IP  (SCHEBE.  EQ.  0)  GO  TO  300 

C  CR  BAT  E  THE  INVENTORY  SUBSYSTEM  VIA  STB ATIPIED  SAMPLING  AND  PUT 
C  TUB  SUBSTSTEB  IN  N,  Kl,  Li,  ...  ,BSL1 

IP  (SCHEME. EC* 1)  CALL  STS  AMP (NI , K, L, RO, MU , S D , D , WT , BSL , N , K 1 , L 1 , 

t  1  HO  1 ,HU 1 , SD1 , D1 , WT 1 , BSL 1 , SAMP , 

c  HUMST,STBA, NUMB, ITEM) 

IP  (SCHEME. EQ. 2)  C ALL  BNDSTS  (NI , K, L. HO, HU , SD , D,WT , BSL , W , K 1 , L 1 , 
r  HOI, HU1,SD1,D1,WT1,BSL1, SAMP, 

G  NUHST,STEA,  NUMB, ITEM, SEED) 

C 

DO  200  II  =  1 #  H 

AL  =  LBAL+  (II-1.  DO)  *INCBAL 

SLI  =  AL 

NOB  =  N 

DO  30  III  =  1 # ■ 

30  BS LB  (III)  =  ES1 1  (III) 

C  COMPUTE  LAHECA-STAB  (LSI)  USING  THE  GKD  ALGCBITHM  WITH  THE 
C  INVENTORY  SUBSYSTEM 

CALL  GKD (FUH, L 1, HC1 ,HU 1,SE1 , WT1, D1 ,SLI,LITS, BIGS, BSLB,EPS, 

G  NS1G  ,BAXPN,LSI) 

LS  (II)  =  LSI 

WRITE  (3,40)  SCHEBE,  (STRA(I)  ,1  =  2, NOMST) 

40  FORM  AT ( *  FESUL1S  FOB  INVENTORY  SU BS YSTEM ' //'  SCHEBE  ',12// 

&  •  STRATA  *,15(815/)/) 

WRITE{3,41)  (NUME  (I)  ,1=2,NUMST) 

41  FORMAT  ( '  NUMBER  *,15(815/)/) 

WRITE  (3, 42)  (ITEM  (I)  ,1=1, SAMP) 


noon  nnnno 
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4 2  FORMATf  ITER  *,15(815/)/) 

8BITE  (3,4  3)  1SEED 
4  3  FORM  AT ( '  SEED  =  *,120/) 

CCHPUTB  THE  (UNIFORM  5 EBV ICE- LET  EL)  POWER  APPBOXISATION  POLICIES 
IND  COMPARE  THEIR  ASSOCIATED  COSTS  WITH  THE  COSTS  ASSOCIATED  WITH 
ThB  GKD  POLICIES 

CALL  COUP AB ( H , K 1 , L 1 , HO  1 , HD  1 , SD 1 , WT 1 , BIGS , LIT S , TITLE , LSI , AL) 

IP  (ENTIRE. EQ. 0)  GO  TO  200 

COSPOT B  THE  OPTIMAL  POLICIES  ASSOCIATED  WITH  THE  MULTIPLIER 
LAHBDA-STAB  (LSI) 

CALL  BJGSLS  (N I , L , HO , BU , SD , D, BT , BSL, LSI , E PS , NSIG, H AX FN , B IG S, LITS) 
C 

WHITE  (3,05) 

45  PORSAT  (*  1  ') 

WBITE(3,50)  SCHEME,  (STBA(I)  ,I=2,WUHST) 

50  FOBS AT  ( *  RESULTS  FUR  INVENTORY  SYSTEM'//*  SCHEME  *,I2// 

e  *  STRATA  ',15(815/)/) 

WRITE  (3,51)  (NUMB  (1)  ,1  =  2,  NUMST) 

51  FORMAT  ('  NUMBER  ',15(815/)/) 

WRITE  (3,52)  (ITEM  (1)  ,1=1, SAMP) 

52  FORMAT  ( '  ITEM  ',15(815/)/) 

WRITE  (3,  53)  ISEED 

53  FORMAT  ( '  SEED  =  ',120/) 

C 

C  COMPUTE  THE  (UNIFORM  SERVICE-LEVEL)  POWER  APPROXIMATION  POLICIES 
C  AND  COMPARE  THEIR  ASSOCIATED  COSTS  WITH  THE  COSTS  ASSOCIATED  WITH 
C  THE  BIGSLS  POLICIES 
C 

CALL  COM PAR (NI , K , L, HO, HU , SD, WT, BIGS , LITS, TITLE , LSI , AL) 

200  CONTINUE 
STOP 
C 

300  DO  350  II  =  1 , H 

AL  =  LBAL-i  (II-  1.  DO)  *INCRAL 
SLI  =  AL 
NUH  =  NI 

DO  320  III  =  1 , N I 
320  BSLB  (III)  =  BSL  (III) 

C 

C  CCHPUTB  LAHEEA-STAR  (LSI)  USING  THE  GKD  ALGORITHM  WITH  THE 
C  ENTIRE  INVEN70RI  SYSTEM 
C 

CALL  GKD(NUH,l,HC,HU,SD,WT,D,SLI,LITS, BIGS, BSLB, EPS, 

S  NSIG,HAIFN,LSI) 


C 

C 
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C  COMPUTE  THE  (UNIFORM  SEBVICE- LEVEL)  POWEB  IPPBOYIHATION  POLICIES 
C  AND  COMP1RB  THEIB  ASSOCIATED  COSTS  WITH  THE  COSTS  ASSOCIATED  WITH 
C  THE  GKD  POLICIES 
C 

CALL  COHPAB (El ,K , I, HO, BU, SD, NT, BIGS, L ITS,T1TLE#LSI, AL) 

350  CONTINUE 
STOP 
END 
C 

SUBROUTINE  ST  SAMP (N I ,K , L, HO, HO , S D, D, HT,BSL,N,K  1,L1, HOl,HU1, 

&  SD1,D1,HT1,BSL1,SAHP,NUHST,STBA,NUME,  IT  EH) 

C 

C  STSAMP  :  STB ATIPIED  SAMPLING 
C 

C  THIS  SUBROUTINE  TAKES  A  STB  ATI FIED  SAMPLE  FROM  THE  INVENTORY 
C  SYSTEM  ACCORDING  TO  THE  POLLOHI NG  PARAHETEBS: 

C 

C  S  ASP 
C  NUHST 
C  STB  A  (  1) 

C  STB  A  (1+  1 ) 

C  NUMB  { 1) 

C  NUMB(I+1) 

C 
C 

C  FOE  A  GIVEN  STRATA  MITH  NN  ITEMS  F BOH  WHICH  H  ARE  TO  BE 
C  SAMPLED,  THE  ITEMS  NN/2H,  3NN/2M,  5NN/2H,  .  .  .  ,  (2H-1)NN/2H 
C  APB  TAKEN,  WBEFE  THESE  FRACTIONS  ARE  BOUNDED-OFP  TO  THE  NEAREST 
C  INTEGER 

C  THESE  ITEMS  ABE  DENOTED 
C 

C  ITEM  (1)  »  ITEM  (2)  ,  .  .  .,  ITEM  (H) 

C 

C  THE  CONSTRAINT  HEIGHTS  (WT) ,  THE  UNIT  HOLDING  COSTS  (HO), 

C  AND  THE  BEPLENISBHENT  SBT-UP  COSTS  (K) 

C  ARE  BECOMPUTED  TC  BEFLECT  THE  IDEA  OF  REPLACING  ITEMS  NOT  IN 
C  THE  SAMPLE  BY  THOSE  ITEMS  NEAR  THEM  IN  THE  SAMPLE  IIE,  THOSE 
C  IN  THE  SAME  ST  BATA) 

C  IN  PARTICULAR,  TRESE  COSTS  ARE  MULTIPLIED  BY  THE  TOTAL 
C  VALUE  OF  THE  STRATA  (VALUE  =  H*HU)  DIVIDED  BY  THE  TOTAL  VALUE 
C  OF  THE  SAMPLE  FRCH  THAT  STB  AT  A 

C  THE  HEIGHT  POR  EACH  ITEM  IN  A  GIVEN  STRATA  IS  THE  SAME,  BEING 
C  THE  TOTAL  HEIGHT  OF  THE  STRATA  DIVIDED  BY  THE  NUMBER  OF  ITEMS 
C  SAMPLED  FROM  THE  STRATA 
C 

REAL* 8  K  (512)  ,L  (512)  ,HO  (512)  ,MU  (512)  ,SD  (512)  ,HT  (512)  , 

6  D  (512)  ,  BS  L  (512) 

REAL*  8  K1  (512)  ,L1  (512)  ,H01  (512)  ,MU  1  (512)  ,SD1  (512)  , 

S  HT  1  (512) , D1  (512) , BSL1  (512) 

BEAL* 8  HEIGHT, RATIO, XM, XNSTR, HALF, ZERO, HWT,HSAHWT 


=  NUMBER  OF  ITEMS  IN  THE  SAMPLE 
=  NUMBER  OF  STRATA  +  1 
=  0 

=  LAST  ITEM  IN  STRATA  I,  I  =  1,  2,  .  .  .  ,  NUHST-1 

=  0 

=  BUHBBR  OF  ITEMS  IN  THE  SAMPLE  F BOH  STRATA  I, 

1  =  1,  2,  ...  ,  NUHST-1 


n  n  n  n 
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INTBGEB  Si  BP,  HUH  ST,  STUB  (512)  ,MUBB(512)  ,ITEB  (5 12)  ,  PI  BST,  L  AST 
IMTBGER  ■I,B,I,J,IJ,1JIJ,JJ,N0NST,B,STBAI,NSTB 
C 

DATA  HALF,ZEFO/0. 5DO,O.ODO/ 

C 

I  *»  S AS P 
IJ  =  0 

DO  60  I  *  2, BOSS! 

STBAI  =  ST B A  (1-1) 

B  *=  BOBB  (I) 

IS  =  H 

IF  B  =  0,  THEM  BC  IIEH  IK  THE  SASPLE  IS  TO  BE  TAKEN  FBOB  THIS 
STRATA 

IF  (B.EQ.O)  GO  TO  60 

■  STB  =  STB  A  (I) -SIFAI 
XM5TB  =  BST  B 
BATIO  =  IBSTB/IS 
HUT  =  ZEBO 
HEIGHT  =  ZEBC 
DO  20  J  =  1 ,  N STF. 

JJ  =  STFAI+J 
BUT  =  HWT+HO  (JJ)  *BU  (JJ) 

20  HEIGHT  =  HEIGHT+HT  (JJ) 

HEIGHT  =  HEIGBT/IB 
BSAMHT  =  ZEBC 
FIBST  =  IJ+1 
DO  40  J  =  1 , B 
IJ  =  IJ+1 

IJIJ  =  STBAI-)  (J-HALE)  *BATIO+HALF 
ITES  (IJ)  =  IJIJ 
FI  (I J)  =  K  (IJIJ) 

L 1  (I J)  *  L(IJIJ) 

HOI(IJ)  =  HO  (IJIJ) 

BUI  (I J)  =  HU  (IJIJ) 

SDI(IJ)  =  SD  (IJIJ) 

D  1  (I  J)  =  D  (IJIJ) 

BSLI(IJ)  =  BS  L  (IJIJ) 

HSASHT  =  HSASW1  +  H01  (IJ) *BU1 (IJ) 

HTI(IJ)  =  HEIGHT 
40  CONTINUE 

l  AST  =  PIFST-JB-I 
HUT  =  HUT/HSABHT 
DO  SO  J  =  FIBST, LAST 
HOI  ( J)  =  HUT  *HO  1  (J) 

50  K  1  (J)  =  HUT*K  1  (J) 

60  CONTINUE 
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C 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


RETURN 

END 


SOB  ROUTINE  BND5TS  (HI  ,K  ,  L  ,  HO,  BO,  S  D,  D,  NT  ,  BSL.  N,  K  1 ,  L 1 ,  HO  1,  HU  1 , 

6  SD1, D1,WT1,BS  LI, SAMP, NUH  ST, STRA, NUMB, ITEM, 

6  SEED) 

BNDSTS  :  BIN  DOM  STB ITIFIED  SAMPLING 


THIS  SOBFODTINE  TAKES  A  B  AN  DON  STRATIFIED  SAMPLE  FBOH  THE 
INVENTORY  STSTEM  ACCCB DING  TO  THE  FOLLOWING  PABAHETEBS: 


SAMP 

NUMST  = 

STB  A ( 1 )  = 
STB  A  (I-f- 1 )  = 
NUMP(I) 

NOME*  (1  +  1 )  * 

SEED  *  SEED 
BECD 


NUMBER  OF 
NUHBEE  C P 
0 

LIST  ITEM 
0 


ITEMS  IN  THE  SAMPLE 
STRATA  +  1 

IH  STRATA  I,  I  *  1,  2, 


NUMST- 1 


NDHBZF  CF  ITEMS  IN  THE  SAMPLE  FBCM  STBATA  I, 

I  =  1,  2,  .  .  .  ,  NUMST- 1 

FOB  PSEODO-BANDOM  N0N3EB  GENEBATOB  GGOD  (IHSL) 
IBE  SEED  IN  (I.ODO, 2147483674. 0D0) 


FOB  I  GIVEN  STBATA  WITH  NN  ITEMS  FBOM  WHICH  H  ARE  TO  BE 
SAMPLED,  THE  ITEMS  IBE  SAMPLED  BIN  DOflLY  ACCORDING  TO  A  DISCRETE 
UNIFORM  DISTRIBUTION 
THESE  ITEMS  ABE  DENOTED 


ITEM  (1)  ,  ITEM  (2)  ,  .  .  ITEM  (N) 


THE  CC1 STB  A I  NT  HEIGHTS  (WT) ,  THE  UNIT  HOLDING  COSTS  (HO), 

AND  THE  REPLENISHMENT  SET-UP  COSTS  (K) 

ARE  RECOMPUTED  TC  REFLECT  THE  IDEA  OF  REPLACING  ITEMS  NOT  IN 
TfcB  SAMPLE  BY  THOSE  ITEMS  NEIB  THEM  IN  THE  SIMPLE  (IE,  THOSE 
IN  THE  SAME  STBATA) 

IN  PARTICULAR,  THESE  COSTS  ARE  MULTIPLIED  BY  THE  TOTAL 
VALUE  OF  THE  STBATA  (VALUE  =  H*HU)  DIVIDED  BY  THE  TOTAL  VALUE 
OF  THE  SAMPLE  FBCM  THAT  STBATA 

THE  WEIGHT  FOB  EACH  ITEM  IN  A  GIVEN  STBATA  IS  THE  SAME,  BEING 
THE  TOTAL  WEIGHT  CF  THE  STRATA  DIVIDED  BY  THE  NUMBER  OF  ITEMS 
SAMPLED  FBOH  THE  STRATA 


R  EAL*8 

6 

RBAL»8 

& 

R  BAL*  8 


K  (512)  ,L  (512)  ,HO  (512)  ,HU  (512)  ,SD(512)  ,NT  (512)  , 

D  (512)  , ES  L  (512) 

K  1(512)  ,L  1(512)  ,HC  1(512)  ,HU1  (512)  ,SD1  (512)  , 

NT  1 (512) , C 1 (512) , BSL 1  (512) 

WEIGHT, PATIO, XM,YNSTR, ZERO, HWT,HSAMWT, SEE D, DSE ED 


I  NT  EGER  SAMP,  NUMST,  STR  A  (512),NUMB(512),ITEM(512),FIRST/LAST 
INTEGER  NI,N,I,J,IJ,IJIJ,JJ, NO HST, H, STB  A I , NSTfl , IN  DEL ,ON E, 

(  I  E  (  1) 


n  n  n  n 
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C 

DATA  XEBO,OIE/0. 0D0# 1/ 
C 

DSEED  =  SEED 
M  *  SAHP 
1J  =  0 

DO  60  I  *  2,MUHST 
STP  k  I  =  SIBA(I-ONE) 
fi  =  NOHB  (I) 

III  *  H 


IF  B  =  0,  TBEN  NC  ITEB  IN  THE  SAMPLE  IS  TO  BE  TAKEN  FROM  THIS 
STB ATA 


IF  (B.EQ.O)  GC  TO  60 
C 

BSTB  =  STB  A  (I) -STBAI 

XMSTR  =  MSTB 

EATIO  =  INSTB/XB 

BBT  *  ZERO 

¥BIGHT  =  ZEBO 

DO  20  J  =  ONE , MSTB 

JJ  *  STBAI-fJ 

BBT  *  HWT+HO  (JJ)  *HU  (JJ) 

20  BRIGHT  *  MElGHI-fNT  (JJ) 

WEIGHT  =  WEIGHT/XB 
HSANk'T  *  ZEBO 
FIBST  =  IJ+OHE 
DO  40  J  =  C  NE  ,  H 

25  CALL  GGUD  (DSEED, N3TB,ONE#IR) 

IJIJ  =  STEAI-)IB  (ONE) 

DO  30  INDEX  =  CNE , I J 

IP  (IJIJ. EC. ITEB  (INDEX) )  GO  TO  25 

30  CONTINUE 

IJ  =  IJ-fO  ME 
ITEB  (IJ)  *  IJIJ 
K1(IJ)  =  K  (IJIJ) 

L  1  (I  J)  =  L(IJIJ) 

HOI (I J)  =  BO (IJIJ) 

BOI(IJ)  =  BU  (IJIJ) 

S D 1  (IJ)  =  SD  (IJIJ) 

D  1  (I  J)  =  D  (IJIJ) 

ESLI(IJ)  =  BSL(IJIJ) 

HSABWT  =  HSAH  AT-)  HCl  (IJ)  »BU1  (IJ) 
WTI(IJ)  *  WEIGHT 

40  CONTINUE 

LAST  =  FIRST-)  B-CNE 
HHT  =  HWT/HSABWT 
DO  50  J  *  FIBST, LAST 
HOI (J)  =  HWT* HCl  (J) 
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50  K1  (J)  «=  HNT*IM  (J) 

60  CONTI  HOE 
RETURN 
END 
C 

SOB BOOT INI  El GS  LS  (N 1 , L , HO , H 0 , SD , D , NT , BSL , LSI , EPS , NSIG ,HAX  FN, 
S  BIGS , LITS) 

C 

C  THIS  SUBROUTINE  COMPUTES  THE  OPTIMAL  (S,S)  POLICIES  ASSOCIATED 
C  fclTH  THE  KUHN-TOCKER  CCNSTBAINT  MULTIPLIER  LAMBDA-STAR  (LSI) 

C 

c 

COMMON  /FINDTS/AI,BI,EI,FI,GI,HI,LAH 
C 

REAL*8  LITS  (512)  ,BIGS(512) 

B  B  AL*8  L  (51 2)  ,  BO  (51 2),  MU  (512)  ,SD(512)  ,  NT  (512)  , 

S  D(512)  ,ESl  (512) 

BEAL* 8  EPS, LSI, LI, HOI, HU  I , SDI , D I , MSI , NTI ,TSI , L AH, 

S  BHI,QI,GAI,DEI, II , BI , El, FI , GI , HI , V IB  1 , V AR2 , ET I ,TU I , 

S  KSI,LPTDI,BSLI, X,  TUI,  LB,  UB, ZERO, XONL, HALF, CONST 

C 

INTEGER  NI,NSIG,HAXFN,HAX,ONE,CNTR 

C 

REAL*8  LAMTH 
EXTERNAL  LAMTH 
C 

ZERO  =  0, 

ONE  =  1 
XONE  =  1. 

HALF  =  .5 

CONST  =  1.700436904 
C 

DO  60  I  =  ONE,NI 
HOI  =  HO  (I) 

M 01  =  HD  (I) 

LI  =  L(I) 

MSI  =  LI*HUI 
SDI  =  SD  (I) 

W T I  =  NT  (I) 

DI  =  D  (I) 

BHI  =  MUI/(DH  ((HUIi(SDI*SDI/HUI))*HALF)) 

QI  =  XONE-RHI 

GII  =  CONST/ (LI*SDI) 

DEI  =  DEX  F  (G  A I *DI ) 

A I  =  D I*  HOI 
V  A B 1  =  XCNE+DEI 
BI  =  AI* (QI+VAR1) 

El  =  AI* (DEI4Q1*VAR1) 

FI  =  AI*QI*DEI 
VAB1  =  D1*RHI*GAI 


noon 


VAF2  =  QI*  (DEI-XCNE) 

GI  *  MTI  •  (VAB1-»VAB2) 

HI  *  *TI*  (V  AB  1*DEI*  VAB2) 
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C 

C  DETERHI NE  BSL  BI  THE  LABGEST  f ALOE  OF  THETA,  TO,  AHD  HENCE  THE 
C  SHAILEST  VALUE  OF  BIG  S,  SO  THAT  D  (L IHBD A) /D  (THET A)  LE  0, 

C  GUARANTEEING  THAT  L A S3 C A (THETA)  IS  INVEBTIBLE  OH  (0,  TO) 

C  IF  THE  0 SEE- S0PP1I E D  BSL  IS  GBEATEB  THAN  THIS,  IT  BBCOHES  THE 

C  BSL  OSED 
C 

ETI  =  GI*  (BI*HI*HI-GI* (EI*BI-FI*GI) )“AI*HI*HI*HI 
C 

C  IF  ETA (I)  IS  LESS  THAN  0,  LANBDA (THET A)  HAT  NOT  BE  CONVEX,  AND 
C  THEBE  IS  HO  GUABANTEE  THAT  THE  GKD  ALGOEITHH  BILL  CONVEEGE  TO 

C  A  LOCAL  HIN1HUH - THE  USEE  SHOULD  EXEECISE  CABE 

C 

IF  (ETI.  LT.  ZEBC)  1BITE(3,30)  I, ETI 
30  FOBS  AT  (  '  ETA  ('  ,16,  *  )  =*,E20.6,'  LT  O') 

C 

CMTF  =  -1 

35  CHTB  *  CITB+CNE 

TUI  =  (DEXF (GAI*  (HSI-CNTE) ) ) /DEI 
ESI  =  (FI/HI)  -  (AI/(GI*TUI*TOI)) 

T  AE  1  =  GI+HI*TUI 

LPTUI  =  KSI- (ETI/(GI*BI*VAR1*VAF1) ) 

IF  (LPTOI.GE.  ZERO)  GO  TO  35 
C 

BSLI  =  CNTB+DI 
X  =  BSL  (I ) 

BSL  (I)  =  ESL1 
IF  (X.GT.BSL1)  BSL  (I)  *  I 
TUI  =  (DEXP  (GAI*  (HSI-  (BSL  (I)-DI)  )  )  )  /DEI 
C 

LAH  =  LSI 
H  AX  =  HAXFN 

LB  =  ZEBO 
03  =  TO I 
C 

C  IF  LAHTH(UB)  IS  GEEATEB  THAN  ZERO,  THEN  THE  FUNCTION 
C  L  AifcDA  (THETA)  FOB  THIS  IT  EH  IS  ALWAYS  GREATER  THAN  LAHBDA-STAF, 
C  AND  SC  BIGS  POB  THIS  1TEH  IS  SET  AT  ITS  LOWEB-BOUND  BSL 
C 

IP  (LAHTH  (UB)  ,GT.  ZEBO)  GO  TO  50 
CALL  ZBKEKT  (LAHTH,EPS,NSIG,LB,OB,HAX,IER) 

TSI  «=  UB 

BIGS  (I)  =  HS1- (  (DLOG  ITS1)  ) /GAI) 

BECAUSE  OF  FOUNE-OPP  ERROR ,  IT  HAT  BE  THAT  BIG  S  LT  BSL. 

IF  SO,  SET  BIG  S  =  BSL. 
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IP  (BIGS(I).IT.BSL(I))  BIGS(I)  *  BSL(I) 

LITS(I)  *  BIGS  (I )  -DI 
GO  TO  60 

50  BIGS  (I)  «  BSL(I) 

UTS  (I)  '  BIGS  (I )  -DI 

60  CONTINUE 
RETURN 

END 

BEAL  FUNCTION  LANTH*8  (THETA) 

C 

COMMON  /FINDTS/AI, B1,EI, FI, GI,H1,LAM 
RE AL*8  AI,BI,EI,FI,GI,HI,LIM,THETA 

LIHTH  =  AI-fTHETA*  (B1-LAM*GI+THETA*  (EI-LAH*HI-fTHETA*PI)  ) 

RETURN 

END 

SUBBOUT INE  CCMFAB(NI,K,L,H0,HU,SD,HT,BIGS,LITS,TITLE,L5I,AL) 

C 

c  THIS  SOBBOUTINE  DOES  THE  FOLLOWING: 

C  (1)  EVALUATES  EXACTLY  THE  HOLDING  AND  F EPLEN IS H HE  NT  COST,  AND 
C  THE  SERVICE-LEVEL  FOR  THE  POLICIES  (LITS,BIGS)  UNDER  THE 

C  ASSUMPTION  OF  N EG AT IV E- BINOMIAL  DEMAND  (IN  EXACT) 

C  (2)  COMPUTES  APPROXIMATELY  OPTIMAL  (S,S)  POLICIES  WITH  UNIPORM 
C  SEB  VICE-  LEVELS  (IN  POWAPP) 

C  (3)  EVALOATBS  THE  POLICIES  DESCRIBED  IN  (2)  EXACTLY  AS  IN  (1) 

C 

C  (4)  PRINTS  THE  RESULTS  OP  STEPS  (1)  AND  (3) 

C 

REAL*8  LITS  (512)  ,PIGS(512) 

PEAL* 3  K  ( 5 1 2)  , L  (5 12)  #HO  (512)  ,HU  (512)  , SD  (512)  , NT  (5  12) 

BEAL* 8  AL,HGFE,HPO*,TGKD, TPOH , LSI , LITS I , BIGS  I 
B  EAL*8  KI , LI , BOI , MUI ,SDI , SLI , HCI , BCI , AVSL, PAVSL 
C 

ISTBGER  NI  ,TITLE  (8) 

AVSL  =  0, 

HGKD  =  0. 

TGKD  =*  0, 

DO  70  I  *  1 , NI 
AI  =  K(I) 

LI  =  L ( I ) 

HOI  =  H 0  (I ) 

MUI  =  MU(I) 

S DI  =  SD  (I ) 

BIGSI  =  BIGS  (I) 

LITSI  =  LITS  (I) 

CALL  EXACT (Ml, LI, HOI, MUI ,SDI, LITSI, BIGSI, SLI, HCI, RCI, IER) 
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IF  (IBB.  EC*  1)  GO  TO  300 
A  YSL  =  A?5L+SL1*WT (1) 
H3KD  *  HGFD+HCI 
TGKD  *  TGKC+HCI+BCl 
70  CONTliOE 


C 


c 

c 

c 


c 


c 


PAYS  L  =  0. 

HPOB  =  0. 

TPOW  =  0. 

DO  75  I  =  1,NI 

A I  *  K  (I ) 

LI  =  MD 

HOI  =  HO  (I) 

HOI  =  HU  (I) 

SOI  *  SO  Cl) 

CALL  POtiPPfKl, Ll,HCI,BOI,Si>l,AYSl,LITSI, BIGSI) 

CALL  EXACT (Kl,IZ,HCi«HUl , SOI , LITSI , BI GS I , SL I , HCI, BCI , IEE ) 
IF  (I  EE. EC. D  GO  TO  300 
PAYSL  =  PAVSL+SLIFHT  (I) 

H POl  =  HP0S*HCI 
TPCi  *  TPOW4HCHBCI 
75  CO  NT  IStl  E 


PRINT  SUBMABY 


NBITE(3,98)  (TITLE(J)  ,J=1,8)  ,LSI 
98  F3FBAT(/T2,8A4,'  LAMBDA  ST 

WRITE(3,100)  Al,  AVSL,PAVSL,HGKD, 
100  FOFBAT(T3,' SPECIFIED  ALPHA  = 

5  T3, 'ACTUAL  ALPHA  (GKD)  =' 

6  T3,  'ACTUAL  ALPHA  (PON)  = 

&  T3, 'HOLDING  COST  (GKD)  =' 

5  T3,  •  BCLDING  COST  (POH)  = 

6  T3,  '  TOTAL  COST  (GKD)  =' 

S  T3, 'TOTAL  COST  (PON) 

BBTUBB 


AH  =  '  ,  F 1 
HPOH  ,  TGKD 
• ,T27,F1 1 
,T27 ,  F 1 1 • 
• ,T27, F1 1 
#T27, F 11. 
• ,T27,F1 1 
,T27,  F  1 1 . 
•,T27,F11 


2.4,//) 

,TPOH 

.4,//, 

4,//, 

•  4 , // , 

4,//, 

.4,//, 

4,//, 

.4) 


300  WHITE  (3 , 3  10)  I 

310  FOBH  AT ('  BIG  S  OR  D  >  2000  FOE 
G  •  E  X  ACT  CANNOT  EVALUATE 


ITEM* , 16,/ 

THE  (S,S)  POLICY',//) 


EETORN 

END 


C 

SOBHOUTINE  EXACT  (K , L, H , HU , SD, LS , BS , SL, HC, BC, I E B) 

C 

C  THIS  SUBROUTINE  COMPUTES  THE  SERVICE-LEVEL,  HOLDING  COST  AND 
C  REPLENISHMENT  COST  EXACTLY  FOR  THE  POLICY  (LS,BS)  UNDER  THE 
C  ASSUMPTION  OF  NEGATIVE  BINOMIAL  DEMAND. 

C 


172 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


K  *  SET-OP  COST 

L  «=  LEIDTIHE  4  1 

H  *  OMIT  HOLDING  COST 

HU  *  BEAN  DEBAND 

SD  *  DEBAND  STANEIBD  DEVIATION 

LS  =*  LITTLB  S 

BS  =  BIG  S 

SL  *  STEADI-STATE  S EB V ICE- LEVEL  (FREQUENCY  OF  PERIODS  THAT 
NO  BACKOBDEB  IS  PLACED) 

HC  =  EIPECTED  HOLDING  COST 

FC  =  EXPECTED  B EPL E N IS HB ENT  COST 


BEAL*  8  K  ,1 , H , HD, SD,  LS,  BS  ,SL,  HC,  BC,  P,  Q,  B,  BS,PROB  (2000)  .HALF, 

&  NB  (200  0)  ,  fl  (2000)  ,  HI ,  BHO,  V  AH  ,SUfl,  PEOBI  ,  HD,  ZE30,  XONZ,  IT  HO 

INTEGEB  S,D,CE  1,032, ONE, TWO, IBR 

DATA  ONE, TWO, ZERO, XONE,XTWO,HALF/1, 2,0. DO , 1 . DO , 2 . DO , . 5D 0/ 

S  AND  D  ABE  BOUNDED  TO  THE  NEAREST  INTEGER.  HE  USE  THE  REORDER 
QUANTITY  D  =  BIG  S  -  LITTLE  S  -  1  BECAUSE  HE  OSE  THE  VERSION  OF 
AN  (S,S)  POLICY  THAT  BEQD1BES  THAT  AN  ORDER  BE  PLACED  WHEN  THE 
INVENTORY  POSITION  DROPS  STRICTLY  BELOH  LITTLE  S. 


PHOB(I)  =  PBOB  (NE  RV  =  1-1) 

N B  (I)  =  PBOB  (NB  4  .  .  .  4  NB  RV  =  1-1) 

- 

H(I)  =  PBOB (RENEWAL  PNC  OF  NB  BV  =  1-1) 

IEB  =  0  EXACT  NAS  ABLE  TO  EVALUATE  (S,S)  POLICY 
1  OTHERHISE 


IEB  =  0 
S  =  BS4HALF 
D  =  (BS-LS-ONE)  4HALF 
Q  =  HU/(SD*SD) 

P  =  XONE-C 
R  =  HU*  Q/P 
BS  =  L*R 
VAR  =  Q**B 
PROB(ONE)  = 

NB(ONB) 

RHO 
H  (ONE) 

UB1  =  S+ONE 
IP  (D.GT.S) 

IP  (UB1.  GT.  2000)  GO 
DO  20  I  =  THO  ,UB  1 
VAR  =  IONE/ (I-XONE) 

PBOB  I  =  PROB  (I-ONE)  *T*  (B-)I-XTHO)  *V  AR 
PROB(I)  «  PBOBI 


VAR 

C**RS 

XCNE/  (XCBE-V AB) 
V  Afc*BHO 


DB1  = 


D40NE 
TO  100 


nnn  nnnnnn 
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NB  (I )  =  NB(I-CNE)  *P*  (BS-f  I-XTMOJ  *VAR 
SUfl  *  ZERO 
032  =  I-ONE 
DO  10  J  »  CNE , UB  2 

io  sun  *  soh+pbob  (I-J40HE)  *a  (J) 

H (I )  =  (PBOBI-fSUP)  *BHO 
20  CONTINUE 

OBI  =  D+ONE 

IP  (OBI.  Cl.  2000)  GO  TO  100 
BD  =  ZERO 
DO  30  I  *  CNE , OB  1 
30  PD  =  BD+H  (I) 

FBO  =  XONE/ (XCNE-fHD) 

COPPDTE  EXPECTED  B EFLENISHB ENT  COST 

BC  =  K*BH0 

COBPUTE  EXPECTED  BCLDIBG  COST 

IP  (D.GT.S)  OE 1  =  S+ONE 

IP  (OBI.  GT.  2000)  GO  TO  100 

HC  =  ZERO 

DO  40  I  =  ONE # OB  1 

OB2  =  S-I4TMC 

PI  =  If  (X) 

DO  40  J  c  CBE,UB2 
HC  =  BC+ (0B2-J)  *NB  (J)  *BI 
40  CONTIIOE 

SUP  =  ZEBO 
DB2  =  S+ONE 
DO  45  I  =  ONE.OB2 
45  SUB  =  SUH-f  (OB 2- 1)  *NB  (I) 

HC  =  H* (HC+SUE) *  BHO 

COPPDTE  (STEADI-STATE)  SERVICE-LEVEL 

SL  =  ZEBO 
DO  60  1  =  ONE, OBI 
SUP  *=  ZEBO 
UB2  =  S-I+TNC 
DO  50  J  *  CNE, 0B2 
50  SUP  =  SUfl-(NB(J) 

60  SL  =  SLfSUP*P  (I) 

SUP  =  ZEBO 
OBI  *  S f C HE 

IP  (UB1. GT. 2000)  GO  TO  100 
DO  70  I  =  ONE, OBI 
70  SUP  *=  SUH4NB  (I) 

SL  =  (SL-fSUH)  *RH0 


n  n 
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C 

100 


c 


RETURN 

3  EF  =  1 
RETURN 

END 


SUBROUTINE  PCWAPP  (K , L, H, HU, S D, ALP HA, LITS , BIGS ) 

C  THIS  SUBROUTINE  COMPUTES  LITTLE  S  AND  BIG  S  USING  THE  REVISED 

C  POWER  APPBOXIHATION  HETHOD  (HOSIER ( 198 1) ,  TECH  REPORT  *18). 

C  T  BE  SHORTAGE  COST  IS  COBPUTED  USING  THE  POWER  APPROXI H ATION 

C  HETHCD  (EHBHAPDT  (1977,  PP  18,  45),  TECH  BEPORT  *12). 

C 

c  K  «=  SET-UP  COST 

C  L  *=  LEADTIBE  +  1 

C  H  =  UNIT  HOLDING  COST 
C  HU  3  BEAN  DEHAND 

C  SD  =  BEAN  STANDARD  DEVIATION 

C  ALPHA  =  (STEADY-STATE)  SERVICE- LEVEL,  THE  FREQUENCY  OF  PERIODS 
C  WITHOUT  E ACKLCGS 

C  I ITS  =  LITTLE  S 

C  EIGS  =  BIG  S 

C 

PEAL* 8  P, H, K, HU, D,L,VL,SrL,Z,V,SD, LITS, BIGS, ONE, 

S  C1,C2,C3,C4,C5,C6,C7,C8,C9, CIO, ALPHA 

C 

DATA  Cl/1.3/ 

DATA  C2/.494/ 

DATA  C3/.506/ 

DATA  C4/.  116/ 

DATA  CS/. 973/ 

DATA  C6/. 183/ 

DATA  C7/1.063/ 

DATA  C8/2.192/ 

DATA  CIO/. 0695/ 

DATA  ONE/1./ 

C 

V  =  SD*SD 

VL  =  V*  L 

SDL  =  DSQBT  (VL) 

C 

P  =  H*  (ALPUA-C10) /(CNE-ALPHA) 

Z  =  Ch'E)  (VL/  (HU*HU)  ) 

D  =  Cl*  (HU*  *C2)  *  (  (K/H)  **C3)  *  (Z**C4) 

Z  =  DSQBT  (D*  11/  (P* SDL)  ) 

C9  =  (C6/Z)  -fC7-C8*Z 

LITS  =  C5*L*HU-J5DL*C9 
BIGS  *  LITS-)D 
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RITUFN 

END 

C 

SUBROUTINE  GKD(Nl,L,HO,HU,SD,W,DD, AL, LITS , BIGS #BSL, EP # NS , HI , LS) 
C 

C  GKD  (GENERALIZED  KNAPSACK  DUALITY)  ALGORITHM 
C 

C  NOTE  :  THIS  SUBBCUTINE  CHANGES  THE  VALUES  OP  THE  INPUT 
C  PARAMETER  BSL 

C 

C  NT  =  HEIGHTS  (BUST  BE  CONVEX) 

C  RB  =  BHO 

C  GA  =  GABMA 

C  DE  =  DELTA 

C  A  =  A 

C  B  =  B 

C  E  =  E 

C  F  =  F 

C  G  =  G 

C  B  =  H 

C  MS  *  BU-STAB  =  L  *  HU 

C  ET  =  ETA 

C  KS  =  KSI 

C  TU  =  THETA-UPPEP 

C  LPTO  *  LAMBDA- PFlflE  (  THETA-UPPER  ) 

C  LL  *  LAMBDA-LOWER 
C  MLL  =  MAI{  LAHEEA-LCHEB  ) 

C  FBLL  =  P(  H LL  ) 

C  LU  =  LAMBDA-UPPER 

C  H  LU  =  MAY(  L ABBC  A-UPPER  ) 

C  T S  *  THETA-STAB 

C  LS  =  LAMBDA- STAR 

C  BSL  =  LOWER  BOUND  CN  BIG  S 
C  MSL(I)  =  SEEVICE-LEVEL  POR  ITEB  I 
C  CONST  =  15  •  PI  /  (16  ♦  SQET  (3) ) 

C  NI  *  NUMBER  OP  ITEMS  IN  THE  INVENTORY  SISTEH 
C  N  (IT  PROGRAM  BEGINNING)  =  NI 

C  N  (AT  PROGRAM  TERMINATION)  =  NUMBER  OF  ITEMS  WITH  BIG  S  GREATER 

C  THAN  BSL 

C  PERU  :  PERM  (1)  v  .  .  .  ,  PE F H  (N)  ARE  THE  ITEMS  MITH  BIG  S 

C  GREATER  THAN  BSL.  WHFN  NI  IS  GREATER  THAN  *, 

C  PERM  (N-f  1 )  ,  .  .  .  ,  PERM  (K I)  ARE  THE  ITEMS  NITH 

C  BIG  S  =  BSL 

C 

COMMON  /INPUT/  NT, ALPHA, EPS, NSIG,HAXFN,N 
COMMON  /LAM/  A , B , E , F,G , H 

COMMON  /PARAMS/  T U, T S, R H , DE, G A , D , XO N E , Z E BO, PEBfl , ON E 
COMMON  /S  U  BR /  LAM, I 
C 

REAL*  8  LITS  (512)  ,BIGS(S12)  ,BSL(512)  ,N  (512), DD  (512) 
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BEIL* a  L  (512)  ,H0  (512)  ,110  (512)  ,SD  i S T 2 J  , MT (512)  ,  ALPHA 
B  EAL*B  D  (512)  ,»H  (512)  ,GA  (512)  ,DE  (512)  ,  A  (512)  ,  B  (51  2)  ,E  (512)  ,r  (512)  , 
*  C  (512)  ,U  (512)  ,TQ  (512)  ,BS  (512) ,  LL  (512)  ,1,0  (512)  .8SL1512)  , 

C  TS  (512)  ,BLi,|BLL,HLU,LS,ZEROfJUL?,EPS,EP,AL 
BEAL* 9  CONST,  HOI  ,!IUI, LI, SD1,UTI,  AL 1 , BI , E I , WTAT5L, 
t  PI, BUI,  Q1,GA1,DE1,11,  B1  ,EI,ri,CJ,HI,fAF1,VAR2,TUI,TS3,LU, 

C  LPTUI, BSI,IAB,BIGSI,THS,TE5T,SUB,I, 

«  LP  ,  OB,  R  5 1 ,  El  I ,  B5LI  ,TL  I ,  LLI ,  LU1  #  IQNE. 

1  IT  EGER  >,OIE,M,CJtTB,  PE  RH  (512)  , ILB, J , IN D El , B, II 
IHTSGEB  iS,lSIG,BMFN,BAI 

PE  A  1*3  LAHPDA,F,ITBL 
EXTERN A L  F , LT H L 

ZEBO  *  0. 

ONE  *1 
XONE  =  1. 

NILE  =  .5 

CONST  =  1.7001436904 

ALPHA  =  IL 
EPS  =  EP 
ISIG  =  IS 

BAXPN  =  HX 
V  =  NI 

DO  3  I  =  ONE , HI 
3  NT  (I)  =  W (I) 

WHITE  (3,5) 

5  FORK  AT  ( *  1 ' ) 

INITIALIZE  THE  PERBOT AT10N  BATRIX  PERU 

DO  6  I  -  CNE.N1 

6  PERU  (I)  =  I 

INITIALLY  SET  LAHBDA-STAR  =  0,  SINCE  THIS  IS  HOW  TO  INTERPRET 
LASEDA-STAB  IP  A  LI  POLICIES  ABE  SUCH  THAT  BIGS  =  BSL 

LS  =*  ZEBO 


COBPUTE 

THE  CONSTANTS 

DO  2  0 

I  =  CNE, NI 

HOI  = 

HO  (I ) 

BUI  = 

BD  (I) 

LI  = 

L  ( I ) 

B  SI  = 

LI*  BOI 

ns(i) 

=  BSI 
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SDI  =  SD  (I ) 

BTI  =  NT  (I) 

DI  *  DD  (I ) 

D(I)  *  DI 

BUI  *  HDI/(DI-(  (  (H U1-)  (SDI*SDI/HUI))*HALF)  ) 

BB  (I)  *  BUI 

Cl  *  IONE-BHI 

Gil  *=  COB  ST/  (LI*SDI) 

G A  (I )  -  Gil 

DEI  =  DEIP  (GII*DI) 

DE (I)  =  DEI 
AI  =  DI* HOI 
i(I>  *  hi 

¥  IB  1  *  ICNE+DEI 
BI  *  il*  (Cl*ViBl) 

B  (I)  -  BI 

El  *  il*  (DEIKI'ViBI) 

E  ( I)  =*  El 

PI  *  iI*QI*DEI 

F (I)  *  PI 

VAF1  =  DI*BHI*GAI 

ViF2  =  Cl*  (DEI-XCBE) 

GI  =  1TI*  (VABl-fVAB2) 

G  (I>  *  GI 

HI  =  UTI*  (VAB1*DEI+VAB2) 

H  (I)  =  HI 
C 

C  DETEBHINE  BSL  BY  THE  LiBGEST  ViLOE  OF  THETA,  TU,  AND  HENCE  THE 
C  SflALLEST  VALUE  OF  BIG  S,  SO  THAT  D  (L  AH  BD  A) /D  (TH  ET  A)  LE  0, 

C  GUIP  INTEBING  THAT  L  A  HBD  A  (TH  ET  A)  IS  INVEBTIBLE  ON  (0,  TO) 

C  IP  THE  USEB-SOPPLIED  BSL  IS  GBEATER  THAN  THIS,  IT  BECOMES  THE 
C  ESL  USED 
C 

ETI  =  GI* (BI*HI*HI-GI*  (EI*HI-FI*G1) ) -AI*HI*HI*HI 
C 

C  IF  ETA(I)  IS  LESS  THAN  0,  LAMBDA  (THETA)  BAY  NOT  BE  CONVEI,  AND 
C  THERE  IS  NO  GUARANTEE  THAT  THE  GKD  ALGORITHM  BILL  CONVERGE  TO 

C  A  LOCAL  MINIMUM - THE  USES  SHOULD  EXERCISE  CABE 

C 

IP  (ETI.  LT.  ZEFC)  NRITE(3,10)  I,  ETI 
10  FORMAT  (*  ETA(',I6,')  =  ',E20.6,»  LT  O') 

C 

CITE  *  -1 

7  CMTB  =  CNTR-fCNE 

TUI  =  (DEXP  (GAI*  (MSI-CNTP)  )  ) /DEI 
KSI  =  (PI/HI) - (AI/ (GI*TUI*TUI) ) 

V  AR  1  *  GI-)  H  1*TUI 

LPTUI  =  KSI- (ETI/ (GI*HI*VAB1*¥AF1) ) 

IP  (LPTUI. GE. ZERO)  GO  TO  7 
C 


n  n  n 
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BSLI  *  CNTH+DI 
X  *=  BSL (I) 

BS1 (I )  *  BSLI 
IT  (I.GT.BSLI)  BSL(I)  «=  X 
TUI  =  (DEXP  (GAI*  (HSI-(B5L  (I)  -Dl)  ))  )  /DEI 
TO(I)  *=  TUI 
C 

LLI  =  LAMBDA  (TUI  ,  I) 

LL  (I)  =  LLI 
20  CONTINUE 
C 

C  COMPUTE  HLL  =  BAX  (LABBDA-LOBEB)  =  H AX  (LAMBDA  (TU) ) ,  SO  ALL 
C  FUNCTIONS  LAHBIA  (THETA)  ABE  INVERTIBLE  ON  [ HLL , I NFI HITT ) 

C 

22  BLL  «  ZERO 

DO  30  II  =  ONE  ,  N 
I  «=  PERM  (II ) 

LLI  =  LL  (I ) 

IP  (BLL. LI. LLI)  BLL  =  LLI 
30  CONTINUE 
C 

PULL  =  P(HLL) 

C 

C  P(LAKBDA)  IS  INCRESING  ON  [  BLL,  IN  FINITT )  BITH  P(INFINITr) 

C  GBEATBB  THAN  0.  CHECK  IF  P (ALL)  fS  IFSS  THAN  OB  TOO *L  Tn 

c  o,  *0  nr*’-'  r (7.%-rrrj  -  *  r<r  r*rrrr  -  •**  r"  [ ti r ,rrrrr — r; 

c 

IF  (PULL. LL. NEMO)  GO  To  50 
C 

C  P  (BLL)  IS  GBEATEB  THAN  0.  LOCATE  THAT  ITEB  ILB  =  PERM  (INDEX) 

C  BITH  THE  LARGEST  VALUE  OF  LL  =  LAMBDA  LOWER,  AND  SET 

C  BIG  S  =  BSL,  IE,  SET  IS(=  THETA-STAB)  =  TU  (=  THETA-UPPER) 

C 

INDEX  =  ONE 

ILB  ^  PERM  (INDEX) 

DO  42  I  =  C BE , N 
J  =  PERM  (I) 

IF  (LL  (J)  .GT.  LL  (ILB)  )  INDEX  *  I 
ILB  =  PEBM  (INDEX) 

42  CONTINUE 
C 

C  CHANGE  N  TO  N  -  1,  AND  BEHOVE  ITEM  ILB  FROM  THE  FIRST  N  ENTRIES 

C  OP  PERH  AND  PUT  IT  IN  PERH(N+1) 

C 

V  =  N-ONE 

IF  N  *  0,  ALL  ITEMS  ARE  SUCH  THAT  BIG  S  =  BSL 

IF  (N.  EQ.  0)  GC  TO  85 
DO  43  I  =  INDEX, N 
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43  PEEK  (I)  * 

C 

PBBH  (B+  1 ) 
SET  IS  *  10 


PEFH  (141) 

-  ILB 

FOB  IT EB  ILB 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 


IS  (ILB)  *  TO  (ILB) 

DETERHINE  THE  SEBV1 CE-LEVEL  ALI  FOB  IT  EH  ILB,  END  BECOBPOT  E  THE 
A  V  EE  EC  E  SEBTICE-LEVEL  (WHICH  IS  AG  LI  N  DENOTED  ALPHA)  FOB  THE 
OTBEB  ITEBS  SO  IB  AT  THE  OVEBALL  AV  EB AGE  S EBV ICE-LEVEL  IS  STILL 
THE  OBIGIN  ALLY  SPECIFIED  ALPHA 

BOTE  :  IHEB  BIG  S  IS  SET  TO  BSL  FOB  THE  ITEH  ILB,  THE  ITEH 

IS  MO  LONGER  OSED  IM  THE  COHPUTATIOH  OF  THE  SOLUTION  LS 
TO  P  (LA  BEE  A)  =  0.  BATHER,  THIS  COHPUTATIOH  IS  DONE  ONLY 
FOB  THOSE  SUHH AMDS  OF  P  THAT  CORRESPOND  TO  ITEBS  BITH 
BIG  S  GREATEB  THAN  BSL.  HENCE  BE  MEED  ONLY  DETEBHIBE 
THE  AVEBAGE  SERVICE-LEVEL  FOB  THESE  ITEBS  SO  THAT  THE 
OVERALL  SFBV ICE-LEVEL  IS  THE  SPECIFIED  ALPHA 

DEI  =  DE (ILB) 

G  A I  *  G  A  (ILB) 

T5I  =  TS  (ILB) 

B HI  =  B H  (ILB) 

01  =  IONE-BRI 
DI  *  D  (ILB) 

BTI  **  HT  (ILB) 

A  LI  «  QI*  (XCME+DLOG  (  (XONE-fTSI)/ (XONE+ 

&  DEI*TS1)  )  /  (GAI*DI)  )4  (RH1/ (XONE+TSI)  ) 

ALPHA  «  (ALPHA-HTI*ALI)/ (XOME-BTI) 

IF  (ALPBA.LT. IONE)  GO  TO  22 
WRITE  (3,47) 

47  FORH AT ('  THE  MEW  AVEBAGE  SERVICE- LEVEL  IS  NOT  LESS  TRAN*,//, 
G  •  ONE ,  SO  THE  PBOBLEH  IS  INFEASIBLE',//) 

GO  TO  82 

50  VAH1  =  -PBLL 


C 

C  SINCE  P(HLL)  LE  0,  SOLVE  P(LAHBDA)  =  0  FOR  LAHBDA  =  LS 

C  01  [ BLL, INFINITY) .  FIRST  CHECK  IF  P(HLL)  =  0  WITHIN  EPSILON; 

C  IF  MOT,  DETERBINE  BLU  LESS  THAN  INFINITY  SO  THAT  P(BLU)  IS 

C  GREATEB  THAN  OF  EQU AL  TO  0.  THEN  SOLVE  P  [LAHBDA)  =  0  FOB 

C  LAHBDA  =  LS  OM  (BLL, BLU) 

C 

LS  *  BLL 

IF  (VAB1.LE.EPS)  GO  TO  81 
C 

BLU  =  BLL 
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DO  80  II  =  ORE,  N 
I  =  PERN (II) 

TLI  «=  HALF*TU  (I) 

FBI  «  BH  (I) 

QI  *  XONE-BHI 
CM  *  G  A  (I) 

DI  =  D (I) 

DEI  «  DE  (I) 

FIND  LU  SUFFICENTLY  L  ABG  E  SO  THAT  THE  I-TH  COMPONENT  OF  THE 
FUNCTION  P  IS  MCN-NEGATIVE  (HERE  TEST  IS  THE  I-TH  COHPONENT 
OF  THE  FUNCTION  P)  ,  THEN  SET  HLU  =  HAI (  LU  ) 

60  TEST  =  (  (XCNE-BHI)  •  (XONE-fDLOG  (  (XONE+TLI)  /  (XONE+ 

6  DEI*TLI )  )  /  (GAI*DI)  )  +(BHI/  (IONE+TLI)  ) -ALPHA) 

IF  (TEST. GT. ZEBO)  GO  TO  70 
TLI  =  HALMTLI 
CO  TO  60 

70  LUI  *  LAMBDA  (TLI  ,1) 

LU  (I )  =L0 1 

IF  (LUI. GT. HLU)  HLU  =  LUI 

80  CONTINUE 

SOLVE  P(LAHBDA)  =  0  FOB  LIHBD A  =  LS  IN  [HLL,HLU] 

HAI  *  H AXFN 
LB  =  HLL 
UB  =  HLU 

CALL  ZBBENT  (P  ,  E PS , NSI G, LB,U B ,H AX  ,  IER) 

LS  *  UB 

81  CONTINUE 

COMPUTE  THE  (S,S)  FOLICIES  FOR  THOSE  ITEHS  NITH  BIG  S  GREATE 
THAN  BSL 

82  DO  84  II  =  ONE,  N 
I  =  PERM  (II) 

G1I  =  GA  (I) 

TSI  =  TS  (I) 

MSI  =  HS  (I) 

DI  =  D (I) 

BIGS  (I)  =  MSI- {  (DLOG  (TSI)  ) /GAI) 

BECAUSE  OF  BOUND-OFF  ERROR,  IT  HAY  BE  THAT  BIG  S  LT  BSL. 

IF  SO,  SET  BIG  S  =  ESL. 

IF  (BIGS  (I)  .  LT.  BSL  (I)  )  BIGS(I)  =  BSL(I) 

84  LITS(I)  =  BIGS  ( 1 )  -  DI 
C 
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C  COMPOTE  THE  (S,S)  POLICIES  FOR  THOSE  ITBHS  11TH  BIG  S  EQUAL 
C  TO  BSL 
C 

IF  (*I.EQ.S)  GO  TO  B8 

85  B  *  If ONE 

DO  86  II  =  H,NI 
I  =  PERM (II) 

BIGS  (I)  *  BSL  (I) 

86  LITS(I)  *  BIGS  (I)-D  (I) 

C 

88  1TAVSL  =  ZEBO 
C 

C  COMPOTE  THE  S EB V I C E-L EVEL  AND  SHORTAGE  COST  FOB  ILL  THE  ITEMS 
C  II  THE  IBVENTOBY  SYSTEfl,  AND  ITAVSL  =  WEIGHTED  AVERAGE 
C  SEB VICE-LEVEL  (HCPEFULLY  THIS  EQUALS  IPPBOI  THE  ORIGINAL  ALPHA) 
C 

DO  90  II  =  GNE,N1 
1  =  PERM  (II) 

DBI  =  DE  (I) 

GAI  =  GA  (I) 

TSI  *  TS  (I) 

BHI  =  BH  ( I) 

QI  =  IONE-BHI 
HSI  =  BS  (I) 

DI  =  D (I) 

BIGSI  =  BIGS  (I) 

THS  =  DEIP  (-GAI*  (EIGSI-HSI)  ) 

ALI  *  QI*  (ICNEf DLOG  (  (XONEfTSI)/ (IONE+ 

6  DEI *1  SI )  )/(GAl*Dl))  +  (BHl/(XONE+TSI)  ) 

BSL(I)  =  ALI 

BTAVSL  =  VTAVSLf ALI*hT (I) 

90  CONTINUE 
C 

BETURI 

EID 


BEIL 

F0 

NOTION 

L 

An 

EDA* 

8  (THETA 

,1) 

COH3C 

N 

/L  Afl  / 

A, 

B, 

E,F, 

G  ,  U 

BEAL* 

8 

A  (512) 

(5 

12), 

E  (5  12)  ,F 

(512)  , G  (512) 

,H  (512) 

REAL* 

8 

THETA 

INTEG 

ER 

I 

LABBD 

A 

*  (A  (I 

)4 

TH 

ETA* 

(B(I)4TH 

ETA*  (E  (I)  +TH 

ET  A  *F  (I)  )  )  ) / 

&  (THE 

TA 

•  (G  (1) 

41 

HE 

T  A  *H 

(I)  )  ) 

BRTUR 

M 

END 
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C 


C 


BEIL  FU  NOTION  P*B  (LAN) 


COHNC*  /INFO!/  NT. ALPHA, IPS, NSIG, HAXFM, N 

COHMOM  /PABAHS/  TU , T S, BH, DE, G A, Df XONE, ZEBO, PEB H, ON E 

COMMON  /SOB 8/  XLAH,I 


B  E AL*8  NT  (512)  ,  ALPHA 
B E AL* 8  D  (512)  ,FH  (512)  ,GA  (512)  ,D 
BEAL*8  LAM, TSI , TUI, BHI,QI , LTHL, 


T  Lt  O  A  f  < 


c 

c 

c 


IBTBGEB  0NE,I,1I,HAX,N, PERM  (512) , HAXFN 
EXTEBNAL  LTHL 

XL AN  =  LAN 

P  =  XEBO 


C 

c 

c 

c 

c 


c 

c 

c 

c 

c 


TO  EVALUATE  P  (LAN)  ,  EVALUATE  THETA,  THE  SOLUTION  IN  [ 0 , TU ]  TO 
LTNL(THBTl)  =  LAHBE  A  (THETA)  -  LAN  *  0 
P IB  ST  CHECK  IF  1TNL  (TU)  =  0 


DO  20  II  =  ONE,  N 
I  =  PEB  N ( II ) 

TUI  *  TO  (I) 


SINCE  LTHL  (TUI)  IS  LESS  THAN  OR  EQUAL  TO 
LTNL(TUI)  IS  GBEATEB  THAN  OB  EQUAL  TO  0, 
1NDBED  LTHL  (TUI)  =  0 


0,  IF  IN  THE  COMPUTER 
THIS  IMPLIES  THAT 


TEST  =  LTML  (TUI) 

IP  (TEST.  LT.  ZEBO)  GO  TO  5 
TSI  =  TUI 
TS  (I )  =  TSI 
CO  TO  10 
5  LB  =  ZEBO 
OB  =  TUI 
MAX  =  MAXFN 

CALL  ZPALSE (LT ML , EPS , N SI G , LB, UB ,TSI , H AX, 1 EB) 

TS  (I )  =  TSI 
10  BHI  =  BH  (I) 

Q1  ='  XONE-RHI 

P=  P+KT  (I)  •  (QI*  (XCNE+DLOG(  (XONE+TSI)  /  (XONE+ 

C  DP  ( I)  *TS  I)  )  /  (GA  (I)  (I)  )  H  (R  HI /(XONE+TSI))  -  ALPHA) 

20  CONTINUE 


C 


C 


BETOBN 

END 
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C 

c 


c 


c 

c 


BBAL  FUNCTION  L?NL*8  (THETA) 

conbon  /Lin/  a,b,e,f,g,h 
COBNON  /SUBB/  LAN ,1 

§B AL*8  A(512)#B(512),E(512)#F(512)#G(512),H(512) 

B  B AL*8  THETA, LIN 

INTEGER  I 

LTHL  =  A  (I)  + 

&  THETA*  (fc  (I)  ”LAB*G  (I)  +THETA*  (E  (I) -L AN* H  (I)fTHET A* P  (I)  )) 

BBTUBH 
END 
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C 

C 

c 

c 

c 

c 


PBOuR AS  BISSNLSI 


c 

THE 

INPUT  IS  FBO 

c 

c 

TIT 

LE  =  USER  - 

C 

NI 

*  RUBBER  OF 

c 

K  = 

SET-UP  ORDE 

c 

L  = 

LE  A  DTI H  E  4 

c 

HO 

*  UNIT  HOLDI 

c 

HU 

=  SEAN  DEHAN 

c 

SD 

*  DEHAND  STA 

c 

WT 

*  CONSTBAINT 

c 

c 

NSIG  *  NURBER  0 

c 

ROUTINES 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


HULTI  -  ITEH  (S,S) 

BOB  NIL  DERAND  DISTRIBUTION  APPBOXIHATION 
LOW EB  BOUNDS  ON  BIG  S 
SINGLE-IT  EH  POLICY  OUTPUT 


SUPPLIED  TITLE  (LE  32  ALPHA  CHARACTERS) 

ITERS  (HI  LE  512) 

BING  COST 

1 

NG  COST 
D 

NEABE  DEVIATION 
WEIGHT 

F  SIGNIFICANT  DIGITS  DESIRED  FBOH  THE  IRSL 
ZBBENT  AND  ZFALSE 
EPS  =  1/(10  ••  BSIG) 

BAXFN  =  HAXIHUB  NURBER  OP  ITERATIONS  ALLOWABLE  BY  USER  FOB 
THE  IRSL  ROUTINES  ZBBENT  AND  ZFALSE 
IER  =  (OUTPUT)  ERROB  PARARETER  FBOH  IRSL  ROUTINES  AND  EIACT 

TBE  FOLLOWING  ABE  INPUT  TO  OB  CORPUTED  BY  GKD 

LIT S  =  LITTLE  S 
BIGS  *  BIG  S 

D  =  FEOBDEB  QUANTITY  BIG  S  -  LITTLE  S 
BSL  =  D  +  B*L*  HU 

*  LOWER  BOUND  FOB  BIG  S  (RAY  BE  INCREASED  BY  GKD) 

LS  *  L  AHBD A-STAB 

OUT  =  0  OTHERWISE 

1  WRITE  OUTPUT  TO  THE  EXTERNAL  FILE 
FORTRAN  UNIT  NURBER  10 

WITH  THE  FOLLOWING  DCB  SPECIFICATIONS: 


C  RECFH  =  FE 

C  LRECL  =  160 

C  BLKS1ZE  =  1600 

C 

C  THE  FOLLOWING  ARE  CCHPUTED  BY  EHRHARDT'S  POWER  RETHOD 
C 

C  PLUS  =  LITTLE  S 

C  PBIGS  =  BIG  S 

C 

C  THE  FOLLOWING  ARE  CORPUTED  EXACTLY 
C 

C  HC (I)  =  EXPECTED  HOLDING  COST  FOR  GKD  POLICY  FOR  ITER  I 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 


FC(I)  =  EXPECTED  REPLENISHMENT  COST  FOR  GKD  POLICY  FOB  ITEM  I 

SI  (I)  »  SEBYICE  LEVEL  FOB  GKD  POLICI  FOB  ITEH  I 

AVSL  =  VEIGHTE D-AVERAGE  SEBVICE  LEVEL  FOR  GKD  POLICIES 

PHC(I)  *  EXPECTED  HOLDING  COST  FOB  POVEB  POLICY  FOR  ITEH  I 

PEC  (I)  *=  EXPECTED  REPLENISHMENT  COST  FOR  POSER  POLICI  FOE  ITEH  I 

PSL(I)  *  SEBVICE  LEVEL  FOR  POVEB  POLICY  FOR  ITEH  I 

PAVSL  *=  WEIGHTED- AVERAGE  SEBVICE  LEVEL  FOB  POVEB  POLICIES 


INPUT  THE  SYSTEH 


P  A  F  AHETEBS  AS  FOLLOWS: 


TITLE 

LS  NS I G  HAIFN  NI  OUT 
WT(1)  D(1)  K  (  1)  L  (1 )  HO  ( 1) 

VT(2)  D(2)  K(2)  L  (2)  HO(2) 


HU  (  1)  SD  ( 1 )  B  (1) 
HU  ( 2)  SD  (2 )  B  (2) 


•  • 

VT (N)  D (N) 


K  (N) 


L  (N)  HO  ( N) 


•  •  • 

HU  ( N)  SD  (N)  B  (N) 


THE  INPUT  FOBHAT  IS: 


XXXXXXXIIXXXXXXXXXX 
IXIXIXXX. xxxxxx  xxxx 
XXX.  IXXX  XXX.  xxxx 
XXX. XXXX  XXX. xxxx 


xxxxxxxxxxxxx 

XX  XXXXXX  XXXXXX  X 
.  .  .  XXX. 

.  .  .  XXX. 


IXXI 

xxxx 


XXX. xxxx 
XXX.  xxxx 


xxx.xixx  XXX. xxxx 


XXX. 


XXXX 


XXX. 


XXXX 


ALL  DATA  STARTS  IN  CCLUHN  1  (EICEPT 
AND  THEBE  IS  A  SINGLE  SPACE  BETWEEN 


TITLE  STARTS  IN  COLUHN  2), 
EACH  DATA  ENTBY 


*•»»•«»«***»••»••*•»•***»»•*•*»******•*•*******•*•**************** 


PEAL* 8  LITS  (512)  ,  BIGS  (51 2)  , BSL  (51 2)  ,  D  (51 2)  ,  PLITS  (5 1 2)  ,PBIGS(512) 
BEAL* 3  K  (512)  #L(512)  ,HO(512)  ,  HU  (512)  ,SD(512)  ,  UT(512)  ,B 
BEAL* 8  AL#  HHISS, HPOW.TH1 SS,TPOW, EPS 

F  E AL*  3  KI,LI,HCI#HUI,SDI ,SLI,HCI,RCI#  AVSL,PAVSL,SUH 
PEAL*  8  SL(512),HC(512),RC  (512)  ,PSL  (512)  ,PHC  (512)  ,  PRC  (512) 

PEAL* 8  LSI,LITSI,BIGSI#D1 .TEST (51 2) 

INTEGER  SI 

INTEGER  I SI G, RAXFN, I TEH, TITLE (0) ,OUT 
BEAD  (1, 3)  (TITLE  (J)  ,J=1,8) 


non 
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3  F0RMAT(8I4) 

READ  ( 1,  5)  LSI, NSIG, HAXFN,N1, OUT 
5  FORMAT  (F15.  6,  II,  3  (16, IX)  ,11) 

EPS  =  10.**  (-NSIG) 

DO  10  I  «  1,NI 

R  BAD  (1,8)  HT(I),D(I),K(1),L(I),  HO  (1)  ,  HU  (I)  ,  SD  (I)  ,  B 
BSL(I)  =  D(I)4R*L  (I)  *H0(I) 

8  FORMAT (8 (F8. 4, IX)  ) 

10  CONTINUE 

n  AK  B  THE  HEIGHTS  HT  INTO  EQUIVALENT  CONVEX  HEIGHTS  HT 

sun  *  0. 

DO  15  I  =  1 , N I 
1  5  SOn  =  SUI-fWT  (I) 
sun  *  i.ODO/sun 
DO  16  I  =  1 , N I 
16  BT  (I)  =  SUH*W7  (I) 

C 

C  COIPOTB  THE  OPTIMAL  POLICIES  ASSOCIATED  WITH  THE  MULTIPLIER 
C  LA1BDA-STAB  (LSI) 

C 

CALL  BIGSLS (NI , L , HO, flU , S D, D, HT, BSL , LS I, E PS, NSIG, H AX FN , BIG S, L ITS ) 
C 

BBT7E(3,45) 

4  5  POPS  AT  ('I') 

C 

C  CCNPUTB  THE  (UNIFORM  SERVICE- LEVEL)  POWER  APPROXIBATION  POLICIES 
C  AND  COMPARE  IREIB  ASSOCIATED  COSTS  WITH  THE  COSTS  ASSOCIATED  WITH 
C  TtB  BIGSLS  POLICIES 
C 

CALL  COMP AR (NI , K , L, HO, MU ,SD, HT, BIGS, L ITS, TITLE, LSI ,OUT) 

STOP 

END 

C 

SUBROUTINE  BIGSLS  (NI , L , HO , MU , SD , D , WT , BSL, LSI , EPS, NSIG ,M AXPN, 

G  BIGS ,LI?S) 

C 

C  THIS  SUBROUTINE  COMPUTES  THE  OPTIMAL  (S,S)  POLICIES  ASSOCIATED 
C  WITH  THE  KUHN-TUCKER  CONSTRAINT  MULTIPLIER  LAMBDA-STAR  (LSI) 

C 
C 

COMMON 
C 

B  E  A  L*  8 
B  R AL*  3 
G 

R  E AL *  8 
G 


/FINDTS/AI,B1,EI,FI,GI,HI,LAH 
LITS  (512)  ,BIGS(512) 

L  (512)  ,HO  (512)  ,MU  (512)  ,SD  (512)  ,HT  (512)  , 

D  (512)  ,  BSL  (512) 

EPS,LSI,LI,U01,HUI,SDI,DI,MSI,WTI,TS1,LAH, 
RHI,QI,GAI,DEI, AI,BI,E1, FI,GI,HI,VAfa 1 , V AR2 , ET1 ,TUI , 
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G  KS1,LPTUI,BSLI,X,TUI,LB,UB,ZEF0,X0NE, HALF,  COHST 

BBAL*8  HC(512)  ,RC  (512)  ,SL  (512)  ,  PHC  (5 1  2)  ,  PBC  (5 12)  ,  PS  L  (51 2) 
BEAL*8  PLITS  (512)  ,PBIGS (512) 

C 

INTEGER  II, NSIG, SAXFN,HAX,ONE,CITa 
C 

B  B  A  L*  8  LANTH 
EXTERNAL  LA  STB 
C 

ZERO  =»  0. 

ONE  =  1 
XONS  =  1. 

HALF  ■=  .  5 

CONST  =  1. 700436904 
C 

DO  60  I  =  ONE, NI 
HOI  =  HO  (I) 

SJI  =  HD  (I) 

LI  =  L  { I ) 

SSI  =  LI*  HUI 
SDI  =  SD  (I) 

SiTI  =  NT  (I) 

DI  =  D (1) 

PHI  =  HU  1/ ( DI4  (  (  fl  U  I-f  (SDI*SDI/H0I)  )  *HALF)  ) 

QI  =  XONE-BHI 
G II  =  CONST/ (LI*SDI) 

DEI  =  DEXP  (GAI+DI) 

AI  =  DI*  HOI 

V  AB 1  =  XCNE+DEI 
BI  =  AI*  (QI-fVARI) 

El  =  AI*  (DEI4QI*V AR1) 

FI  =  AI*  QI*DEI 
X  AB 1  =  DI *B  HI*GAI 

V  AB  2  =  QI*  (DEI-XONE) 

G I  =  *TI*  (VAR14VAB2) 

HI  =  HTI*  (VAP1*DEI4VAE2) 

C 

C  DETERMINE  BSL  BT  THE  LARGEST  VALUE  OF  THETA,  TU,  AND  HENCE  THE 
C  SMALLEST  VALUE  OF  BIG  S,  SO  THAT  D  (L AH BDA) /D  (TH ET A)  LE  0, 

C  GUARANTEEING  THAT  LAH3DA  (THETA)  IS  INVERTIBLE  ON  (0,  TU) 

C  IF  THE  USER-SUPPLIED  BSL  IS  GREATER  THAN  THIS,  IT  BECOHES  THE 
C  BSL  USED 
C 

ETI  =  GI* (BI*HI*HI-GI*  (El *HI- FI*GI) ) -AI*HI*H1*HI 
C 

C  IF  FTA(I)  IS  LESS  THAN  0,  LAMBDA  (THETA)  HAY  NOT  BE  CONVEX,  AND 
C  THHPB  IS  NO  GUARANTEE  THAT  THE  GKD  ALGORITHM  WILL  CONVERGE  TO 

C  A  LOCAL  MINIBUS - THE  USER  SHOULD  EXERCISE  CARE 

C 

IP  (ETI.  LT.ZEBO)  1*FITE(3,30)  I, ETI 
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30  r3B11T(»  ET»(',I6,»)  *  1 ,  E20.  6 , 1  LT  0») 

C 

CNTS  =  -1 

35  C  NT  R  =  CNTB+CNE 

TUI  =  (DEXP  (GAI*  (HSI-CWTR) ) ) /DEI 
RSI*  (PI/HI)  -  (AI/  (GI*TOI*TUI)) 

YAB1  =  GI+HI*TUI 

LPTUI  *  KS1-  (ETI  /  (G1*HI*¥AR1*VAE1) ) 

IP  (LPTUI. GE. EEFC)  GO  TO  35 
C 

BSLI  =  CNTB+DI 
X  =  B5L  (I ) 

BSL(I)  =  ESLI 
IP  (I.GT.BSLI)  BSL(I)  =  I 
TUI  *  (DBXP(GAI*  (HSI-  (BSL  (I)-DI)  )  )  )  /DEI 
C 

LAB  =  LSI 
BAX  =  BAXPN 

LB  =  ZERO 
OB  =  TO I 
C 

C  IP  LIETH(OD)  IS  GREATER  TUAN  ZERO,  THEN  THE  FUNCTION 
C  LAMBDA  (THETA)  POF  THIS  IT  EB  IS  ALB  A IS  GREATER  THAN  LAHBDA-STAR, 
C  AND  SC  BIGS  POB  THIS  ITEM  IS  SET  AT  ITS  LOWER-DOUND  BSL 
C 

IP  (LAHTH  (UB)  .GT.ZEEO)  GO  TO  50 

CALL  ZBRENT (LAHTH, EPS, NS  I G, LB ,0B , H AX , IEE) 

TSI  =  UB 

BIGS  (I)  *  HSI- (  (DLOG  (TSI)  ) /GAI) 

C 

C  BECAOSE  OF  ROUND-OPF  ERROR,  IT  HAY  BE  THAT  BIG  S  LT  BSL. 

C  IF  SO,  SET  BIG  S  =  ESL. 

C 

IP  (BIGS  (I)  .  LT.  BSL  (I) )  BIGS  (I)  =  BSL(I) 

LITS(I)  *  BIGS  (I)  -DI 
GO  TO  60 

50  BIGS  (I)  =  BSL  (I) 

LITS(I)  =  BIGS  (I)-DI 
60  CONTINUE 
RETURN 
END 

BEAL  FUNCTION  L  A  HTH*8  (THETA) 

C 

COB  HON  /PINDT5/A1,BI,EI, PI, GI, HI, LAN 
BEAL* 8  A I , BI , El , FI , GI, HI , LAB, THETA 
C 

LAHTH  =  AI+ThEl  A*  (BI-LAH*GI(THET  A*  (El-LAH*bI  +  THETA*FI)  ) 

RETURN 

END 
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SUBROUTI ME  CCBPAB  (NI ,K , L, HO, HU, SD, WT, BIGS, LITS, TITLE, LSI, OUT) 
C 

C  THIS  SUBROUTINE  DOES  TBE  FOLLOWING: 

C 

C  (1)  EVALUATES  EXACTLY  THE  HOLDING  AND  B EPLEN ISHBENT  COST,  AND 
C  THE  SERVICE-LEVEL  FOB  THE  POLICIES  (LITS,BIGS)  UNDER  THE 

C  ASSUNPTICN  OF  N EGATIV E- B1NOBI A L  DEBAND  (IN  EXACT) 

C 

C  (2)  CONP’JTES  APFROXIBATELY  OPTIBAL  (S,S)  POLICIES  WITH  UNIPOBB 
C  SERVICE- LEVELS  (IN  POWAPP) 

C 

C  (3)  EVALUATES  THE  POLICIES  DESCRIBED  IN  (2)  EXACTLY  AS  IN  (1) 

C 

C  (4)  PRINTS  TBE  RESULIS  OF  STEPS  (1)  AND  (3) 

C 

1 EA L*8  LITS  (512)  ,EIGS(512) 

BEAL* 8  K  (512)  ,L  (512)  ,BO  (512)  ,BU  (512)  ,SD  (512)  ,1T  (512) 

BEAL* 3  HGKD, B POW , TGKD, TPOW , LS I, LIT5I,BIGSI 
BEAL* 8  KI,LI,HOI,BUl,SDI,SLI,HCI,RCI, AVSL,PAVSL 
BEIL* 8  HC  (512)  ,PC  (512)  ,SL  (512)  ,PLITS  (512)  ,PBIGS  (512)  , 

&  PHC(512)  ,  PFC  (512)  ,PSL  (512) 

C 

INTEGER  NI, TITLE  (8)  , OUT 
A  VS  L  *  0. 

R jKD  —  0  • 

TGKD  =  0. 

DO  70  I  =  1 , NI 

K1  =  K  (I ) 

LI  =  L  (I ) 

HJI  *  HO  (I) 

BUI  =  BU  (I) 

SDI  *  SD  (I) 

B1GSI  =  BIGS (I) 

LIT5I  =  LITS  (I) 

CALL  El  ACT ( KL , LI , HOI, B UI ,S DI , LIT SI , BI GS I , SL I , HCI , RCI , IEE ) 

IF  (IEB.  EQ.  1)  GO  TO  300 
A  VS  L  =  AVSL+SLI*WT  (I) 

BGKD  =  HGKD-fHCl 
TGKD  =  TGKD-fHCI+RCI 
HC (I)  -  HCI 
B  C  (I )  =  RCI 
S  L  (I )  =  SLI 
70  CONTINUE 
C 

PAVSL  *  0. 

UPON  =  0. 

TPOW  =  0. 

DO  75  I  =  1, NI 
K I  =  K (I) 

LI  =  L  (I ) 


non 
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boi  *  no  (i) 

BUI  «  BU  (I) 

SDI  *  SD  (I) 

CUL  PON APP  (KI  ,LI ,  HC'l,  BUI,SDI,AVSL,LITSI,BIGSI) 

Cm  EX  ACT  (KI,II,HOI,HOI,SDI,LITSI,BIGSI  , SLI ,HCI  ,  BCI,IEB) 

IF  (IER.  IQ.  1)  GO  TO  300 
PAKSL  *  PAVSL4SLI*UT  (I) 

BPOK  *  HPON+HCI 
TPON  =  TPOU4BCI4 BC1 
PLUS  (I)  •  LI15I 
PBIGS(I)  =  BIGSI 
PHC(I)  «  BCI 
PRC  (I)  =  BCI 

PSL(I)  =  SLI 
75  COITINUE 

PRIBT  SUHHAB1 

SRITE  (3,98)  (TITLE  (J)  ,  J=  1 , 8)  ,  LSI 
98  FORHAT(T2,8A4,  '  LA3BDA-STAP.  =  ',Fl5.6//> 

B RITE (3, 100)  A VSL, PAVSL,HGKD,HPOB,TGKD, TP CK 
100  FOE3AT(T3,' ACTUAL  ALPHA  (GKD)  = • ,T27, FI 1 . 4, //, 

G  T3  , '  ACT  UAL  ALPHA  (PON)  = • ,T27, PI  1 . 4,//, 

G  T3,'RCLDING  COST  (GKO)  =  *  ,T27,  P 1 1 .  4,  //, 

G  T3,  ' HCIDI NG  COST  (POi)  = • ,T27 , P 1 1 . 4 ,//, 

&  T3, 'TOTAL  COST  (GKD)  =  '  ,T27 , F 1 1 . 4, //, 

6  T3,'TCTAL  COST  (PON)  = ' ,T27,F1 1 .4,////, 

G  T4,'I',T8,'K* ,T12,'Li1',T19,'H',T25, 'HU', 

G  T32, 'SE',T40,'GKD' ,T50,« PONEB • ,T65,' HOLD' , 

6  T85, 'EEP' ,T99,'SER  LE? • # T1 16 , • N*NT ' ,//, 

&  T  39, 'LS',T44,'BS',750,'LS',T55,'BS', 

&  T62, »GKD» ,T70,' POKER' ,T80, 'GKD • ,T88, ' POKES • , 

G  T98,'GRC',T105, 'PONEB' ,//) 

DO  130  I  =  1,« 

SUB  =  NI*NT  (I) 

K  RITE  (3,  120)  I,K  (I)  ,  L  (I)  ,  RO(I),HU(I)  ,  SD  (I)  ,  LITS  (I)  ,  BIGS  (I)  , 

C  PLUS  (I)  ,PBIGS(I)  ,HC  (I)  ,PHC(I)  ,BC  (I)  ,PBC  (I)  , 

8  SL(I),PSL(I)  ,  SUH 

KRITE  (10,  120)  I,K  (I)  ,L{I)  ,HO(I)  ,HU  (I)  ,SD  (I)  ,  LITS  (I)  ,  BIGS  (I)  , 
t  PLUS  (I)  ,PBIGS(I)  ,UC  (I)  ,PHC  (I)  ,BC  (I)  ,PBC(I), 

G  SL  (I)  ,PSL  (I)  ,SUH 

120  rJBMAT(I4,P5.0,P4.0,P9.4,P4.0,P9.4,2(F5.0,F6.0)  ,4F9.2, 

&  2F8.  4, FI 2. 4,//) 

130  CONTINUE 
RETURN 

300  KRITE  (3,310)  I 

310  FOB.1AT  (•  BIG  S  OR  D  >  2000  FOR  ITEB',16,/ 

G  '  EXACT  CANNOT  EVALUATE  THE  (S,S)  POLICY',//) 
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C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


RETURN 

END 

SUBROUTINE  EXACT (K, L, H , H U , SD, LS , BS #SL, RC , BC# IEB) 

THIS  SUBBOOT I  ME  COHFUTES  THE  SERVICE-LEVEL,  BOLDING  COST  AND 
FEPLENISHHENT  CCST  EXACTLI  FOP  THE  POL1CT  (LS,BS)  UNDER  THE 
ASSUHPTION  OF  NEGATIVE  BINOHIAi  DEHAND. 

K  *  SET-UP  COST 

L  «=  LEADTIHE  4  1 

H  «=  UNIT  HOLDING  COST 

HO  =  HEIN  DEHANC 

SD  =  DESAND  STAND  ABC  DEVIATION 

LS  =  LITTLE  S 

BS  =  BIG  S 

SL  =  ST  BA  DI-ST AT  E  SERVICE-LEVEL  (FREQUENCY  OF  PERIODS  THAT 
NO  BACKORDER  IS  PLACED) 

BC  =  EXPECTED  HOLDING  COST 

EC  =  EXPECTED  F EPLEN1SBHENT  COST 

E EAL* 8  K,L,H,HU,SD,LS,BS,SL,HC,  RC,P,Q,B,BS,PROB(2000)  ,HALF, 

&  NB  (200  0)  ,  S  (2000)  ,  HI ,  RHO,  V  AR ,S0  H  ,  PROBI ,  HD,  Z EROf  ZONE,  IT  NO 

INTEGER  S,D,tB1,  UB2,ONE,TWO,IER 

DATA  ONE,TWO,ZERO,ICNE,XTWO,HALF/1,2,O.DO,1.DO,2.DO,.  5D0/ 

S  AND  D  ABE  ROUNDED  TO  THE  NEAREST  INTEGER.  NE  USE  THE  REORDER 
QUANTITY  D  =  BIG  S  -  LITTLE  S  -  1  BECAUSE  WE  USE  THE  VERSION  OF 
AN  (S, S)  POLICY  THAT  BEQUIRES  THAT  AN  ORDER  BE  PLACED  WHEN  THE 
INVENTORY  POSITION  DROPS  STRICTLY  BELOW  LITTLE  S. 

PR03  (I)  =  PROB  (NB  PV  =  1-1) 

N B (I)  =  PBOB  (NB  4  •  •  •  4  NB  RV  =  1-1) 

I _ I 

L 

B  (I)  =  PBOB  (RENEWAL  FNC  OF  NB  RV  =  1-1) 

I  EE  =  0  EXACT  WAS  ABLE  TO  EVALUATE  (S,S)  POLICY 
1  OTHERWISE 

IER  =  0 
S  =  BS+HALF 
D  *  (BS-LS-ONE) +HALF 
Q  =  H'J/(SD*SD) 

P  =  XONE-C 

R  =  H  U*  Q/P 

BS  =  L*  B 

VAR  =  Q**P 

PROB  (ONE)  =  VAB 

NB  (ONE)  =  Q**RS 


nnn  nnnnnn 
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BHO  *  IONE/  (XONE-VAR) 

II  (ONE)  *  ?AB*BHO 
OBI  =  S+ONE 

IP  (D.GT.S)  OBI  =  D+OHE 
IF  (OB  1.GT.  2000)  GO  TO  100 
DO  20  I  =  TWO, OBI 
VAR  =  ZONE/  (I-XCNE) 

PBOBI  =  PBOB  (I-ONE)  *P*  (B  +  I-ITBO)  ♦VAR 
PBOB(I)  =  PBCBI 

N 3  (I)  =  NB(I-CNE)  *P*  (RS+I-XTWO)  *VAR 
SUM  =  ZEBO 
0B2  =  I-ONE 
DO  10  J  =  0NE,0B2 
10  SOU  =  SOH+PBOB  (I-J+CNE)  *H  (J) 

H(I)  =  (PBOBI+SUn) ♦BHO 
20  CONTINUE 

0B1  =  D+OHE 

IP  (OBI,  Gl.  2000)  GO  TO  100 
no  =  ZERO 
DO  30  I  =  ONE, 0B1 
30  no  =  no+n  (i) 

BhO  =  ZONE/ (XCNE4HD) 

COHPUTE  EXPECTED  R  EPLENISHHENT  COST 

RC  =  K*RHO 

COMPUTE  EXPECTEE  HOLDING  COST 

IP  (D.GT.S)  UE1  =  S-tONE 

IP  (UB1.GI.  2000)  GO  TO  100 

HC  =  ZEBO 

DO  40  I  =  ONE, OBI 

OB2  =  S-I+TWO 

HI  =  «{I) 

DO  40  J  «  C NE , OB2 
HC  =  HC+  (UB2-J)  *NB(J)  ♦HI 
40  CONTINUE 

SUM  =  ZERO 
UB2  =  S+ONE 
DO  45  I  =  ONE ,OB2 
4 S  sun  =  sun-f  (032-1)  ♦NB(I) 

HC  =  H*  (HC+SUN)  ♦RHO 

COMPUTE  (STEADT-STATE)  SERVICE-LEVEL 

SL  =  ZERO 
DO  60  I  =  ONE, OBI 
sun  =  ZERO 
UB2  =  S-  I-f  TWO 


n  n 
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DO  50  J  *  CNE,UB2 
50  SUH  =  SUHiNB(J) 

60  5L  *  SL4SUH*H  (I) 
sun  =  ZEBO 
OBI  =  S-fONE 

IF  (UB1.GI.  2000)  GO  TO  100 
DO  70  I  *  ONE ,UB 1 
70  son  =  sun+NB(i) 

SL  =  (SL4SUH)  *RHO 
RETURN 
C 

100  IER  =  1 
RETURN 
C 

END 


SUBROUTINE  POiAPF  (K,L, H,  HU , SD, ALPHA, LITS , BIGS) 

C 

C  THIS  SDBPOOTINI  CONFUTES  LITTLE  S  AND  BIG  S  USING  THE  REVISED 

C  POWER  A  PPROXI HATIC  N  HETHOD  (HOSI ER  { 1 98 1 ) ,  TECH  BEPORT  #18). 

C  THE  SHORTAGE  CCS!  IS  COHPUTED  USING  THE  POWER  APPROXIHATION 

C  HETHOD  (EHBHABDT  (1977,  PP  18,  4 5),  TECH  REPORT  #12). 

C 

C  K  =  SET-UP  COST 
C  L  *=  LE  A  DTI  HE  4  1 
C  H  =  UNIT  HOLDING  COST 
C  HU  =  HEAN  DEHAND 

C  SD  =  HEAN  STANDARD  DEVIATION 

C  ALPHA  =  (STEADY-STATE)  SERVICE-LEVEL,  THE  FREQUENCY  OP  PERIODS 
C  WITHOUT  BACKLOGS 

C  LITS  =  LITTLE  S 

C  BIGS  =  BIG  S 

C 

REAL* 8  P,H, K, HU, D, L , VL , S DL, Z , V, SD, LITS , BIGS , ONE, 

G  C1,C2,C3,C4,C5,C6,C7,C8,C9, CIO, ALPHA 

C 


DATA 

Cl/1 

.3/ 

DATA 

C2/. 

494/ 

DATA 

C3/» 

506/ 

DATA 

C4/. 

116/ 

DATA 

C5/. 

973/ 

DATA 

C6/. 

183/ 

DATA 

C7/1 

.  063/ 

DATA 

C8/2 

.  192/ 

DATA 

CIO/ 

.0695/ 

DATA 

ONE/ 

1.  / 

V  = 

SD*  S 

D 

VL  = 

¥*L 

SDL  *  DSQBT (VL) 


P  =  H*(UPHl-C10)/(CNE-lLPHi) 

Z  =  OME-f  (VL/(RU*flO)  ) 

D  =  Cl*  (HU**C2)  *  (  (K/H)  **C3)  •  (Z**C4) 
Z  =  DSQBT  (D*B/ (P*SDL) ) 

C9  =  (C6/Z)-fC7-C8*Z 
LITS  »  C5*L*8D-fSDL*C9 
BIGS  =  LITS+D 
RBTURN 
END 
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